379 double a222 =
mos_exec(model,here,ms,0.0,0.0,0.0);
381 double a322 =
mos_exec(model,here,ms,delta,0.0,0.0);
382 double a232 =
mos_exec(model,here,ms,0.0,delta,0.0);
383 double a223 =
mos_exec(model,here,ms,0.0,0.0,delta);
385 double a122 =
mos_exec(model,here,ms,-delta,0.0,0.0);
386 double a212 =
mos_exec(model,here,ms,0.0,-delta,0.0);
387 double a221 =
mos_exec(model,here,ms,0.0,0.0,-delta);
389 double a233 =
mos_exec(model,here,ms,0.0,delta,delta);
390 double a323 =
mos_exec(model,here,ms,delta,0.0,delta);
391 double a332 =
mos_exec(model,here,ms,delta,delta,0.0);
393 double a333 =
mos_exec(model,here,ms,delta,delta,delta);
395 double a123 =
mos_exec(model,here,ms,-delta,0.0,delta);
396 double a132 =
mos_exec(model,here,ms,-delta,delta,0.0);
397 double a213 =
mos_exec(model,here,ms,0.0,-delta,delta);
398 double a312 =
mos_exec(model,here,ms,delta,-delta,0.0);
399 double a321 =
mos_exec(model,here,ms,delta,0.0,-delta);
400 double a231 =
mos_exec(model,here,ms,0.0,delta,-delta);
403 double a022 =
mos_exec(model,here,ms,-(delta+delta),0.0,0.0);
404 double a202 =
mos_exec(model,here,ms,0.0,-(delta+delta),0.0);
405 double a220 =
mos_exec(model,here,ms,0.0,0.0,-(delta+delta));
407 cd->
d2_p2 = ((a322 - a222)*dd - (a222 - a122)*dd)*dd;
408 cd->
d2_q2 = ((a232 - a222)*dd - (a222 - a212)*dd)*dd;
409 cd->
d2_r2 = ((a223 - a222)*dd - (a222 - a221)*dd)*dd;
410 cd->
d2_pq = ((a332 - a322)*dd - (a232 - a222)*dd)*dd;
411 cd->
d2_qr = ((a233 - a232)*dd - (a223 - a222)*dd)*dd;
412 cd->
d2_pr = ((a323 - a322)*dd - (a223 - a222)*dd)*dd;
414 cd->
d3_p3 = (((a322 - a222)*dd - (a222 - a122)*dd)*dd -
415 ((a222 - a122)*dd - (a122 - a022)*dd)*dd)*dd;
416 cd->
d3_q3 = (((a232 - a222)*dd - (a222 - a212)*dd)*dd -
417 ((a222 - a212)*dd - (a212 - a202)*dd)*dd)*dd;
418 cd->
d3_r3 = (((a223 - a222)*dd - (a222 - a221)*dd)*dd -
419 ((a222 - a221)*dd - (a221 - a220)*dd)*dd)*dd;
421 cd->
d3_p2r = (((a323 - a223)*dd - (a223 - a123)*dd)*dd -
422 ((a322 - a222)*dd - (a222 - a122)*dd)*dd)*dd;
424 cd->
d3_p2q = (((a332 - a232)*dd - (a232 - a132)*dd)*dd -
425 ((a322 - a222)*dd - (a222 - a122)*dd)*dd)*dd;
427 cd->
d3_q2r = (((a233 - a223)*dd - (a223 - a213)*dd)*dd -
428 ((a232 - a222)*dd - (a222 - a212)*dd)*dd)*dd;
430 cd->
d3_pq2 = (((a332 - a232)*dd - (a322 - a222)*dd)*dd -
431 ((a322 - a222)*dd - (a312 - a212)*dd)*dd)*dd;
433 cd->
d3_pr2 = (((a323 - a223)*dd - (a322 - a222)*dd)*dd -
434 ((a322 - a222)*dd - (a321 - a221)*dd)*dd)*dd;
436 cd->
d3_qr2 = (((a233 - a223)*dd - (a232 - a222)*dd)*dd -
437 ((a232 - a222)*dd - (a231 - a221)*dd)*dd)*dd;
439 cd->
d3_pqr = (((a333 - a233)*dd - (a323 - a223)*dd)*dd -
440 ((a332 - a232)*dd - (a322 - a222)*dd)*dd)*dd;