skade eaction ffusion ne-dimensional System - ZIB

More documents

Recommendations

Info

verboseP scaling masteps nt tstart eal tend eal timestep eal globtol spacetol eal timetol eal itertol eal itermasteps nt showsol nt atoli eal rtoli eal Important da prnted. A he momnt, always true. ep size control: 0 absolute, 1 for relative. If scaling is set to be 1, so all corresponding parameters atoli, rtoli il,2,...,number of unknowns have to be et too. umber of maximal tim teps starting tim final time initial ime ep accuracy demanded for the integrator (internal pacetol and timetol ar et automatically) desred space accurac desired time accuracy (if spacetol and timetol are selected, hen internal globtol is set) desred relative accuracy of he ierativ olver for the linear sysems standard is 0.0001 maximal number of iterations, standard is 1000 graphic is plotted every Int seps, standard is 1 i=l,...,number of unknowns, necessary input for every unknown to control the time step size relatively (scaling 1). The tolerance atol should ndicate the absolute value at which the corresponding component is essntiall insignificant. The tolerance rtol affects the relative accurac with respect to the maximal value in time.
haper r o g e r only This chapter gives some information to get a feeling for the internal structur of the program Of course, it cannot mention all details. Furthermore for a better understanding it is necessary o read al header and c-files of K ARD OS very carefully. 5.1 Defining a time problem To include all data of a initial boundary value problem two new data struc ures are added. struct timeProblem char *name; int (*Parabolic)(...); int (*Laplace)(...); int (*onvection)(...); int (*elmholtz)(...); int (*ource)(...); int (*Jacobian)(...); real (*InitialFunc)(...); int (*auchy)(...); int (*irichlet)(...); int (*ol)(...); struct timeProblemType char *name; int maxNoOfTimeProblems; int noOfEquations; The access to these srucures is realized by he pointers timeProblem *actTimeProblem, *timeProblems; timeProblemType *theTimeProblem; Furthermore, there are some functions t announc ertain problms within he program int etTimeProblem(char *name, char **varames int *Parabolic, 10
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: haper C o n d ng All specal command
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 and 22: haper le In this chapte how how to
Page 23 and 24: eal valO, vail; real *fVal, int equ
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?