skade eaction ffusion ne-dimensional System - ZIB

More documents

Recommendations

Info

This system of nonlinear equations modls a cerain planktonit predator prey suations in which crowding is a facor. The KARDOS user defines in user. he following funtions: static int EcoParabolic(real x, int classA, real t, real *fVals, int equation, int variable) if (equation!=variable) return F_IGN0RE; switch (variable) { case 0: fVals[0] = 1.0; break; case 1: fVals[0] = 1.0; break; return F_C0ISTAIT; static int EcoLaplace(real x, int classA, real t, real *fVals, int equation, int variable) if (equation!=variable) return F_IGN0RE; switch (variable) { case 0: fVals[0] = 0.0125; break; case 1: fVals[0] = 1.0; break; return F_C0ISTAIT; static int EcoConvection(real x, int classA, real t, real *fVal, int equation, int variable) return F_IGI0RE; static int EcoHelmholtz(real x, int classA, real t, real *fVal, int equation, int variable) return F_IGI0RE; static int EcoSource(real x, int classA, real t, EDG* ed,
eal valO, vail; real *fVal, int equation, real (*Func)(real,EDG*,int,timeIntegMethod*)) switch (equation) { case 0: valO = Func(x,ed,0,actTimeIntegtor); vail = Func(x,ed,l,actTimeIntegtor); fVal[0] = ((35.0+16.0*val0-val0*val0)/9.0-vall)*val0; break; case 1: valO = Func(x,ed,0,actTimeIntegtor); vail = Func(x,ed,l,actTimeIntegtor); fVal[0] = (val0-1.0-0.4*vall)*vall; break; return F_VARIABLE; static int EcoJacobian(real x, int classA, real t, EDG* ed, real *fVal, int equation, int variable, real (*Func)(real,EDG*,int,timeIntegMethod*)) real valO, vail; switch (equation) { case 0: if (variable==0) { valO = Func(x,ed,0,actTimeIntegtor); vail = Func(x,ed,l,actTimeIntegtor); fVal[0] = (35.0+32.0*val0-3.0*val0*val0)/9.0-vall; else if (variable==l) { vail = Func(x,ed,l,actTimeIntegtor); fVal[0] = -vail; else ZIBStdOut("error in EcoJacobian\n"); break; case 1: if (variable==0) { vail = Func(x,ed,l,actTimeIntegtor); fVal[0] = vail;
Page 1 and 2: Tehnicl Repot TR 9-9 (July 193) Jen
Page 3 and 4: Contents Introducti efng a problem
Page 5 and 6: Yet this frm not suited to be imple
Page 7 and 8: To cntrol the pace dicretiztion, we
Page 9 and 10: haper efining roblem At first the u
Page 11 and 12: haper C o n d ng All specal command
Page 13 and 14: haper r o g e r only This chapter g
Page 15 and 16: int nonNegativity ontrol of the non
Page 17 and 18: eal spaceTolFac paceTolglobTol, san
Page 19 and 20: espping
Page 21: haper le In this chapte how how to
Page 25 and 26: fVals[l] = 0.0; /* */ else ZIBStdOu
Page 27 and 28: seldraw 1 timesteppin timestep 0.00

skade eaction ffusion ne-dimensional System - ZIB

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?