sandbox/magfield/src
--- sandbox/magfield/src/OneDOneComponentMagneticFieldMap.cpp 2014-04-04 21:16:00 UTC (rev 3086)
+++ sandbox/magfield/src/OneDOneComponentMagneticFieldMap.cpp 2014-04-04 21:51:04 UTC (rev 3087)
@@ -8,7 +8,7 @@
using namespace std;
OneDimOneComponentMagneticFieldMap::OneDimOneComponentMagneticFieldMap( COORDINATE coord, BVAL bval, const char* filename, double offset, bool isSymmetric )
- :_coord(coord),_bval(bval),_offset(offset),_isSymmetric(false)
+ :_coord(coord),_bval(bval),_offset(offset),_isSymmetric(isSymmetric)
{
cout << "\n-----------------------------------------------------------"
<< "\n Magnetic field"
sandbox/magfield/src
--- sandbox/magfield/src/OneDimOneComponentMagneticFieldMap_t.cpp 2014-04-04 21:16:00 UTC (rev 3086)
+++ sandbox/magfield/src/OneDimOneComponentMagneticFieldMap_t.cpp 2014-04-04 21:51:04 UTC (rev 3087)
@@ -46,8 +46,21 @@
cout << p[0] << " " << p[1] << " "<< p[2] << " : "<< B[0] << " "<< B[1] << " "<< B[2] << " " << endl;
if(abs(B[val]-By[i]) >= epsilon ) cout << B[val] << " " << By[i] << endl;
assert( abs(B[val]-By[i]) < epsilon );
+ p[coord] = -z[i];
+ map.GetFieldValue( p, B );
+ cout << p[0] << " " << p[1] << " "<< p[2] << " : "<< B[0] << " "<< B[1] << " "<< B[2] << " " << endl;
+ if(abs(B[val]-By[i]) >= epsilon ) cout << B[val] << " " << By[i] << endl;
+ assert( abs(B[val]-By[i]) < epsilon );
}
+ p[coord] = max + 1.;
+ map.GetFieldValue( p, B );
+ for(int i=0; i<3; ++i) assert(B[i]==0);
+
+ p[coord] = -p[coord];
+ map.GetFieldValue( p, B );
+ for(int i=0; i<3; ++i) assert(B[i]==0);
+
//clean up
delete[] z;
delete[] By;