12.07.2015 Views

GNU Octave - Local Sector 7 web page

GNU Octave - Local Sector 7 web page

GNU Octave - Local Sector 7 web page

SHOW MORE
SHOW LESS
  • No tags were found...

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

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 23: Differential Equations 193"compute consistent initial condition"Denoting the differential variables in the state vector by ‘Y_d’ and thealgebraic variables by ‘Y_a’, ddaspk can solve one of two initializationproblems:1. Given Y d, calculate Y a and Y’ d2. Given Y’, calculate Y.In either case, initial values for the given components are input, and initialguesses for the unknown components must also be provided as input. Setthis option to 1 to solve the first problem, or 2 to solve the second (thedefault default is 0, so you must provide a set of initial conditions thatare consistent).If this option is set to a nonzero value, you must also set the "algebraicvariables" option to declare which variables in the problem are algebraic."use initial condition heuristics"Set to a nonzero value to use the initial condition heuristics options describedbelow."initial condition heuristics"A vector of the following parameters that can be used to control the initialcondition calculation.MXNIT Maximum number of Newton iterations (default is 5).MXNJ Maximum number of Jacobian evaluations (default is 6).MXNH Maximum number of values of the artificial stepsizeparameter to be tried if the "compute consistent initialcondition" option has been set to 1 (default is 5).Note that the maximum number of Newton iterationsallowed in all is MXNIT*MXNJ*MXNH if the "computeconsistent initial condition" option has been set to 1and MXNIT*MXNJ if it is set to 2.LSOFFSTPTOLEPINITSet to a nonzero value to disable the linesearch algorithm(default is 0).Minimum scaled step in linesearch algorithm (default iseps^(2/3)).Swing factor in the Newton iteration convergence test. Thetest is applied to the residual vector, premultiplied by theapproximate Jacobian. For convergence, the weighted RMSnorm of this vector (scaled by the error weights) must be lessthan EPINIT*EPCON, where EPCON = 0.33 is the analogoustest constant used in the time steps. The default is EPINIT= 0.01."print initial condition info"Set this option to a nonzero value to display detailed information aboutthe initial condition calculation (default is 0).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!