pwSetIntegrationOptions

 

 oldValue = pwSetIntegrationOptions(optionName, newValue)

 A set of options can be set for the numerical integration.
 Each integrator has its own set of parameters, often these
 sets overlap, i.e. most integrators have the option 'MaxStep'.
 Within this function you can set the option values. If you
 switch to a different integrator where the option is not
 available, the specified value will be ignored.
 If the oldValue is empty, than the default value has been used before.

 Possible values (with default in brackets) for the integrators:

Arguments for pwSetIntegrationOptions

 RADAU5 (the default integrator):

  RelTol (1e-3):               relative tolerance
  AbsTol (1e-6):               absolute tolerance
  eps (1e-16):                 rounding unit
  rho (0.9):                   safety factors for step size control
  StepSizeMinSelection (0.2)
  StepSizeMaxSelection (8.0)
                               parameters for step size selection
                               StipSizeMinSelection <= tauNew/tauOld <= StepSizeMaxSelection
  MaxStep (tEnd-tStart):       maximal step size
  InitialStep (1e-6):          initial step size guess
  MaxNumberOfSteps (100000):   maximal number of allowed steps
  FreezeStepSizeLeftBound (1.0)
  FreezeStepSizeRightBound (1.2)
      If FreezeStepSizeLeftBound<=tauNeu/tauOld<=FreeStepSizeRightBound,
      then the step size is not changed. This saves LU-decompositions
  StepSizeStrategy (1)
      switch for step size strategy; possible values are 1 and 2
      1: mod. predictive controller (Gustaffson)
      2: classical step size control
  MaxNewtonIterations (7)
      the maximal number of Newton iterations for the solution of
      the implicit system in each step.
  StartNewtonWithZeros (0):    flag for starting value for Newton's method.
       =0: extrapolated collocation solution is taken as starting value
      ~=0: zero starting values are used
  NewtonStopCriterion (min(0.03,sqrt(RelTol(1)))): stopping criterion for Newtons's method
  RecomputeJACFactor (0.001)
      decides whether the jacobian should be recomputed.
      The higher this value (like 0.1) the less often
      the jacobian is evaluated. If evaluation of the jacobian
      is cheap, take a smaller value.
      If RecomputeJACFactor<0 then the Jacobian is evaluated
      after every accepted step.


 Other integrators: Under construction (contact support potterswheel.de)

Example

 
 oldValue = pwSetIntegrationOptions('RelTo', 1e-5);

See also

odeset