global.stepch.h

Go to the documentation of this file.

#if(MAXTIMEORDER==4)
#define UFSET(theCUf,thedt,theUi,theUf,thedUriemann,thedUgeom,thedUnongeom) ((theCUf[0])*(theUi) + (theCUf[1])*(theUf) + (theCUf[2])*(thedt)*((thedUriemann)+(thedUgeom))  + (theCUf[NUMPREDTCUFS+0])*(thedt)*(thedUnongeom[0]) + (theCUf[NUMPREDTCUFS+1])*(thedt)*(thedUnongeom[1]) + (theCUf[NUMPREDTCUFS+2])*(thedt)*(thedUnongeom[2]) + (theCUf[NUMPREDTCUFS+3])*(thedt)*(thedUnongeom[3]) )
#define UCUMUPDATE(theCunew,thedt,theUi,theUf,thedUriemann,thedUgeom,thedUnongeom) ((theCunew[0])*(theUi) + (theCunew[1])*(thedt)*((thedUriemann)+(thedUgeom)) + (theCunew[2])*(theUf) + (theCunew[NUMPREDTCUFS+0])*(thedt)*(thedUnongeom[0]) + (theCunew[NUMPREDTCUFS+1])*(thedt)*(thedUnongeom[1]) + (theCunew[NUMPREDTCUFS+2])*(thedt)*(thedUnongeom[2]) + (theCunew[NUMPREDTCUFS+3])*(thedt)*(thedUnongeom[3]) )
#elif(MAXTIMEORDER==5)
#define UFSET(theCUf,thedt,theUi,theUf,thedUriemann,thedUgeom,thedUnongeom) ((theCUf[0])*(theUi) + (theCUf[1])*(theUf) + (theCUf[2])*(thedt)*((thedUriemann)+(thedUgeom))  + (theCUf[NUMPREDTCUFS+0])*(thedt)*(thedUnongeom[0]) + (theCUf[NUMPREDTCUFS+1])*(thedt)*(thedUnongeom[1]) + (theCUf[NUMPREDTCUFS+2])*(thedt)*(thedUnongeom[2]) + (theCUf[NUMPREDTCUFS+3])*(thedt)*(thedUnongeom[3]) + (theCUf[NUMPREDTCUFS+4])*(thedt)*(thedUnongeom[4]) )
#define UCUMUPDATE(theCunew,thedt,theUi,theUf,thedUriemann,thedUgeom,thedUnongeom) ((theCunew[0])*(theUi) + (theCunew[1])*(thedt)*((thedUriemann)+(thedUgeom)) + (theCunew[2])*(theUf) + (theCunew[NUMPREDTCUFS+0])*(thedt)*(thedUnongeom[0]) + (theCunew[NUMPREDTCUFS+1])*(thedt)*(thedUnongeom[1]) + (theCunew[NUMPREDTCUFS+2])*(thedt)*(thedUnongeom[2]) + (theCunew[NUMPREDTCUFS+3])*(thedt)*(thedUnongeom[3]) + (theCunew[NUMPREDTCUFS+4])*(thedt)*(thedUnongeom[4]) )
#else
#error "Need to setup UFSET for other MAXTIMEORDER"
#endif


#define dUfromUFSET(theCUf,thedt,theUi,theUf,theUfnew) ( ((theUfnew) - ((theCUf[0])*(theUi) + (theCUf[1])*(theUf)))/((theCUf[2])*(thedt)) )