public class ExpandableStatefulODE extends Object
One typical use case is the computation of the Jacobian matrix for some ODE. In this case, the primary set of equations corresponds to the raw ODE, and we add to this set another bunch of secondary equations which represent the Jacobian matrix of the primary set.
We want the integrator to use only the primary set to estimate the
errors and hence the step sizes. It should not use the secondary
equations in this computation. The integrator
will
be able to know where the primary set ends and so where the secondary sets begin.
FirstOrderDifferentialEquations
,
JacobianMatrices
Constructor and Description |
---|
ExpandableStatefulODE(FirstOrderDifferentialEquations primary)
Build an expandable set from its primary ODE set.
|
Modifier and Type | Method and Description |
---|---|
int |
addSecondaryEquations(SecondaryEquations secondary)
Add a set of secondary equations to be integrated along with the primary set.
|
void |
computeDerivatives(double t,
double[] y,
double[] yDot)
Get the current time derivative of the complete state vector.
|
double[] |
getCompleteState()
Get the complete current state.
|
FirstOrderDifferentialEquations |
getPrimary()
Get the primary set of differential equations.
|
EquationsMapper |
getPrimaryMapper()
Get an equations mapper for the primary equations set.
|
double[] |
getPrimaryState()
Get primary part of the current state.
|
double[] |
getPrimaryStateDot()
Get primary part of the current state derivative.
|
EquationsMapper[] |
getSecondaryMappers()
Get the equations mappers for the secondary equations sets.
|
double[] |
getSecondaryState(int index)
Get secondary part of the current state.
|
double[] |
getSecondaryStateDot(int index)
Get secondary part of the current state derivative.
|
double |
getTime()
Get current time.
|
int |
getTotalDimension()
Return the dimension of the complete set of equations.
|
void |
setCompleteState(double[] completeState)
Set the complete current state.
|
void |
setPrimaryState(double[] primaryState)
Set primary part of the current state.
|
void |
setSecondaryState(int index,
double[] secondaryState)
Set secondary part of the current state.
|
void |
setTime(double time)
Set current time.
|
public ExpandableStatefulODE(FirstOrderDifferentialEquations primary)
primary
- the primary set of differential equations to be integrated.public FirstOrderDifferentialEquations getPrimary()
public int getTotalDimension()
The complete set of equations correspond to the primary set plus all secondary sets.
public void computeDerivatives(double t, double[] y, double[] yDot) throws MaxCountExceededException, DimensionMismatchException
t
- current value of the independent time variabley
- array containing the current value of the complete state vectoryDot
- placeholder array where to put the time derivative of the complete state vectorMaxCountExceededException
- if the number of functions evaluations is exceededDimensionMismatchException
- if arrays dimensions do not match equations settingspublic int addSecondaryEquations(SecondaryEquations secondary)
secondary
- secondary equations setpublic EquationsMapper getPrimaryMapper()
getSecondaryMappers()
public EquationsMapper[] getSecondaryMappers()
getPrimaryMapper()
public void setTime(double time)
time
- current timepublic double getTime()
public void setPrimaryState(double[] primaryState) throws DimensionMismatchException
primaryState
- primary part of the current stateDimensionMismatchException
- if the dimension of the array does not
match the primary setpublic double[] getPrimaryState()
public double[] getPrimaryStateDot()
public void setSecondaryState(int index, double[] secondaryState) throws DimensionMismatchException
index
- index of the part to set as returned by addSecondaryEquations(SecondaryEquations)
secondaryState
- secondary part of the current stateDimensionMismatchException
- if the dimension of the partial state does not
match the selected equations set dimensionpublic double[] getSecondaryState(int index)
index
- index of the part to set as returned by addSecondaryEquations(SecondaryEquations)
public double[] getSecondaryStateDot(int index)
index
- index of the part to set as returned by addSecondaryEquations(SecondaryEquations)
public void setCompleteState(double[] completeState) throws DimensionMismatchException
completeState
- complete current state to copy data fromDimensionMismatchException
- if the dimension of the complete state does not
match the complete equations sets dimensionpublic double[] getCompleteState() throws DimensionMismatchException
DimensionMismatchException
- if the dimension of the complete state does not
match the complete equations sets dimensionCopyright © 2003–2016 The Apache Software Foundation. All rights reserved.