CG  Version 25
bcdivMaxwell.h
Go to the documentation of this file.
1  b3u=c11*a11
2  b3v=a12*c11
3  b2u=(2*c22*a12s*a11**2+c1*a11**2*a22+c11r*a11**2*a22-2*c22*a11s*a12*a11-c1*a11*a21*a12-c11r*a11*a21*a12-a11r*c11*a21*a12+a11r*c11*a11*a22)/(-a21*a12+a11*a22)
4  b2v=(-2*c22*a11s*a12**2+2*c22*a12s*a11*a12+a12*c1*a11*a22+a12*c11r*a11*a22-a12r*c11*a21*a12+a12r*c11*a11*a22-c11r*a21*a12**2-c1*a21*a12**2)/(-a21*a12+a11*a22)
5  b1u=(-4*c22*a11s*a12*a11r-c22*a11ss*a11*a22+c1r*a11**2*a22+4*c22*a12s*a11*a11r-c1r*a11*a21*a12+2*c22*a11s**2*a22-2*c22*a12s*a21*a11s-2*c22*a11s*a12*a21s+2*c22*a12s*a11*a21s-c2*a11*a11s*a22+c22*a11ss*a21*a12+a12*c2*a21*a11s-a11r*c1*a21*a12+a11r*c1*a11*a22)/(-a21*a12+a11*a22)
6  b1v=(-2*c22*a21*a12s**2+2*c22*a11s*a12s*a22+2*c22*a12s*a11*a22s-2*c22*a11s*a12*a22s+c22*a12ss*a21*a12+4*c22*a12s*a11*a12r-c22*a12ss*a11*a22-4*c22*a11s*a12*a12r+a12*c1r*a11*a22+a12*c2*a21*a12s-a12r*c1*a21*a12+a12r*c1*a11*a22-c1r*a21*a12**2-c2*a11*a12s*a22)/(-a21*a12+a11*a22)
7  bf =-(-2*c22*a12s*a11*a21rs*uex+2*c22*a12s*a21*a21ss*uex+3*c22*a22s*vss*a11*a22-c22r*uss*a11**2*a22-2*c22*a11s*a22**2*vss-c22*a21**2*usss*a12+c22*a22**2*vsss*a11+2*c22*a12s*a21**2*uss-c2r*us*a11**2*a22+c2*a11*a22**2*vss-3*c22*a21ss*us*a21*a12-2*c22*a12s*a11*a11rr*uex+3*c22*a21ss*us*a11*a22-3*c22*a22ss*vs*a21*a12+2*c22*a11s*a12*a22r*vs+3*c22*a22ss*vs*a11*a22+c2r*us*a11*a21*a12+2*c2*a11*a22s*vs*a22-c22*a22*vsss*a21*a12-3*c22*a21s*uss*a21*a12+3*c22*a21s*uss*a11*a22-c22*a12rss*uey*a21*a12+c22*a12rss*uey*a11*a22-c22*a11rss*uex*a21*a12+2*c22*a12s*a21*a22ss*uey+2*c22*a12s*a21*a22*vss+4*c22*a12s*a21*a21s*us+c22*a11rss*uex*a11*a22-2*c22*a12rs*vs*a21*a12+2*c22*a12rs*vs*a11*a22-3*c22*a22s*vss*a21*a12-c22*a21sss*uex*a21*a12+c22*a21sss*uex*a11*a22-2*c22*a11rs*us*a21*a12+2*c22*a12s*a21*a11rs*uex+2*c22*a12s*a21*a11r*us+2*c22*a11s*a12*a21r*us+2*c22*a11s*a12*a21rs*uex-2*c22*a12s*a11*a21r*us-2*c22*a12s*a11*a22rs*uey-2*c22*a12s*a11*a22r*vs+2*c22*a12s*a21*a12r*vs+2*c22*a11s*a12*a22rs*uey+2*c22*a11rs*us*a11*a22+c2*a11*a21*uss*a22+c2*a11*a11rs*uex*a22-2*c22*a11s*a12rs*uey*a22+c22*a22sss*uey*a11*a22-2*c22*a12s*a11*a12rr*uey+4*c22*a12s*a21*a22s*vs+2*c22*a12s*a21*a12rs*uey-2*c22*a11s*a21ss*uex*a22-2*c22*a11s*a22ss*uey*a22+2*c22*a11s*a12*a11rr*uex-4*c22*a11s*a22s*vs*a22-2*c22*a11s*a11rs*uex*a22-2*c22*a11s*a21*uss*a22-4*c22*a11s*a21s*us*a22-2*c22*a11s*a11r*us*a22+2*c22*a11s*a12*a12rr*uey+c22*a21*usss*a11*a22-c22*a22sss*uey*a21*a12-2*c22*a11s*a12r*vs*a22+c22r*uss*a11*a21*a12-a12*c22r*vss*a11*a22-a12*c2r*vs*a11*a22-a12*c2*a21*a12rs*uey-2*a12*c2*a21*a22s*vs-a12*c2*a21*a11rs*uex-a12*c2*a21*a22*vss-2*a12*c2*a21*a21s*us-a12*c2*a21*a21ss*uex-a12*c2*a21*a22ss*uey-a12*c2*a21**2*uss+c2r*vs*a21*a12**2+c22r*vss*a21*a12**2+c2*a11*a12rs*uey*a22+c2*a11*a22ss*uey*a22+c2*a11*a21ss*uex*a22+2*c2*a11*a21s*us*a22)/(-a21*a12+a11*a22)
8 
9 ! -- Here are the approximations for urs, vrs from the divergence
10 ! ursm =-(2*a22s*vs*a22-a12*a11*urr-a12*a22r*vs-a12*a22s*vr-a12*a22rs*uey-a12*a21r*us-a12*a21s*ur-a12*a21rs*uex-2*a12*a11r*ur-a12*a11rr*uex+a12rs*uey*a22+a12r*vs*a22+a12s*vr*a22+a22ss*uey*a22+a21ss*uex*a22+2*a21s*us*a22+a21*uss*a22+a22**2*vss-a12**2*vrr+a11rs*uex*a22+a11r*us*a22-2*a12*a12r*vr+a11s*ur*a22-a12*a12rr*uey)/(-a21*a12+a11*a22)
11 ! vrsm =(-a11**2*urr-2*a11*a12r*vr-a11*a12rr*uey+a21*a12r*vs+a21*a12rs*uey+2*a21*a22s*vs+a21*a11s*ur+a21*a11r*us+a21*a11rs*uex+a21*a12s*vr-a11*a22r*vs-a11*a22s*vr-a11*a22rs*uey-a11*a21r*us-a11*a21s*ur-a11*a21rs*uex-a11*a12*vrr+a21*a22*vss+2*a21*a21s*us+a21*a21ss*uex+a21*a22ss*uey-2*a11*a11r*ur+a21**2*uss-a11*a11rr*uex)/(-a21*a12+a11*a22)