CalculateConservationError.c

Go to the documentation of this file.

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

int CalculateConservationError(
                                void *s, 
                                void *m  
                              )
{
  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);
}