--- 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/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;