TimeForwardEuler.c

Go to the documentation of this file.

#include <basic.h>
#include <arrayfunctions.h>
#include <simulation_object.h>
#include <timeintegration.h>

int TimeForwardEuler(
                      void *ts 
                    )
{
  TimeIntegration* TS = (TimeIntegration*) ts;
  SimulationObject* sim = (SimulationObject*) TS->simulation;
  int ns, nsims = TS->nsims;
  _DECLARE_IERR_;

  for (ns = 0; ns < nsims; ns++) {

    HyPar* solver = &(sim[ns].solver);
    MPIVariables* mpi = &(sim[ns].mpi);

    double* rhs = TS->rhs + TS->u_offsets[ns];
  
    /* Evaluate right-hand side */
    IERR TS->RHSFunction( rhs,
                          solver->u,
                          solver,
                          mpi,
                          TS->waqt  );   CHECKERR(ierr);
  
    /* update solution */
    _ArrayAXPY_(  rhs,
                  TS->dt,
                  solver->u,
                  TS->u_sizes[ns] );

    _ArrayScaleCopy1D_( solver->StageBoundaryIntegral,
                        TS->dt,
                        solver->StepBoundaryIntegral,
                        TS->bf_sizes[ns] );
  }

  return(0);
}