SparseGridsInitializationWrapup.cpp

Go to the documentation of this file.

#include <arrayfunctions.h>
#include <sparse_grids_simulation.h>

extern "C" int VolumeIntegral (double*,double*,void*,void*);

int SparseGridsSimulation::InitializationWrapup()
{
  int ierr;

  for (int n = 0; n < m_nsims_sg; n++) {

    if (!m_rank) {
      printf("Interpolating initial solution to sparse grids domain %d.\n", n);
    }
    interpolate( &m_sims_sg[n], m_sim_fg );

    HyPar* solver = &(m_sims_sg[n].solver);
    MPIVariables* mpi = &(m_sims_sg[n].mpi);

    int nvars = solver->nvars;
    int ghosts = solver->ghosts;
    int *dim_local = solver->dim_local;

    /* exchange MPI-boundary values of u between processors */
    MPIExchangeBoundariesnD(  m_ndims,
                              nvars,
                              dim_local,
                              ghosts,
                              mpi,
                              solver->u  );

    /* calculate volume integral of the initial solution */
    ierr = ::VolumeIntegral(solver->VolumeIntegralInitial,
                            solver->u,
                            solver,
                            mpi );
    if (ierr) return ierr;

    if (!mpi->rank) {
      printf("  Volume integral of the initial solution on sparse grids domain %d:\n", n);
      for (int d=0; d<nvars; d++) {
        printf("    %2d:  %1.16E\n",d,solver->VolumeIntegralInitial[d]);
      }
    }

    /* Set initial total boundary flux integral to zero */
    _ArraySetValue_( solver->TotalBoundaryIntegral, nvars, 0 );

  }

  return 0; 
}