HARM
harm and utilities
 All Data Structures Files Functions Variables Typedefs Macros Pages
Macros | Functions | Variables
docs/newcodediff.txt File Reference

Macros

#define FIELDTYPE   BLANDFORDQUAD | #define FIELDTYPE VERTFIELD
 
#define N1   64 |
 
#define N2   64 |
 
#define FIELDSTAGMEM   1
 
#define PRODUCTION   0 | #define PRODUCTION 1
 
#define LIMITDTWITHSOURCETERM   0
 
#define DODISS   1 | #define DODISS 0
 
#define DOLUMVSR   1 | #define DOLUMVSR 0
 
#define DODISSVSR   1 | #define DODISSVSR 0
 

Functions

interpline c interpline jmono
c interpline mono c interpline
para c interpline smono c
interpline smono h interpline
smono_static h interppoint c
interppoint para c jon_interp
c jon_interp_computepreprocess
static c void 		vec2vecortho (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE *gcon, FTYPE *gcov, FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, FTYPE *vec, FTYPE *vecortho)
 
static void vec2vecortho (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE(*gcon)[NDIM], FTYPE(*gcov)[NDIM], FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, FTYPE *vec, FTYPE *vecortho)
 
static void vB2poyntingdensity (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE *gcon, FTYPE *gcov, FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, int vectorcomponent, FTYPE *vecv, FTYPE *vecB, FTYPE *compout)
 
static void vB2poyntingdensity (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE(*gcon)[NDIM], FTYPE(*gcov)[NDIM], FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, int vectorcomponent, FTYPE *vecv, FTYPE *vecB, FTYPE *compout)
 
static void vecup2vecdowncomponent (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE *gcon, FTYPE *gcov, FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, int vectorcomponent, FTYPE *vec, FTYPE *compout)
 
static void vecup2vecdowncomponent (int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE(*gcon)[NDIM], FTYPE(*gcov)[NDIM], FTYPE gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM], int oldgridtype, int newgridtype, int vectorcomponent, FTYPE *vec, FTYPE *compout)
 
static void read_gdumpline (FILE *in, int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE *gcon, FTYPE *gcov, FTYPE *gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM])
 
