ShallowWater2DJacobian.c

Go to the documentation of this file.

#include <mathfunctions.h>
#include <matmult_native.h>
#include <physicalmodels/shallowwater2d.h>

int ShallowWater2DJacobian(
                    double  *Jac, 
                    double  *u,   
                    void    *p,   
                    int     dir,  
                    int     nvars,
                    int     upw   
                   )
{
  ShallowWater2D  *param = (ShallowWater2D*) p;
  static double   R[_MODEL_NVARS_*_MODEL_NVARS_], D[_MODEL_NVARS_*_MODEL_NVARS_], 
                  L[_MODEL_NVARS_*_MODEL_NVARS_], DL[_MODEL_NVARS_*_MODEL_NVARS_];

  /* get the eigenvalues and left,right eigenvectors */
  _ShallowWater2DEigenvalues_      (u,D,param,dir);
  _ShallowWater2DLeftEigenvectors_ (u,L,param,dir);
  _ShallowWater2DRightEigenvectors_(u,R,param,dir);

  int aupw = absolute(upw), k;
  k = 0; D[k] = absolute( (1-aupw)*D[k] + 0.5*aupw*(1+upw)*max(0,D[k]) + 0.5*aupw*(1-upw)*min(0,D[k]) );
  k = 4; D[k] = absolute( (1-aupw)*D[k] + 0.5*aupw*(1+upw)*max(0,D[k]) + 0.5*aupw*(1-upw)*min(0,D[k]) );
  k = 8; D[k] = absolute( (1-aupw)*D[k] + 0.5*aupw*(1+upw)*max(0,D[k]) + 0.5*aupw*(1-upw)*min(0,D[k]) );

  MatMult3(_MODEL_NVARS_,DL,D,L);
  MatMult3(_MODEL_NVARS_,Jac,R,DL);

  return(0);
}