zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix in Operator_Cylinder

pull/1/head
Thorsten Liebig 2010-04-14 11:01:17 +02:00
parent 1d54bda1d0
commit b3e2d9428a
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
FDTD/operator_cylinder.cpp
View File

@ -136,10 +136,13 @@ int Operator_Cylinder::CalcECOperator()
for (pos[1]=0;pos[1]<numLines[1]-CC_closedAlpha;++pos[1]) for (pos[1]=0;pos[1]<numLines[1]-CC_closedAlpha;++pos[1])
{ {
Calc_ECPos(2,pos,inEC); Calc_ECPos(2,pos,inEC);
C+=inEC[0]; // if (pos[2]==0)
G+=inEC[1]; // cerr << inEC[0] << endl;
C+=inEC[0]*0.5;
G+=inEC[1]*0.5;
} }
// cerr << C << " and " << G << endl; if (pos[2]==0)
cerr << C << " and " << G << endl;
vv_R0[pos[2]] = (1-dT*G/2/C)/(1+dT*G/2/C); vv_R0[pos[2]] = (1-dT*G/2/C)/(1+dT*G/2/C);
vi_R0[pos[2]] = (dT/C)/(1+dT*G/2/C); vi_R0[pos[2]] = (dT/C)/(1+dT*G/2/C);
} }