dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #03108
New ODE solver interface
I fixed one final thing in the ODE solver interface. The initial
condition is now set by assigning to a given vector u:
void u0(uBlasVector& u)
{
u(0) = 1.0;
u(1) = 0.0;
}
instead of
real u0(uint i)
{
if ( i == 0 )
return 1.0;
else
return 0.0;
}
This is more convenient for both the ODE solver and the user...
This breaks the updated elasticity solver, the heat equation solver
and the TimeDependentPDE interface.
It also breaks most of the ODE demos and benchmarks but that will be
an easy fix.
I'm very happy with the current ODE solver interface, so this should
be the end of modifications to the ODE solver interface for some time.
/Anders