HyPar  1.0
Finite-Difference Hyperbolic-Parabolic PDE Solver on Cartesian Grids
FPPowerSystemComputeDiffNumber.c
Go to the documentation of this file.
1 #include <stdlib.h>
2 #include <basic.h>
3 #include <mathfunctions.h>
4 #include <arrayfunctions.h>
6 #include <mpivars.h>
7 #include <hypar.h>
8 
9 double FPPowerSystemDissipationFunction(int,void*,double);
10 
11 double FPPowerSystemComputeDiffNumber(void *s,void *m,double dt,double t)
12 {
13  HyPar *solver = (HyPar*) s;
14  FPPowerSystem *params = (FPPowerSystem*)solver->physics;
15 
16  int ndims = solver->ndims;
17  int ghosts = solver->ghosts;
18  int *dim = solver->dim_local;
19 
20  double max_diff = 0;
21  int index[ndims];
22  int done = 0; _ArraySetValue_(index,ndims,0);
23  while (!done) {
24  double dxinv; _GetCoordinate_(0,index[0],dim,ghosts,solver->dxinv,dxinv);
25  double dyinv; _GetCoordinate_(1,index[1],dim,ghosts,solver->dxinv,dyinv);
26  double dissp_x= FPPowerSystemDissipationFunction(0,params,t);
27  double dissp_y= FPPowerSystemDissipationFunction(1,params,t);
28 
29  double local_diff_x = absolute(dissp_x) * dt * dxinv * dxinv;
30  double local_diff_y = absolute(dissp_y) * dt * dyinv * dyinv;
31 
32  if (local_diff_x > max_diff) max_diff = local_diff_x;
33  if (local_diff_y > max_diff) max_diff = local_diff_y;
34 
35  _ArrayIncrementIndex_(ndims,dim,index,done);
36  }
37 
38  return(max_diff);
39 }
#define absolute(a)
Definition: math_ops.h:32
MPI related function definitions.
Contains function definitions for common mathematical functions.
Some basic definitions and macros.
int ndims
Definition: hypar.h:26
double FPPowerSystemDissipationFunction(int, void *, double)
Structure containing all solver-specific variables and functions.
Definition: hypar.h:23
Contains structure definition for hypar.
#define _ArraySetValue_(x, size, value)
int * dim_local
Definition: hypar.h:37
#define _ArrayIncrementIndex_(N, imax, i, done)
#define _GetCoordinate_(dir, i, dim, ghosts, x, coord)
Definition: basic.h:31
void * physics
Definition: hypar.h:266
int ghosts
Definition: hypar.h:52
double FPPowerSystemComputeDiffNumber(void *s, void *m, double dt, double t)
Contains macros and function definitions for common array operations.
double * dxinv
Definition: hypar.h:110