Table of Contents
Dynamic implicit/explicit integration scheme
Description
Can shift from implicit scheme to explicit scheme and reciprocally.
Data set
Old Metafor Version <= 2422
See dynamic implicit and explicit for a description of schemes and the definition of density and initial velocities.
In addition, parameters must be defined to handle the shift (Global Parameters [REMOVED]):
MDE_SHIM
: implicit scheme (=MDE_NDYN
)MDE_SHEX
: explicit scheme
Shifting can be done automatically or manually or both. Another application of this algorithm is to stabilize explicit steps.
New Metafor Version > 2422
ti = IeTimeIntegration(metafor) ti.setExplicitTimeIntegration(tiExp) ti.setImplicitTimeIntegration(tiImp) ti.setCurrentDynamicTimeIntegration(tiImp)
where tiExp and tiImp are the references towards the explicit and implicit time integration scheme respectively.
Automatic Shift (method 1 or 2)
By default, the algorithm chooses the method based on an integration error, the time step and the machine time. The definition MDE_ICPU = 1
must be done (shift based on CPU time, otherwise dummy shift for the battery). The shift is done only if estimated and critical time steps are equal (to a safety factor). When shifting from an explicit alphageneralized method, explicit time steps can be computed with a bifurcation radius equal to zero to increase numerical dissipation. The option is activated for method=1
but not for method=2
(see Implicit/explicit shift management).
Manual Shift (method 3 or 4)
The user can impose shifts (see command .sie
). When shifting from an explicit alphageneralized method, explicit time steps can be computed with a bifurcation radius equal to zero to increase numerical dissipation. The option is activated for method=3
but not for method=4
(see Implicit/explicit shift management).
Balancing
When shifting from explicit to implicit, the first implicit time step is computed based on an integration error (majored by a value equal to the maximal change of step times the explicit step) in automatic shift and is equal to this value (maximal change of step times the explicit step) in manual shift. In addition, the user can impose a number of implicit balancing time steps when shifting from explicit to implicit. An implicit balancing step is computed based on an explicit number of steps corresponding to the maximal change, justifying why this is an integer.
Balancing all explicit time steps (method 5)
All the computation is done with explicit time steps (method5
). All explicit maximal change of explicit steps are balanced by an implicit step.
Implicit/explicit shift management
It is mandatory to instanciate the IeShifManager and then the parameters can adjusted if required :
iesm = IeShiftManager(metafor) iems.setIeShiftMethod(meth) iesm.setNbOfBalanceSteps(nbbs) iesm.setMaxStepIncrease(msin) iesm.setAutomaticShiftSecurity(secu) iesm.setCpuDependency(cpuDep)
^ Parameter ^ default ^ description ^
meth  = IESM_NOSHIFTMETHOD : no shift (default) 

= IESM_AUTOMATICSHIFTMETHODWITHDISS shift done automatically from integration error and CPU times (requires MDE_ICPU=1 ); when shifting to an implicit method, explicit time steps are computed with numerical dissipation 

= IESM_AUTOMATICSHIFTMETHOD : shift done automatically from integration error and CPU times (requires MDE_ICPU=1 ); when shifting to an implicit method, there is no explicit time steps computed with numerical dissipation 

= IESM_MANUALSHIFTMETHODWITHDISS : shift times are set with the command .sie (see below); when shifting to an implicit method, explicit time steps are computed with numerical dissipation 

= IESM_MANUALSHIFTMETHOD : shift times are set with the command .sie (see below); when shifting to an implicit method, there is no explicit time steps computed with numerical dissipation 

= IESM_BALANCEDEXPLICITMETHOD : every msin explicit steps, balanced is done with an implicit step to stabilize nonlinear computing 

nbbs  0  when shifting to implicit, is equal to the number of balancing steps 
msin  6  maximal change of the time step; this value is the ratio of the implicit step to the explicit step for a change to implicit or a balancing, so this is an integer 
secu  1.5  safety coefficient for automatic shift 
cpuDep  False  specify if the automatic criterion for shifting from one time integration to another depends on the User CPU time. (Metafor Version > 2422) 
For manual shift, the user can set times at which method is changed. This is done with the commands:
metafor.getIeShiftManager().pushTimeOfShiftingToExplicit(time) metafor.getIeShiftManager().pushTimeOfShiftingToImplicit(time)
where
time  time at which the shift takes place 