static void read_gdumpline (FILE *in, int ti[], FTYPE X[], FTYPE V[], FTYPE(*conn)[NDIM][NDIM], FTYPE(*gcon)[NDIM], FTYPE(*gcov)[NDIM], FTYPE *gdet, FTYPE ck[], FTYPE(*dxdxp)[NDIM])
 
 DLOOP (jj, kk) ucov[jj]+ = ucon[kk]*gcov[jj][kk]
 
 DLOOPA (jj) DLOOPA(kk) fscanf(in
 
jon_interp_filter c
jon_interp_interpolationitself
c jon_interp_mnewt c
jon_interp_newt c kazfulleos c
kazfulleos decsglobalprivate h
kazfulleos defsglobalprivate h
kazfulleos defs h kazfulleos
eostablesdefs h kazfulleos
global h kazfulleos_set_arrays
c kazfulleos superdefs h
liaison c liaison_set_arrays c
lnsrch c lubksb c ludcmp c
main c math_tools c metric c
void 		gcov_func (struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE *gcovinfunc, FTYPE *gcovpertinfunc)|void gcov_func(struct of_geom *ptrgeom
 obtain gcov in primcoords of whichcoord type metric/coords here ptrgeom is only expected to contain i,j,k,p location More...
 
void set_gcov_ksmetric (FTYPE *V, FTYPE *gcovinfunc, FTYPE *gcovpertinfunc)
 (~t,r,,~) More...
 
void set_gcov_ksmetric (FTYPE *V, FTYPE(*gcovinfunc)[NDIM], FTYPE *gcovpertinfunc)
 
void gcov2gcovprim (struct of_geom *ptrgeom, FTYPE *X, FTYPE *V, FTYPE *gcovinfunc, FTYPE *gcovpertinfunc, FTYPE *gcovinfuncprim, FTYPE *gcovpertinfuncprim)
 
void gcov2gcovprim (struct of_geom *ptrgeom, FTYPE *X, FTYPE *V, FTYPE(*gcovinfunc)[NDIM], FTYPE *gcovpertinfunc, FTYPE(*gcovinfuncprim)[NDIM], FTYPE *gcovpertinfuncprim)
 

Variables

global openmploops h global
openmpsharedlist h global
other h global realdef h
global stepch h global storage
h global structs h FTYPE 		gcov [SYMMATRIXNDIM]
 
FTYPE gcon [SYMMATRIXNDIM]
 
FTYPE gengcov [SYMMATRIXNDIM]
 
FTYPE gengcon [SYMMATRIXNDIM]
 
idealgaseos c image c initbase
boundloop c initbase c 		DODIAGS =1
 
initbase defaultnprlists c nprstart =0
 
 nprend =8
 
 nprlist [8] =8
 
< nprfluxboundstart=0;nprfluxboundend=8;|nprfluxboundstart=0;nprfluxboundend=7;nprfluxboundlist[8]=8;
< initbase.enerregions.cinitbase.gridsectioning.cinitbase.tools.cinit.cdefcoord=9;|defcoord=JET3COORDS;a=0.9375;|a=0.92;R0=0.0;|Rout=40.0;|> 		R0 =0
 
 Rout =10
 
 hslope = 1.04*pow(h_over_r,2.0/3.0)
 
 FLUXB =FLUXCTTOTH
 
 tf = 2000.0
 
init ff h init fishmon c init
fishmon h init grb c 		pr [U3] = 4.0*dxdxp[3][3]*ptrgeom->gcon[GIND(3,3)]*sin(th)*sin(th)
 
init ns backup2 c ur = -geom.gcon[GIND(0,1)]/sqrt(-geom.gcon[GIND(0,0)])
 
 uh = -geom.gcon[GIND(0,2)]/sqrt(-geom.gcon[GIND(0,0)])
 
 up = -geom.gcon[GIND(0,3)]/sqrt(-geom.gcon[GIND(0,0)])
 
init ns h init readdata c przamobl [U1] = (geombl.gcon[GIND(0,1)])/(geombl.gcon[GIND(0,0)])
 
init rebecca c prlocal [U1] = -ptrgeom->gcon[GIND(0,1)]/sqrt(-ptrgeom->gcon[GIND(0,0)])
 
static void FTYPE X []
 
static void FTYPE FTYPE V []
 
static void FTYPE FTYPE FTYPE(* conn )[NDIM][NDIM]
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE gdet
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE 		ck []
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(* 		dxdxp )[NDIM]
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int 		oldgridtype
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int 		newgridtype
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int FTYPE * 		vec
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int FTYPE
FTYPE *static vecortho void
FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int int 		vectorcomponent
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int FTYPE
FTYPE *static vecortho void
FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int int FTYPE * 		vecv
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int FTYPE
FTYPE *static vecortho void
FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int int FTYPE
FTYPE * 		vecB
 
static void FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int FTYPE
FTYPE *static vecortho void
FTYPE FTYPE FTYPE(FTYPE(FTYPE(FTYPE
FTYPE FTYPE(int int int FTYPE
FTYPE FTYPE *compout FTYPE 		newgcov [SYMMATRIXNDIM] =0.0
 
static void int ti []
 
static void int FTYPE FTYPE
FTYPE(FTYPE(FTYPE(FTYPE FTYPE
FTYPE(*dxdxp) 		DLOOPA [NDIM](jj) DLOOPA(kk) fscanf(in
 
static void int FTYPE FTYPE
FTYPE(FTYPE(FTYPE(FTYPE FTYPE
FTYPE(*dxdxp) 		SCANARG [NDIM]
 
jon_interp_filter c
jon_interp_interpolationitself
c jon_interp_mnewt c
jon_interp_newt c kazfulleos c
kazfulleos decsglobalprivate h
kazfulleos defsglobalprivate h
kazfulleos defs h kazfulleos
eostablesdefs h kazfulleos
global h kazfulleos_set_arrays
c kazfulleos superdefs h
liaison c liaison_set_arrays c
lnsrch c lubksb c ludcmp c
main c math_tools c metric c
void int 		getprim
 
jon_interp_filter c
jon_interp_interpolationitself
c jon_interp_mnewt c
jon_interp_newt c kazfulleos c
kazfulleos decsglobalprivate h
kazfulleos defsglobalprivate h
kazfulleos defs h kazfulleos
eostablesdefs h kazfulleos
global h kazfulleos_set_arrays
c kazfulleos superdefs h
liaison c liaison_set_arrays c
lnsrch c lubksb c ludcmp c
main c math_tools c metric c
void int int 		whichcoord
 
jon_interp_filter c
jon_interp_interpolationitself
c jon_interp_mnewt c
jon_interp_newt c kazfulleos c
kazfulleos decsglobalprivate h
kazfulleos defsglobalprivate h
kazfulleos defs h kazfulleos
eostablesdefs h kazfulleos
global h kazfulleos_set_arrays
c kazfulleos superdefs h
liaison c liaison_set_arrays c
lnsrch c lubksb c ludcmp c
main c math_tools c metric c
void int int FTYPE FTYPE(* 		gcovinfunc )[NDIM] =localgcov[GIND(j,k)]
 

Detailed Description

~~~~ advance.c DLOOPA(jj) dphidt += GLOBALMETMACP0A3(conn,i,j,k,jj,TT,TT)*(ptrgeom->gcov[GIND(jj,TT)]); | DLOOPA(jj) dphidt += GLOBALMETMACP0A3(conn,i,j,k,jj,TT,TT)*(ptrgeom->gcov[jj][TT]); phi = -(1.0+ptrgeom->gcov[GIND(TT,TT)])*0.5; | phi = -(1.0+ptrgeom->gcov[TT][TT])*0.5; frac = fabs(ptrgeom->gcon[GIND(TT,TT)]*dphidt); | frac = fabs(ptrgeom->gcon[TT][TT]*dphidt); bcucof.c bcuint.c bin2txt.c boundmpi.c boundmpiint.c bounds.c bounds.ff.c pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[GIND(3,3)])/(geom.gcov[GIND(3,3)]); | pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[3][3])/(geom.gcov[3][3]); prface[B3] = rprim[B3]*(geomc.gcov[GIND(3,3)])/(geomf.gcov[GIND(3,3)]); | prface[B3] = rprim[B3]*(geomc.gcov[3][3])/(geomf.gcov[3][3]); pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[GIND(3,3)])/(geom.gcov[GIND(3,3)]); | pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[3][3])/(geom.gcov[3][3]); prface[B3]=MACP0A1(prim,ri,rj,rk,B3)*fabs(geomc.gcov[GIND(3,3)])/(geomf.gcov[GIND(3,3)]); | prface[B3]=MACP0A1(prim,ri,rj,rk,B3)*fabs(geomc.gcov[3][3])/(geomf.gcov[3][3]); pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[GIND(3,3)])/(geom.gcov[GIND(3,3)]); | pl = B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(rgeom.gcov[3][3])/(geom.gcov[3][3]); bounds.fishmon.c boundsflux.c bounds.grb.c boundsint.c bounds.ns.backup2.c bounds.ns.backup.c bounds.ns.c MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(sqrt(fabs(geom.gcov[GIND(2,2)]))*pow(V[1],3.0))/(sqrt(fabs(rgeom.gcov[GIND(2,2)]))*pow(rV[1],3.0)); | MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*(sqrt(fabs(geom.gcov[2][2]))*pow(V[1],3.0))/(sqrt(fabs(rgeom.gcov[2][2]))*pow(rV[1],3.0)); v1=prface[B1]*sqrt(fabs(fgeom.gcov[GIND(1,1)]))*pow(fV[1],3); // close to constant for dipole | v1=prface[B1]*sqrt(fabs(fgeom.gcov[1][1]))*pow(fV[1],3); // close to constant for dipole v2=MACP0A1(prim,ri,rj,rk,B1)*sqrt(fabs(rgeom.gcov[GIND(1,1)]))*pow(rV[1],3); | v2=MACP0A1(prim,ri,rj,rk,B1)*sqrt(fabs(rgeom.gcov[1][1]))*pow(rV[1],3); MACP0A1(prim,i,j,k,B1) = newv/(sqrt(fabs(geom.gcov[GIND(1,1)]))*pow(V[1],3)); | MACP0A1(prim,i,j,k,B1) = newv/(sqrt(fabs(geom.gcov[1][1]))*pow(V[1],3)); v1=prface[B1]*sqrt(fabs(fgeom.gcov[GIND(1,1)]))*pow(fV[1],3); // close to constant for dipole | v1=prface[B1]*sqrt(fabs(fgeom.gcov[1][1]))*pow(fV[1],3); // close to constant for dipole v2=MACP0A1(prim,ri,rj,rk,B1)*sqrt(fabs(rgeom.gcov[GIND(1,1)]))*pow(rV[1],3); | v2=MACP0A1(prim,ri,rj,rk,B1)*sqrt(fabs(rgeom.gcov[1][1]))*pow(rV[1],3); v3=MACP0A1(prim,ri+1,rj,rk,B1)*sqrt(fabs(rrgeom.gcov[GIND(1,1)]))*pow(rrV[1],3); | v3=MACP0A1(prim,ri+1,rj,rk,B1)*sqrt(fabs(rrgeom.gcov[1][1]))*pow(rrV[1],3); MACP0A1(prim,i,j,k,B1) = newv/(sqrt(fabs(geom.gcov[GIND(1,1)]))*pow(V[1],3)); | MACP0A1(prim,i,j,k,B1) = newv/(sqrt(fabs(geom.gcov[1][1]))*pow(V[1],3)); up2face+=prface[U1+jj-1]*prface[U1+kk-1]*fgeom.gcov[GIND(jj,kk)]; | up2face+=prface[U1+jj-1]*prface[U1+kk-1]*fgeom.gcov[jj][kk]; Bp2face+=prface[B1+jj-1]*prface[B1+kk-1]*fgeom.gcov[GIND(jj,kk)]; | Bp2face+=prface[B1+jj-1]*prface[B1+kk-1]*fgeom.gcov[jj][kk]; up2cent+=MACP0A1(prim,i,j,k,U1+jj-1)*MACP0A1(prim,i,j,k,U1+kk-1)*geom.gcov[GIND(jj,kk)]; | up2cent+=MACP0A1(prim,i,j,k,U1+jj-1)*MACP0A1(prim,i,j,k,U1+kk-1)*geom.gcov[jj][kk]; Bp2cent+=MACP0A1(prim,i,j,k,B1+jj-1)*MACP0A1(prim,i,j,k,B1+kk-1)*geom.gcov[GIND(jj,kk)]; | Bp2cent+=MACP0A1(prim,i,j,k,B1+jj-1)*MACP0A1(prim,i,j,k,B1+kk-1)*geom.gcov[jj][kk]; prface[B1] = GLOBALMACP0A1(panalytic,i,j,k,B1)*sqrt(fabs(geomc.gcov[GIND(1,1)]))*pow(Vc[1],3)/(sqrt(fabs(geomf.gcov[GIND(1,1)]))*pow(Vf[1],3)); // set from offset of analytical solution | prface[B1] = GLOBALMACP0A1(panalytic,i,j,k,B1)*sqrt(fabs(geomc.gcov[1][1]))*pow(Vc[1],3)/(sqrt(fabs(geomf.gcov[1][1]))*pow(Vf[1],3)); // set from offset of analytical solution prface[B1] = MACP0A1(prim,i,j,k,B1)*sqrt(fabs(geomc.gcov[GIND(1,1)]))*pow(Vc[1],3)/(sqrt(fabs(geomf.gcov[GIND(1,1)]))*pow(Vf[1],3)); | // prface[B1] = MACP0A1(prim,i,j,k,B1)*sqrt(fabs(geomc.gcov[1][1]))*pow(Vc[1],3)/(sqrt(fabs(geomf.gcov[1][1]))*pow(Vf[1],3)); prface[B2] = MACP0A1(prim,i,j,k,B2)*(sqrt(fabs(geomc.gcov[GIND(2,2)]))*pow(Vc[1],3.0))/(sqrt(fabs(geomf.gcov[GIND(2,2)]))*pow(Vf[1],3.0)); | prface[B2] = MACP0A1(prim,i,j,k,B2)*(sqrt(fabs(geomc.gcov[2][2]))*pow(Vc[1],3.0))/(sqrt(fabs(geomf.gcov[2][2]))*pow(Vf[1],3.0)); bounds.nsold.c bounds.rebecca.c bounds.sasha.c bounds.tools.c dualfprintf(fail_file,"JUST BEFORE INFLOWCHECK: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) *sqrt(geom[U1].gcov[GIND(1,1)]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[GIND(2,2)]),MACP0A1(prim,i,j,k,U3) *sqrt(geom | // dualfprintf(fail_file,"JUST BEFORE INFLOWCHECK: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) *sqrt(geom[U1].gcov[1][1]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[2][2]),MACP0A1(prim,i,j,k,U3) *sqrt(geom[U1].gcov[ DLOOP(jj,kk) dualfprintf(fail_file,"gcov[%d][%d]=%21.15g\n",jj,kk,geom[U1].gcov[GIND(jj,kk)]); | // DLOOP(jj,kk) dualfprintf(fail_file,"gcov[%d][%d]=%21.15g\n",jj,kk,geom[U1].gcov[jj][kk]); dualfprintf(fail_file,"JUST BEFORE LIMIT: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) *sqrt(geom[U1].gcov[GIND(1,1)]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[GIND(2,2)]),MACP0A1(prim,i,j,k,U3) *sqrt(geom[U1].g | // dualfprintf(fail_file,"JUST BEFORE LIMIT: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) *sqrt(geom[U1].gcov[1][1]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[2][2]),MACP0A1(prim,i,j,k,U3) *sqrt(geom[U1].gcov[3][3]) dualfprintf(fail_file,"JUST AFTER LIMIT: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) *sqrt(geom[U1].gcov[GIND(1,1)]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[GIND(2,2)]),MACP0A1(prim,i,j,k,U3) *sqrt(geom[U1].gc | // dualfprintf(fail_file,"JUST AFTER LIMIT: i=%d j=%d k=%d prim[U1]=%21.15g prim[U2]=%21.15g prim[U3]=%21.15g\n",i,j,k,MACP0A1(prim,i,j,k,U1) sqrt(geom[U1].gcov[1][1]),MACP0A1(prim,i,j,k,U2) *sqrt(geom[U1].gcov[2][2]),MACP0A1(prim,i,j,k,U3) *sqrt(geom[U1].gcov[3][3])) Bd3=0.0; SLOOPA(jj) Bd3 += MACP0A1(prim,ri,rj,rk,B1+jj-1)(ptrrgeom[B1+jj-1]->gcov[GIND(3,jj)]); | Bd3=0.0; SLOOPA(jj) Bd3 += MACP0A1(prim,ri,rj,rk,B1+jj-1)*(ptrrgeom[B1+jj-1]->gcov[3][jj]); Bd3ri2=0.0; SLOOPA(jj) Bd3ri2 += MACP0A1(prim,ri2,rj,rk,B1+jj-1)*(ptrri2geom[B1+jj-1]->gcov[GIND(3,jj)]); | Bd3ri2=0.0; SLOOPA(jj) Bd3ri2 += MACP0A1(prim,ri2,rj,rk,B1+jj-1)*(ptrri2geom[B1+jj-1]->gcov[3][jj]); Bd3ri3=0.0; SLOOPA(jj) Bd3ri3 += MACP0A1(prim,ri3,rj,rk,B1+jj-1)*(ptrri3geom[B1+jj-1]->gcov[GIND(3,jj)]); | Bd3ri3=0.0; SLOOPA(jj) Bd3ri3 += MACP0A1(prim,ri3,rj,rk,B1+jj-1)*(ptrri3geom[B1+jj-1]->gcov[3][jj]); gcon03=ptrgeom[B3]->gcon[GIND(0,3)]; | gcon03=ptrgeom[B3]->gcon[0][3]; gcon13=ptrgeom[B3]->gcon[GIND(1,3)]; | gcon13=ptrgeom[B3]->gcon[1][3]; gcon23=ptrgeom[B3]->gcon[GIND(2,3)]; | gcon23=ptrgeom[B3]->gcon[2][3]; gcon33=ptrgeom[B3]->gcon[GIND(3,3)]; | gcon33=ptrgeom[B3]->gcon[3][3]; gcov01=ptrgeom[B3]->gcov[GIND(0,1)]; | gcov01=ptrgeom[B3]->gcov[0][1]; gcov02=ptrgeom[B3]->gcov[GIND(0,2)]; | gcov02=ptrgeom[B3]->gcov[0][2]; gcov11=ptrgeom[B3]->gcov[GIND(1,1)]; | gcov11=ptrgeom[B3]->gcov[1][1]; gcov12=gcov21=ptrgeom[B3]->gcov[GIND(1,2)]; | gcov12=gcov21=ptrgeom[B3]->gcov[1][2]; gcov22=ptrgeom[B3]->gcov[GIND(2,2)]; | gcov22=ptrgeom[B3]->gcov[2][2]; gcov03=ptrgeom[B3]->gcov[GIND(0,3)]; | gcov03=ptrgeom[B3]->gcov[0][3]; gcov13=ptrgeom[B3]->gcov[GIND(1,3)]; | gcov13=ptrgeom[B3]->gcov[1][3]; gcov23=ptrgeom[B3]->gcov[GIND(2,3)]; | gcov23=ptrgeom[B3]->gcov[2][3]; pl=B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*fabs((ptrrgeom[pl]->gcov[GIND(3,3)])/(ptrgeom[pl]->gcov[GIND(3,3)])); | // pl=B3; MACP0A1(prim,i,j,k,pl) = MACP0A1(prim,ri,rj,rk,pl)*fabs((ptrrgeom[pl]->gcov[3][3])/(ptrgeom[pl]->gcov[3][3])); MACP0A1(prim,i,j,k,pl) = ftemp*fabs((ptrrgeom[pl]->gcov[GIND(3,3)])/(ptrgeom[pl]->gcov[GIND(3,3)])); | MACP0A1(prim,i,j,k,pl) = ftemp*fabs((ptrrgeom[pl]->gcov[3][3])/(ptrgeom[pl]->gcov[3][3])); broydn.c coord.c coord.h copyandinit_functions.c definit.h defs.general.h defs.grmhd.h defs.h@ defs.jon_interp.h defs.liaison.h defs.user.h defs.user.sasha.h diag.c dudp_calc_3vel.c ducov_dv[i][j] += geom->gcov[GIND(i,k)] * ducon_dv[k][j]; | ducov_dv[i][j] += geom->gcov[i][k] * ducon_dv[k][j]; dbcov_dv[i][j] += geom->gcov[GIND(i,k)] * dbcon_dv[k][j]; | dbcov_dv[i][j] += geom->gcov[i][k] * dbcon_dv[k][j]; dudp_calc.c SLOOP(j,k) dgdvi[j]+=1.0/gamma*ptrgeom->gcov[GIND(j,k)]*pr[U1+k-1]; | SLOOP(j,k) dgdvi[j]+=1.0/gamma*ptrgeom->gcov[j][k]*pr[U1+k-1]; alpha=1.0/sqrt(-ptrgeom->gcon[GIND(TT,TT)]); | //alpha=1.0/sqrt(-ptrgeom->gcon[TT][TT]); SLOOPA(j) betacon[j]=ptrgeom->gcon[GIND(TT,j)]*alpha*alpha; | SLOOPA(j) betacon[j]=ptrgeom->gcon[TT][j]*alpha*alpha; DLOOP(j,k) n += ptrgeom->gcov[GIND(j,k)]*vcon1[j]*vcon2[k] ; | DLOOP(j,k) n += ptrgeom->gcov[j][k]*vcon1[j]*vcon2[k] ; DLOOP(j,k) n += ptrgeom->gcon[GIND(j,k)]*vcov1[j]*vcov2[k] ; | DLOOP(j,k) n += ptrgeom->gcon[j][k]*vcov1[j]*vcov2[k] ; DSLOOP(j,k) duudud[j][k] = ptrgeom->gcon[GIND(j,k)] + ptrgeom->gcon[GIND(j,TT)]*(-ucon[k]/ucon[TT]) ; | DSLOOP(j,k) duudud[j][k] = ptrgeom->gcon[j][k] + ptrgeom->gcon[j][TT]*(-ucon[k]/ucon[TT]) ; DSLOOP(j,k) dudduu[j][k] = ptrgeom->gcov[GIND(j,k)] + ptrgeom->gcov[GIND(j,TT)]*dutdui[k]; | DSLOOP(j,k) dudduu[j][k] = ptrgeom->gcov[j][k] + ptrgeom->gcov[j][TT]*dutdui[k]; DSLOOP(j,k) for(l=0;l<NDIM;l++) dbdiduj[j][k] += ptrgeom->gcov[GIND(j,l)]*dbiduj[l][k]; | DSLOOP(j,k) for(l=0;l<NDIM;l++) dbdiduj[j][k] += ptrgeom->gcov[j][l]*dbiduj[l][k]; dump.c ptrftemp=(FTYPE*)(&localgcon[GIND(0,0)]); | ptrftemp=(FTYPE*)(&localgcon[0][0]); ptrftemp=(FTYPE*)(&localgcov[GIND(0,0)]); | ptrftemp=(FTYPE*)(&localgcov[0][0]); dump_ener.c dumpgen.c eos.c extractnonfunc.c f2c.h fail.c fdjac.c fixup.c alpha = 1. / sqrt(-ptrgeom->gcon[GIND(0,0)]); | // alpha = 1. / sqrt(-ptrgeom->gcon[0][0]); alpha = alpha = 1./sqrt(-ptrgeom->gcon[GIND(TT,TT)]) ; | // alpha = alpha = 1./sqrt(-ptrgeom->gcon[TT][TT]) ; alpha = 1./sqrt(-ptrgeom->gcon[GIND(TT,TT)]) ; | //alpha = 1./sqrt(-ptrgeom->gcon[TT][TT]) ; alpha = alpha = 1./sqrt(-ptrgeom->gcon[GIND(TT,TT)]) ; | // alpha = alpha = 1./sqrt(-ptrgeom->gcon[TT][TT]) ; probs[U1+i-1] = ptrgeom->gcon[GIND(TT,i)]/ptrgeom->gcon[GIND(TT,TT)] ; | probs[U1+i-1] = ptrgeom->gcon[TT][i]/ptrgeom->gcon[TT][TT] ; gradient[i]=2.0*(ptrgeom->gcov[GIND(0,i)]); | gradient[i]=2.0*(ptrgeom->gcov[0][i]); gradient[i]+=2.0*ucon[j]*ptrgeom->gcov[GIND(i,j)]; | gradient[i]+=2.0*ucon[j]*ptrgeom->gcov[i][j]; alpha = 1./sqrt(-ptrgeom->gcon[GIND(0,0)]) ; | alpha = 1./sqrt(-ptrgeom->gcon[0][0]) ; betacon = ptrgeom->gcon[GIND(0,1)]*alpha*alpha ; | betacon = ptrgeom->gcon[0][1]*alpha*alpha ; betacon = ptrgeom->gcon[GIND(0,2)]*alpha*alpha ; | betacon = ptrgeom->gcon[0][2]*alpha*alpha ; betacon = ptrgeom->gcon[GIND(0,3)]*alpha*alpha ; | betacon = ptrgeom->gcon[0][3]*alpha*alpha ; SLOOP(j,k) vsq += ptrgeom->gcov[GIND(j,k)]*pr[U1+j-1]*pr[U1+k-1] ; | SLOOP(j,k) vsq += ptrgeom->gcov[j][k]*pr[U1+j-1]*pr[U1+k-1] ; flux.c fluxcompute.c fluxct.c fluxctstag.c flux.mergedc2ea2cmethod.c fluxvpot.c fmin.c freespace.c gaussj.c generatenprs.c global.bounds.h global.comploops.h global.depnmemonics.h global.dump.h global.fieldmacros.h global.funcdeclare.h int interpX_gcov(FTYPE *X, struct of_compgeom (*compgeom)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE (*gcovgrid)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3][SYMMATRIXNDIM], FTYPE (*gcovpertgrid)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3 | // int interpX_gcov(FTYPE *X, struct of_compgeom (*compgeom)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3], FTYPE (*gcovgrid)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SHIFTSTORE3][NDIM][NDIM], FTYPE (*gcovpertgrid)[NSTORE1+SHIFTSTORE1][NSTORE2+SHIFTSTORE2][NSTORE3+SH int interpX_gcov(FTYPE *X, struct of_compgeom PTRDEFMETMACP1A0(compgeom,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3), FTYPE PTRDEFMETMACP1A2(gcovgrid,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3,NDIM,NDIM), FTYPE PTRDEFMETMACP1A1(gcovpertgrid,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3,NDIM), FTYPE *gcov, FTYPE *g | int interpX_gcov(FTYPE *X, struct of_compgeom PTRDEFMETMACP1A0(compgeom,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3), FTYPE PTRDEFMETMACP1A2(gcovgrid,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3,NDIM,NDIM), FTYPE PTRDEFMETMACP1A1(gcovpertgrid,FILL,N1M+SHIFT1,N2M+SHIFT2,N3M+SHIFT3,NDIM), FTYPE (*gcov)[NDIM], extern int assignmetricstorage_new(struct of_compgeom *mygeom, FTYPE **localgcov, FTYPE **localgcon, FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); | extern int assignmetricstorage_new(struct of_compgeom *mygeom, FTYPE (**localgcov)[NDIM], FTYPE (**localgcon)[NDIM], FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); extern int assignmetricstorage_old(int loc, int i, int j, int k, FTYPE **localgcov, FTYPE **localgcon, FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); | extern int assignmetricstorage_old(int loc, int i, int j, int k, FTYPE (**localgcov)[NDIM], FTYPE (**localgcon)[NDIM], FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); extern int assignmetricstorage_oldlast(int loc, int i, int j, int k, FTYPE **localgcov, FTYPE **localgcon, FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); | extern int assignmetricstorage_oldlast(int loc, int i, int j, int k, FTYPE (**localgcov)[NDIM], FTYPE (**localgcon)[NDIM], FTYPE **localgcovpert, FTYPE **localgdet, FTYPE **localgdetvol, FTYPE **localalphalapse, FTYPE **localbetasqoalphasq, FTYPE **beta, FTYPE **localeomfunc); FTYPE (*localgcov)[NDIM]; \ < FTYPE (*localgcon)[NDIM]; \ < FTYPE *localgcov; \ | FTYPE (*localgcov)[NDIM];\ FTYPE *localgcon;\ | FTYPE (*localgcon)[NDIM];\ extern FTYPE gdet_func(int whichcoord, FTYPE *gcov); | extern FTYPE gdet_func(int whichcoord, FTYPE (*gcov)[NDIM]); extern FTYPE gdet_func_singcheck(int whichcoord, FTYPE *V, FTYPE *gcov); | extern FTYPE gdet_func_singcheck(int whichcoord, FTYPE *V, FTYPE (*gcov)[NDIM]); extern void bl_gcov_func(FTYPE r, FTYPE th, FTYPE *gcov); | //extern void bl_gcov_func(FTYPE r, FTYPE th, FTYPE (*gcov)[NDIM]); extern void bl_gcon_func(FTYPE r, FTYPE th, FTYPE *gcon); | //extern void bl_gcon_func(FTYPE r, FTYPE th, FTYPE (*gcon)[NDIM]); extern void gcov_func(int getprim, int whichcoord, FTYPE *X, FTYPE *gcov); | //extern void gcov_func(int getprim, int whichcoord, FTYPE *X, FTYPE (*gcov)[NDIM]); extern void gcon_func(int getprim, int whichcoord, FTYPE *X, FTYPE *gcov, FTYPE *gcon); | //extern void gcon_func(int getprim, int whichcoord, FTYPE *X, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM]); extern void matrix_inverse_metric(int whichcoord, FTYPE *gcov, FTYPE *gcon); | extern void matrix_inverse(int whichcoord, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM]); extern void matrix_inverse(int whichcoord, FTYPE (*genmatrixlower)[NDIM], FTYPE (*genmatrixupper)[NDIM]); | extern void alphalapse_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM], FTYPE *alphalapse); extern void alphalapse_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE *gcov, FTYPE *gcon, FTYPE *alphalapse); | extern void betasqoalphasq_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM], FTYPE *betasqoalphasq); extern void betasqoalphasq_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE *gcov, FTYPE *gcon, FTYPE *betasqoalphasq); | extern void beta_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM], FTYPE alphalapse, FTYPE *beta); extern void beta_func(struct of_geom *ptrgeom, int getprim, int whichcoord, FTYPE *X, FTYPE *gcov, FTYPE *gcon, FTYPE alphalapse, FTYPE *beta); < extern int tetr_func(int inputtype, FTYPE *gcov, FTYPE (*tetr_cov)[NDIM],FTYPE (*tetr_con)[NDIM], FTYPE eigenvalues[]); | extern int tetr_func(int inputtype, FTYPE (*gcov)[NDIM], FTYPE (*tetr_cov)[NDIM],FTYPE (*tetr_con)[NDIM], FTYPE eigenvalues[]); extern int tetr_func_frommetric(FTYPE (*dxdxp)[NDIM], FTYPE *gcov, FTYPE (*tetrcov)[NDIM],FTYPE (*tetrcon)[NDIM], FTYPE eigenvalues[]); | extern int tetr_func_frommetric(FTYPE (*dxdxp)[NDIM], FTYPE (*gcov)[NDIM], FTYPE (*tetrcov)[NDIM],FTYPE (*tetrcon)[NDIM], FTYPE eigenvalues[]); extern void SHOULDNOTREACHHEREEVERBUGYOUHAVE(void); < global.funcdeclare.user.h global.funcdeclare.user.sasha.h global.general.h global.gridsectioning.h global.grmhd.h global.h@ global.inits.h global.jon_interp.h extern void matrix_inverse(int whichcoord, FTYPE *gcov, FTYPE *gcon); | extern void matrix_inverse(int whichcoord, FTYPE (*gcov)[NDIM], FTYPE (*gcon)[NDIM]); extern int tetr_func(int inputtype, FTYPE *gcov, FTYPE (*tetr_cov)[NDIM],FTYPE (*tetr_con)[NDIM], FTYPE eigenvalues[]); | extern int tetr_func(int inputtype, FTYPE (*gcov)[NDIM], FTYPE (*tetr_cov)[NDIM],FTYPE (*tetr_con)[NDIM], FTYPE eigenvalues[]); extern int tetr_func_frommetric(FTYPE (*dxdxp)[NDIM], FTYPE *gcov, FTYPE (*tetrcov)[NDIM],FTYPE (*tetrcon)[NDIM], FTYPE eigenvalues[]); | extern int tetr_func_frommetric(FTYPE (*dxdxp)[NDIM], FTYPE (*gcov)[NDIM], FTYPE (*tetrcov)[NDIM],FTYPE (*tetrcon)[NDIM], FTYPE eigenvalues[]); global.liaison.h global.loops.boundaries.h global.loops.diagnostics.h global.loops.h global.loops.manypoints1d.h global.loops.manypoints.h global.loops.perpoint.h global.mpi_grmhd_grray_liaison.h global.nondepnmemonics.h < define how to access symmetric matrices with size 4x4 without redundant elements < < 1) regexp: gcov[([_>a-zA-Z0-9+-\ ()]+)][([_>a-zA-Z0-9+-\ ()]+)] -> gcov[GIND(\1,\2)] < Then revert: gcov[GIND( <em>NDIM *, *NDIM *)] -> gcov[SYMMATRIXNDIM] since otherwise will be 1 larger than required < also need to catch: localgcon gcon gcovinfunc gcovtovks gcovbhks gcovmcoord gcovmid tmpgcov gcovprim tmpgcon glgen ghgen < 2) Then need to replace any multi-D pointer arg type with simple arg type: < FTYPE (**localgcov)[NDIM] -> FTYPE **localgcov < 3) Then need to get all global variables with [NDIM][NDIM] < 4) Also get *gcov type things: (*gcov)[NDIM] -> *gcov < also need to catch: gcon gcovinfunc gcovinfuncprim gcovprim gconprim gcovmcoord gcovptr gconptr < 5) If making assignment TO something using GIND(), then must control loops. < E.g. tetrad.c:tetr_func_frommetric(): newgcov[GIND(jj,kk)] += < Find some maybe by doing: grep -e "+ \{0,\}=" *.c *.h | grep GIND < < below are non-conditional ways of getting same result as: < #define GCOVI(i,j) (i>=j) ? i : j < #define GCOVJ(i,j) (i>=j) ? j : i < #define GCOVI(i,j) ((i>=j)(i-j) + j) < #define GCOVJ(i,j) ((i>=j)*(j-i) + i) < #define GIND(i,j) GCOVJ(i,j)*4 + GCOVI(i,j) - MAX(GCOVJ(i,j),0) - MAX(GCOVJ(i,j)-1,0) - MAX(GCOVJ(i,j)-2,0) < must multiply assignments by the below so don't duplicate sums < #define GINDASSIGNFACTOR(i,j) (1.0*(i>=j)) < #if(PRODUCTION==0) < #define GINDASSIGNMAC(name,i,j) (i>=j ? name[GIND(i,j)] : SHOULDNOTREACHHEREEVERBUGYOUHAVE()) < #else < #define GINDASSIGNMAC(name,i,j) name[GIND(i,j)] < #endif < < < < #define NDIM 4 /* number of total dimensions. Never changes

Definition in file newcodediff.txt.

Macro Definition Documentation

#define DODISS   1 | #define DODISS 0

Definition at line 291 of file newcodediff.txt.

#define DODISSVSR   1 | #define DODISSVSR 0

Definition at line 293 of file newcodediff.txt.

#define DOLUMVSR   1 | #define DOLUMVSR 0

Definition at line 292 of file newcodediff.txt.

#define FIELDSTAGMEM   1

Definition at line 288 of file newcodediff.txt.

#define FIELDTYPE   BLANDFORDQUAD | #define FIELDTYPE VERTFIELD

Definition at line 264 of file newcodediff.txt.

#define LIMITDTWITHSOURCETERM   0

Definition at line 290 of file newcodediff.txt.

#define N1   64 |

Definition at line 285 of file newcodediff.txt.

#define N2   64 |

Definition at line 286 of file newcodediff.txt.

#define PRODUCTION   0 | #define PRODUCTION 1

Definition at line 289 of file newcodediff.txt.

Function Documentation

DLOOP ( jj  ,
kk   
) = ucon[kk]*gcov[jj][kk]
DLOOPA ( jj  )
void gcov2gcovprim ( struct of_geom *  ptrgeom,
FTYPE *  X,
FTYPE *  V,
FTYPE *  gcovinfunc,
FTYPE *  gcovpertinfunc,
FTYPE *  gcovinfuncprim,
FTYPE *  gcovpertinfuncprim 
)

Definition at line 1580 of file metric.tools.c.

Here is the call graph for this function:

Here is the caller graph for this function:

void gcov2gcovprim ( struct of_geom *  ptrgeom,
FTYPE *  X,
FTYPE *  V,
FTYPE(*)  gcovinfunc[NDIM],
FTYPE *  gcovpertinfunc,
FTYPE(*)  gcovinfuncprim[NDIM],
FTYPE *  gcovpertinfuncprim 
)
jon_interp_filter c jon_interp_interpolationitself c jon_interp_mnewt c jon_interp_newt c kazfulleos c kazfulleos decsglobalprivate h kazfulleos defsglobalprivate h kazfulleos defs h kazfulleos eostablesdefs h kazfulleos global h kazfulleos_set_arrays c kazfulleos superdefs h liaison c liaison_set_arrays c lnsrch c lubksb c ludcmp c main c math_tools c metric c void gcov_func ( struct of_geom *  ptrgeom,
int  getprim,
int  whichcoord,
FTYPE *  X,
FTYPE *  gcovinfunc,
FTYPE *  gcovpertinfunc 
)

obtain gcov in primcoords of whichcoord type metric/coords here ptrgeom is only expected to contain i,j,k,p location

Definition at line 24 of file metric.c.

Here is the call graph for this function:

static void read_gdumpline ( FILE *  in,
int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE(*)  gcon[NDIM],
FTYPE(*)  gcov[NDIM],
FTYPE *  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM] 
)
static
static void read_gdumpline ( FILE *  in,
int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE *  gcon,
FTYPE *  gcov,
FTYPE *  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM] 
)
static
void set_gcov_ksmetric ( FTYPE *  V,
FTYPE *  gcovinfunc,
FTYPE *  gcovpertinfunc 
)

(~t,r,,~)

Definition at line 1492 of file metric.c.

Here is the caller graph for this function:

void set_gcov_ksmetric ( FTYPE *  V,
FTYPE(*)  gcovinfunc[NDIM],
FTYPE *  gcovpertinfunc 
)
static void vB2poyntingdensity ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE(*)  gcon[NDIM],
FTYPE(*)  gcov[NDIM],
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
int  vectorcomponent,
FTYPE *  vecv,
FTYPE *  vecB,
FTYPE *  compout 
)
static
static void vB2poyntingdensity ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE *  gcon,
FTYPE *  gcov,
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
int  vectorcomponent,
FTYPE *  vecv,
FTYPE *  vecB,
FTYPE *  compout 
)
static
static void vec2vecortho ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE *  gcon,
FTYPE *  gcov,
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
FTYPE *  vec,
FTYPE *  vecortho 
)
static
static void vec2vecortho ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE(*)  gcon[NDIM],
FTYPE(*)  gcov[NDIM],
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
FTYPE *  vec,
FTYPE *  vecortho 
)
static
static void vecup2vecdowncomponent ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE(*)  gcon[NDIM],
FTYPE(*)  gcov[NDIM],
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
int  vectorcomponent,
FTYPE *  vec,
FTYPE *  compout 
)
static
static void vecup2vecdowncomponent ( int  ti[],
FTYPE  X[],
FTYPE  V[],
FTYPE(*)  conn[NDIM][NDIM],
FTYPE *  gcon,
FTYPE *  gcov,
FTYPE  gdet,
FTYPE  ck[],
FTYPE(*)  dxdxp[NDIM],
int  oldgridtype,
int  newgridtype,
int  vectorcomponent,
FTYPE *  vec,
FTYPE *  compout 
)
static

Variable Documentation

void int FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE ck[]

Definition at line 367 of file newcodediff.txt.

static void int FTYPE FTYPE FTYPE(* conn

Definition at line 367 of file newcodediff.txt.

DLOOPA

Definition at line 379 of file newcodediff.txt.

idealgaseos c image c initbase boundloop c initbase c DODIAGS =1

Definition at line 236 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE(* dxdxp)[NDIM]

Definition at line 367 of file newcodediff.txt.

FLUXB =FLUXCTTOTH

Definition at line 262 of file newcodediff.txt.

& gcon

Definition at line 216 of file newcodediff.txt.

& gcov

Definition at line 214 of file newcodediff.txt.

gcovinfunc[j][k] =localgcov[GIND(j,k)]

Definition at line 402 of file newcodediff.txt.

void int FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE* gdet

Definition at line 367 of file newcodediff.txt.

FTYPE gengcon

Definition at line 218 of file newcodediff.txt.

FTYPE gengcov

Definition at line 217 of file newcodediff.txt.

jon_interp_filter c jon_interp_interpolationitself c jon_interp_mnewt c jon_interp_newt c kazfulleos c kazfulleos decsglobalprivate h kazfulleos defsglobalprivate h kazfulleos defs h kazfulleos eostablesdefs h kazfulleos global h kazfulleos_set_arrays c kazfulleos superdefs h liaison c liaison_set_arrays c lnsrch c lubksb c ludcmp c main c math_tools c metric c void int getprim

Definition at line 402 of file newcodediff.txt.

hslope = 1.04*pow(h_over_r,2.0/3.0)

Definition at line 261 of file newcodediff.txt.

FTYPE newgcov =0.0

Definition at line 368 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int newgridtype

Definition at line 367 of file newcodediff.txt.

nprend =8

Definition at line 238 of file newcodediff.txt.

nprlist[8] =8

Definition at line 239 of file newcodediff.txt.

nprstart =0

Definition at line 238 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int oldgridtype

Definition at line 367 of file newcodediff.txt.

pr[U3] = 4.0*dxdxp[3][3]*ptrgeom->gcon[GIND(3,3)]*sin(th)*sin(th)

Definition at line 276 of file newcodediff.txt.

init rebecca h init sasha c init sasha h init tools c prlocal = -ptrgeom->gcon[GIND(0,1)]/sqrt(-ptrgeom->gcon[GIND(0,0)])

Definition at line 327 of file newcodediff.txt.

przamobl[U3] = (geombl.gcon[GIND(0,1)])/(geombl.gcon[GIND(0,0)])

Definition at line 323 of file newcodediff.txt.

<nprfluxboundstart=0;nprfluxboundend=8;|nprfluxboundstart=0;nprfluxboundend=7;nprfluxboundlist[8]=8;<initbase.enerregions.cinitbase.gridsectioning.cinitbase.tools.cinit.cdefcoord=9;|defcoord=JET3COORDS;a=0.9375;|a=0.92;R0=0.0;|Rout=40.0;|> R0 =0

Definition at line 258 of file newcodediff.txt.

Rout =10

Definition at line 259 of file newcodediff.txt.

SCANARG

Definition at line 379 of file newcodediff.txt.

tf = 2000.0

Definition at line 263 of file newcodediff.txt.

void int ti[]

Definition at line 378 of file newcodediff.txt.

uh = -geom.gcon[GIND(0,2)]/sqrt(-geom.gcon[GIND(0,0)])

Definition at line 296 of file newcodediff.txt.

up = -geom.gcon[GIND(0,3)]/sqrt(-geom.gcon[GIND(0,0)])

Definition at line 297 of file newcodediff.txt.

init ns backup h init ns c ur = -geom.gcon[GIND(0,1)]/sqrt(-geom.gcon[GIND(0,0)])

Definition at line 295 of file newcodediff.txt.

static void int FTYPE FTYPE V

Definition at line 367 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int int FTYPE* vec

Definition at line 367 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int FTYPE FTYPE* static vecortho void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int int FTYPE FTYPE* vecB

Definition at line 368 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int int vectorcomponent

Definition at line 368 of file newcodediff.txt.

void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int FTYPE FTYPE* static vecortho void FTYPE FTYPE FTYPE( FTYPE( FTYPE( FTYPE FTYPE FTYPE( int int int FTYPE* vecv

Definition at line 368 of file newcodediff.txt.

jon_interp_filter c jon_interp_interpolationitself c jon_interp_mnewt c jon_interp_newt c kazfulleos c kazfulleos decsglobalprivate h kazfulleos defsglobalprivate h kazfulleos defs h kazfulleos eostablesdefs h kazfulleos global h kazfulleos_set_arrays c kazfulleos superdefs h liaison c liaison_set_arrays c lnsrch c lubksb c ludcmp c main c math_tools c metric c void int int whichcoord

Definition at line 402 of file newcodediff.txt.

jon_interp_filter c jon_interp_interpolationitself c jon_interp_mnewt c jon_interp_newt c kazfulleos c kazfulleos decsglobalprivate h kazfulleos defsglobalprivate h kazfulleos defs h kazfulleos eostablesdefs h kazfulleos global h kazfulleos_set_arrays c kazfulleos superdefs h liaison c liaison_set_arrays c lnsrch c lubksb c ludcmp c main c math_tools c metric c void int int FTYPE * X

Definition at line 367 of file newcodediff.txt.

Generated on Fri May 20 2016 15:52:36 for HARM by   doxygen 1.8.3.1