Diss.-ETH 12317 Diss.-ETH 12317 Diss.-ETH 12317 Diss.-E Eine ...

5.5. ADAPTION DES DIENSTES NETZBERECHNUNG 105
anz_vars = 8 ;
anz_pars = 10 ;
anz_funs = 4 ;
state
break;
= PUQMode ;
case
...
PUQMode :
default:
break;
}
}
Die folgenden Zeilen zeigen die erzeugte Methode GetAElem() des
PU-Generators. Diese zentrale Methode erzeugt für alle spezifizierten
Zustände je einen Zugriff auf die schwach besetzte Jacobi-Matrix. Er
wird durch geschachtelte Abfrage mittels case-Anweisungen implementiert.
Der Default-Fall liefert den Wert 0 zurück. Für die Nichtnullen der
Matrix setzt die Adaption das aus den Spezifikationen von Vars, Funs und
StateDefs symbolisch berechnete Jacobi-Element ein.
double PUGen :: GetAElem(int i, int j, int state)
{
double t0;
switch(state) {
case PUMode :
switch(i) {
case 1 :
switch (j) {
case 1 :
t0 = 1.0 ;
return t0;
break;
case 3 :
t0= -P1/(e1*e1+f1*f1)
+ 2.0*(e1*P1+f1*Q1)/
pow(e1*e1+f1*f1,2.0)*e1-G ;
return t0;
break;
Diss.-ETH 12317 Diss.-ETH 12317 Diss.-ETH 12317 Diss.-E