CalculateConservationError.c File Reference

Compute the conservation error. More...

#include <math.h>
#include <basic.h>
#include <mathfunctions.h>
#include <mpivars.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	CalculateConservationError (void *s, void *m)

Detailed Description

Compute the conservation error.

Author

Debojyoti Ghosh

Definition in file CalculateConservationError.c.

Function Documentation

int CalculateConservationError	(	void *	s,
		void *	m
	)

Calculates the error (L2) in conservation by computing the difference between the initial volume integral of the solution, and the sum of the current volume integral and the time integral of the boundary flux from the start of the simulation to the current simulation time.

Parameters

s	Solver object of type HyPar
m	MPI object of type MPIVariables

Definition at line 16 of file CalculateConservationError.c.

{
  HyPar         *solver = (HyPar*) s;
  int           v,nvars = solver->nvars;
  double        error;

  double base[nvars];
  for (v=0; v<nvars; v++) {
    if (absolute(solver->VolumeIntegralInitial[v]) > 1.0)
      base[v] = absolute(solver->VolumeIntegralInitial[v]);
    else base[v] = 1.0;
  }
  
  for (v=0; v<nvars; v++) {
    error =  (solver->VolumeIntegral[v]+solver->TotalBoundaryIntegral[v]-solver->VolumeIntegralInitial[v]) 
           * (solver->VolumeIntegral[v]+solver->TotalBoundaryIntegral[v]-solver->VolumeIntegralInitial[v]);
    solver->ConservationError[v] = sqrt(error)/base[v];
  }
  
  return(0);
}