SparseGridsCalculateError.cpp File Reference

Calculate errors if an exact solution is provided. More...

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

Go to the source code of this file.

Functions
int	ExactSolution (void *, void *, double *, char *, int *)

Detailed Description

Calculate errors if an exact solution is provided.

Author

Debojyoti Ghosh, John Loffeld, Lee Ricketson

Definition in file SparseGridsCalculateError.cpp.

Function Documentation

ExactSolution()

int ExactSolution	(	void *	s,
		void *	m,
		double *	uex,
		char *	fname,
		int *	flag
	)

Read in the exact solution, if available.

Parameters

s	Solver object of type HyPar
m	MPI object of type MPIVariables
uex	Array to hold the exact solution, if available
fname	Filename root from which to read exact solution
flag	Flag to indicate if exact solution was available

Definition at line 16 of file ExactSolution.c.

{
  HyPar  *solver = (HyPar*) s;

  int same_size = 1;
  for (int i=0; i<solver->ndims; i++) {
    if (solver->dim_global[i] != solver->dim_global_ex[i]) same_size = 0;
  }

  if (same_size) {
    ReadArray( solver->ndims,
               solver->nvars,
               solver->dim_global,
               solver->dim_local,
               solver->ghosts,
               solver,
               m,
               NULL,
               uex,
               fname,
               flag);
  } else {
    ReadArraywInterp( solver->ndims,
                      solver->nvars,
                      solver->dim_global,
                      solver->dim_local,
                      solver->dim_global_ex,
                      solver->ghosts,
                      solver,
                      m,
                      NULL,
                      uex,
                      fname,
                      flag);
  }

  return(0);
}