bug fix in current integration

pull/1/head
Thorsten Liebig 2010-05-07 22:58:26 +02:00
parent 5b3ccd3f53
commit 7cc2a7fcf2
1 changed files with 6 additions and 6 deletions

View File

@ -54,27 +54,27 @@ int ProcessCurrent::Process()
}
//x-current
if (m_start_inside[1] && m_start_inside[2])
for (unsigned int i=start[0];i<stop[0];++i)
for (unsigned int i=start[0]+1;i<=stop[0];++i)
current+=Eng->GetCurr(0,i,start[1],start[2]);
//y-current
if (m_stop_inside[0] && m_start_inside[2])
for (unsigned int i=start[1];i<stop[1];++i)
for (unsigned int i=start[1]+1;i<=stop[1];++i)
current+=Eng->GetCurr(1,stop[0],i,start[2]);
//z-current
if (m_stop_inside[0] && m_stop_inside[1])
for (unsigned int i=start[2];i<stop[2];++i)
for (unsigned int i=start[2]+1;i<=stop[2];++i)
current+=Eng->GetCurr(2,stop[0],stop[1],i);
//x-current
if (m_stop_inside[1] && m_stop_inside[2])
for (unsigned int i=start[0];i<stop[0];++i)
for (unsigned int i=start[0]+1;i<=stop[0];++i)
current-=Eng->GetCurr(0,i,stop[1],stop[2]);
//y-current
if (m_start_inside[0] && m_stop_inside[2])
for (unsigned int i=start[1];i<stop[1];++i)
for (unsigned int i=start[1]+1;i<=stop[1];++i)
current-=Eng->GetCurr(1,start[0],i,stop[2]);
//z-current
if (m_start_inside[0] && m_start_inside[1])
for (unsigned int i=start[2];i<stop[2];++i)
for (unsigned int i=start[2]+1;i<=stop[2];++i)
current-=Eng->GetCurr(2,start[0],start[1],i);
// cerr << "ts: " << Eng->numTS << " i: " << current << endl;