00017 {
00018 double x,y,z,bmf,emf;
00019 x = r.x()/m;
00020 y = r.y()/m;
00021 z = r.z()/m;
00022
00023 if(part==1)
00024 {
00025 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx0[i]; }
00026 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx0[i]; }
00027 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx1[i]; }
00028 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx1[i]; }
00029 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx2[i]; }
00030 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx2[i]; }
00031 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx3[i]; }
00032 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx3[i]; }
00033 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx4[i]; }
00034 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx4[i]; }
00035 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx5[i]; }
00036 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx5[i]; }
00037 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx6[i]; }
00038 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx6[i]; }
00039 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx7[i]; }
00040 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx7[i]; }
00041 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx8[i]; }
00042 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z
00043 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y
00044 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y;
00045 b[0] = bmf*tesla;
00046
00047 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy0[i]; }
00048 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy0[i]; }
00049 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy1[i]; }
00050 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy1[i]; }
00051 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy2[i]; }
00052 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy2[i]; }
00053 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy3[i]; }
00054 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy3[i]; }
00055 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy4[i]; }
00056 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy4[i]; }
00057 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy5[i]; }
00058 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy5[i]; }
00059 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy6[i]; }
00060 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy6[i]; }
00061 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy7[i]; }
00062 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy7[i]; }
00063 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy8[i]; }
00064 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z
00065 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y
00066 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y;
00067 b[1] = bmf*tesla;
00068
00069 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz0[i]; }
00070 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz0[i]; }
00071 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz1[i]; }
00072 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz1[i]; }
00073 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz2[i]; }
00074 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz2[i]; }
00075 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz3[i]; }
00076 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz3[i]; }
00077 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz4[i]; }
00078 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz4[i]; }
00079 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz5[i]; }
00080 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz5[i]; }
00081 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz6[i]; }
00082 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz6[i]; }
00083 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz7[i]; }
00084 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz7[i]; }
00085 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz8[i]; }
00086 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z
00087 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y
00088 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y;
00089 b[2] = bmf*tesla;
00090 }
00091
00092 if(part==0)
00093 {
00094 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx0[i]; }
00095 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx0[i]; }
00096 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx1[i]; }
00097 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx1[i]; }
00098 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx2[i]; }
00099 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx2[i]; }
00100 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx3[i]; }
00101 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx3[i]; }
00102 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx4[i]; }
00103 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx4[i]; }
00104 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx5[i]; }
00105 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx5[i]; }
00106 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx6[i]; }
00107 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx6[i]; }
00108 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx7[i]; }
00109 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx7[i]; }
00110 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx8[i]; }
00111 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x
00112 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y
00113 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x;
00114 b[0] = emf*tesla;
00115
00116 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy0[i]; }
00117 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy0[i]; }
00118 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy1[i]; }
00119 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy1[i]; }
00120 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy2[i]; }
00121 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy2[i]; }
00122 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy3[i]; }
00123 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy3[i]; }
00124 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy4[i]; }
00125 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy4[i]; }
00126 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy5[i]; }
00127 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy5[i]; }
00128 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy6[i]; }
00129 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy6[i]; }
00130 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy7[i]; }
00131 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy7[i]; }
00132 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy8[i]; }
00133 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x
00134 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y
00135 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x;
00136 b[1] = emf*tesla;
00137
00138 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz0[i]; }
00139 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz0[i]; }
00140 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz1[i]; }
00141 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz1[i]; }
00142 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz2[i]; }
00143 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz2[i]; }
00144 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz3[i]; }
00145 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz3[i]; }
00146 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz4[i]; }
00147 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz4[i]; }
00148 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz5[i]; }
00149 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz5[i]; }
00150 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz6[i]; }
00151 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz6[i]; }
00152 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz7[i]; }
00153 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz7[i]; }
00154 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz8[i]; }
00155 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x
00156 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y
00157 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x;
00158 b[2] = emf*tesla;
00159 }
00160 }