bc4eMaxwell.h

Go to the documentation of this file.

! ************ Answer *******************
      gIII=-tau1*(c22*uss+c2*us)-tau2*(c22*vss+c2*vs)

      tauU=tau1*u(i1,i2,i3,ex)+tau2*u(i1,i2,i3,ey)

      tauUp1=tau1*u(i1+is1,i2+is2,i3+is3,ex)+tau2*u(i1+is1,i2+is2,i3+is3,ey)

      tauUp2=tau1*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+tau2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)

      gIV=-6*tauU+4*tauUp1-tauUp2 +gIVf

      ttu1=1/4.*(16*c11*tauUp1-30*c11*tauU-c11*tauUp2-c11*gIV+8*c1*dra*tauUp1-c1*dra*tauUp2+c1*dra*gIV-12*gIII*dra**2-12*tau1DotUtt*dra**2)/(-3*c11+c1*dra)
      ttu2=(16*c11*tauUp1-30*c11*tauU-c11*tauUp2-4*c11*gIV+8*c1*dra*tauUp1-c1*dra*tauUp2+2*c1*dra*gIV-12*gIII*dra**2-12*tau1DotUtt*dra**2)/(-3*c11+c1*dra)

      u(i1-is1,i2-is2,i3-is3,ex) = -1.*(-12.*a12*tau1*tau2*Da1DotU*dra+12.*tau2**2*a11*Da1DotU*dra-1.*tau2**2*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)*a11m2+tau2*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)*a12m2*tau1-1.*tau2**2*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)*a11m2+tau2*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)*a12m2*tau1-5.*tau2*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)*a12m2*tau1-1.*a12*tau1*tau2*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)+5.*tau2**2*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)*a11m2+tau2**2*a11*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-5.*tau2*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)*a12m2*tau1+5.*tau2**2*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)*a11m2-1.*a12*tau1*tau2*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+tau2**2*a11*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-8.*a12*tau1*a12m1*ttu1+8.*tau2*a11*a12m1*ttu1-1.*tau2*a11*a12m2*ttu2+5.*a12*ttu1*a12m2*tau1-5.*a12*ttu1*tau2*a11m2+a12*ttu2*tau2*a11m2-10.*tau2*a12*u(i1,i2,i3,ey)*a12m2*tau1+10.*tau2**2*a12*u(i1,i2,i3,ey)*a11m2+10.*tau2**2*a11*u(i1,i2,i3,ex)*a11m2-10.*tau2*a11*u(i1,i2,i3,ex)*a12m2*tau1+8.*a12*tau1*tau2*a11p1*u(i1+is1,i2+is2,i3+is3,ex)-8.*tau2**2*a11*a11p1*u(i1+is1,i2+is2,i3+is3,ex)+10.*tau2*a11*u(i1+is1,i2+is2,i3+is3,ex)*a12m2*tau1-10.*tau2**2*a11*u(i1+is1,i2+is2,i3+is3,ex)*a11m2+8.*a12*tau1*tau2*a12p1*u(i1+is1,i2+is2,i3+is3,ey)-10.*tau2**2*a12*u(i1+is1,i2+is2,i3+is3,ey)*a11m2-8.*tau2**2*a11*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+10.*tau2*a12*u(i1+is1,i2+is2,i3+is3,ey)*a12m2*tau1)/(5.*tau2*a11*a12m2*tau1-5.*tau2**2*a11*a11m2-5.*a12*tau1**2*a12m2+5.*a12*tau1*tau2*a11m2-8.*tau2*a11*a12m1*tau1-8.*a12*tau1*tau2*a11m1+8.*tau2**2*a11*a11m1+8.*a12*tau1**2*a12m1)

      u(i1-is1,i2-is2,i3-is3,ey) = -1.*(8.*tau1*a11m1*a12*ttu1-8.*ttu1*tau2*a11*a11m1+5.*ttu1*tau2*a11*a11m2-1.*tau1*a12*ttu2*a11m2+tau1*a11*a12m2*ttu2-5.*ttu1*a11*a12m2*tau1+10.*a11*u(i1,i2,i3,ex)*a12m2*tau1**2-10.*tau1*tau2*a11*u(i1,i2,i3,ex)*a11m2-10.*a12*u(i1+is1,i2+is2,i3+is3,ey)*a12m2*tau1**2-8.*a12*tau1**2*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+8.*tau1*tau2*a11*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+10.*tau1*tau2*a12*u(i1+is1,i2+is2,i3+is3,ey)*a11m2-10.*a11*u(i1+is1,i2+is2,i3+is3,ex)*a12m2*tau1**2+10.*tau1*tau2*a11*u(i1+is1,i2+is2,i3+is3,ex)*a11m2+8.*tau1*tau2*a11*a11p1*u(i1+is1,i2+is2,i3+is3,ex)-8.*a12*tau1**2*a11p1*u(i1+is1,i2+is2,i3+is3,ex)-10.*tau1*tau2*a12*u(i1,i2,i3,ey)*a11m2+10.*a12*u(i1,i2,i3,ey)*a12m2*tau1**2-5.*tau1*tau2*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)*a11m2-1.*tau1*tau2*a11*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+a12*tau1**2*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+5.*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)*a12m2*tau1**2+a12*tau1**2*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)+5.*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)*a12m2*tau1**2-1.*tau1*tau2*a11*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-5.*tau1*tau2*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)*a11m2-1.*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)*a12m2*tau1**2+tau1*tau2*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)*a11m2+tau1*tau2*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)*a11m2-1.*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)*a12m2*tau1**2+12.*a12*tau1**2*Da1DotU*dra-12.*tau1*tau2*a11*Da1DotU*dra)/(5.*tau2*a11*a12m2*tau1-5.*tau2**2*a11*a11m2-5.*a12*tau1**2*a12m2+5.*a12*tau1*tau2*a11m2-8.*tau2*a11*a12m1*tau1-8.*a12*tau1*tau2*a11m1+8.*tau2**2*a11*a11m1+8.*a12*tau1**2*a12m1)

      u(i1-2*is1,i2-2*is2,i3-2*is3,ex) = -1.*(-40.*tau2*a12m1*tau1*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-40.*tau2*a11m1*a12*ttu1-80.*tau2**2*a11m1*a12*u(i1+is1,i2+is2,i3+is3,ey)-5.*a12*tau1*tau2*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)+5.*tau2**2*a11*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-60.*a12*tau1*tau2*Da1DotU*dra+60.*tau2**2*a11*Da1DotU*dra-8.*tau2**2*a11m1*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)-8.*tau2**2*a11m1*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)+8.*tau2*a12m1*tau1*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)+8.*tau2*a12m1*tau1*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)+40.*tau2**2*a11m1*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+5.*tau2**2*a11*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+40.*tau2**2*a11m1*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-40.*tau2*a12m1*tau1*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-5.*a12*tau1*tau2*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-8.*a12m1*tau1*a12*ttu2+40.*tau2*a11*a12m1*ttu1+8.*tau2*a11m1*a12*ttu2-5.*tau2*a11*a12m2*ttu2+5.*a12*ttu2*a12m2*tau1+80.*tau2**2*a11m1*a12*u(i1,i2,i3,ey)-80.*tau2*a12m1*tau1*a11*u(i1,i2,i3,ex)+80.*tau2**2*a11m1*a11*u(i1,i2,i3,ex)-80.*tau2**2*a11m1*a11*u(i1+is1,i2+is2,i3+is3,ex)-80.*tau2*a12m1*tau1*a12*u(i1,i2,i3,ey)+40.*a12*tau1*tau2*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+80.*tau2*a12m1*tau1*a12*u(i1+is1,i2+is2,i3+is3,ey)-40.*tau2**2*a11*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+80.*tau2*a12m1*tau1*a11*u(i1+is1,i2+is2,i3+is3,ex)+40.*a12*tau1*tau2*a11p1*u(i1+is1,i2+is2,i3+is3,ex)-40.*tau2**2*a11*a11p1*u(i1+is1,i2+is2,i3+is3,ex))/(5.*tau2*a11*a12m2*tau1-5.*tau2**2*a11*a11m2-5.*a12*tau1**2*a12m2+5.*a12*tau1*tau2*a11m2-8.*tau2*a11*a12m1*tau1-8.*a12*tau1*tau2*a11m1+8.*tau2**2*a11*a11m1+8.*a12*tau1**2*a12m1)

      u(i1-2*is1,i2-2*is2,i3-2*is3,ey) = (40.*tau1*a11*a12m1*ttu1-8.*ttu2*a11*a12m1*tau1-40.*tau1*a11m1*a12*ttu1+8.*ttu2*tau2*a11*a11m1+5.*tau1*a12*ttu2*a11m2-5.*ttu2*tau2*a11*a11m2-40.*tau1*tau2*a11*a12p1*u(i1+is1,i2+is2,i3+is3,ey)+40.*a12*tau1**2*a12p1*u(i1+is1,i2+is2,i3+is3,ey)-80.*tau1*tau2*a11m1*a12*u(i1+is1,i2+is2,i3+is3,ey)-80.*tau1*tau2*a11m1*a11*u(i1+is1,i2+is2,i3+is3,ex)-40.*tau1*tau2*a11*a11p1*u(i1+is1,i2+is2,i3+is3,ex)+80.*a12m1*tau1**2*a12*u(i1+is1,i2+is2,i3+is3,ey)+80.*a12m1*tau1**2*a11*u(i1+is1,i2+is2,i3+is3,ex)+40.*a12*tau1**2*a11p1*u(i1+is1,i2+is2,i3+is3,ex)+80.*tau1*tau2*a11m1*a12*u(i1,i2,i3,ey)-80.*a12m1*tau1**2*a12*u(i1,i2,i3,ey)+80.*tau1*tau2*a11m1*a11*u(i1,i2,i3,ex)-80.*a12m1*tau1**2*a11*u(i1,i2,i3,ex)+5.*tau1*tau2*a11*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-5.*a12*tau1**2*a11p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)+40.*tau1*tau2*a11m1*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-40.*a12m1*tau1**2*a11*u(i1+2*is1,i2+2*is2,i3+2*is3,ex)-60.*a12*tau1**2*Da1DotU*dra+60.*tau1*tau2*a11*Da1DotU*dra+40.*tau1*tau2*a11m1*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)+5.*tau1*tau2*a11*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-5.*a12*tau1**2*a12p2*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-40.*a12m1*tau1**2*a12*u(i1+2*is1,i2+2*is2,i3+2*is3,ey)-8.*tau1*tau2*a11m1*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)+8.*a12m1*tau1**2*a12*u(i1+3*is1,i2+3*is2,i3+3*is3,ey)-8.*tau1*tau2*a11m1*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex)+8.*a12m1*tau1**2*a11*u(i1+3*is1,i2+3*is2,i3+3*is3,ex))/(5.*tau2*a11*a12m2*tau1-5.*tau2**2*a11*a11m2-5.*a12*tau1**2*a12m2+5.*a12*tau1*tau2*a11m2-8.*tau2*a11*a12m1*tau1-8.*a12*tau1*tau2*a11m1+8.*tau2**2*a11*a11m1+8.*a12*tau1**2*a12m1)


 ! *********** done *********************