LinearADRJacobian.c

Go to the documentation of this file.

#include <mathfunctions.h>
#include <physicalmodels/linearadr.h>

int LinearADRAdvectionJacobian( double* Jac,  
                                double* u,    
                                void*   p,    
                                int     dir,  
                                int     nvars,
                                int     upw    )
{
  LinearADR *param = (LinearADR*) p;

  if (param->a) {
    *Jac =    (1-absolute(upw))*absolute(param->a[dir])
           +  absolute(upw) * (1+upw) * max(0,param->a[dir]) * 0.5
           -  absolute(upw) * (1-upw) * min(0,param->a[dir]) * 0.5 ;
  } else {
    /* no advection term */
    *Jac = 0.0;
  }

  return 0;
}

int LinearADRDiffusionJacobian( double* Jac,  
                                double* u,    
                                void*   p,    
                                int     dir,  
                                int     nvars  )
{
  LinearADR *param = (LinearADR*) p;

  int v;
  for (v = 0; v < nvars; v++) {
    Jac[nvars*v+v] = -param->d[nvars*dir+v];
  }

  return 0;
}