← Back to team overview

dolfin team mailing list archive

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