From 36436470a411db1a38e96afc4bca05cb43d31810 Mon Sep 17 00:00:00 2001
From: Marina Ganeva <m.ganeva@fz-juelich.de>
Date: Tue, 26 Nov 2013 12:30:33 +0100
Subject: [PATCH] PolarizedDWBAZeroMag works also for python now.
---
.../FunctionalTests/TestPyCore/CMakeLists.txt | 1 +
.../FunctionalTests/TestPyCore/TestPyCore.py | 4 +-
.../TestPyCore/polarizeddwbazeromag.py | 104 ++++++++++++++++++
3 files changed, 108 insertions(+), 1 deletion(-)
create mode 100644 Tests/FunctionalTests/TestPyCore/polarizeddwbazeromag.py
diff --git a/Tests/FunctionalTests/TestPyCore/CMakeLists.txt b/Tests/FunctionalTests/TestPyCore/CMakeLists.txt
index f288b9d9d1c..860daac8990 100644
--- a/Tests/FunctionalTests/TestPyCore/CMakeLists.txt
+++ b/Tests/FunctionalTests/TestPyCore/CMakeLists.txt
@@ -16,6 +16,7 @@ set(list_of_tests
"isgisaxs11.py"
"isgisaxs15.py"
# "mesocrystal1.py"
+ "polarizeddwbazeromag.py"
)
foreach(_test ${list_of_tests})
diff --git a/Tests/FunctionalTests/TestPyCore/TestPyCore.py b/Tests/FunctionalTests/TestPyCore/TestPyCore.py
index 94afa627190..8cea65a1646 100644
--- a/Tests/FunctionalTests/TestPyCore/TestPyCore.py
+++ b/Tests/FunctionalTests/TestPyCore/TestPyCore.py
@@ -18,6 +18,7 @@ import isgisaxs10
import isgisaxs11
import isgisaxs15
import mesocrystal1
+import polarizeddwbazeromag
Tests = {
"IsGISAXS01": isgisaxs01.runTest,
@@ -31,7 +32,8 @@ Tests = {
"IsGISAXS10": isgisaxs10.runTest,
"IsGISAXS11": isgisaxs11.runTest,
"IsGISAXS15": isgisaxs15.runTest,
- "MesoCrystal1": mesocrystal1.runTest,
+# "MesoCrystal1": mesocrystal1.runTest,
+ "PolarizedDWBAZeroMag": polarizeddwbazeromag.runTest,
}
test_info = []
diff --git a/Tests/FunctionalTests/TestPyCore/polarizeddwbazeromag.py b/Tests/FunctionalTests/TestPyCore/polarizeddwbazeromag.py
new file mode 100644
index 00000000000..432028afb4c
--- /dev/null
+++ b/Tests/FunctionalTests/TestPyCore/polarizeddwbazeromag.py
@@ -0,0 +1,104 @@
+# IsGISAXS01 example: Mixture of cylinders and prisms without interference
+import sys
+import os
+import numpy
+import gzip
+
+sys.path.append(os.path.abspath(
+ os.path.join(os.path.split(__file__)[0],
+ '..', '..', '..', 'lib')))
+
+
+from libBornAgainCore import *
+
+
+# ----------------------------------
+# describe sample and run simulation
+# ----------------------------------
+def RunSimulation():
+ # defining materials
+ mAmbience = MaterialManager.getHomogeneousMaterial("Air", 0.0, 0.0 )
+ mSubstrate = MaterialManager.getHomogeneousMaterial("Substrate", 6e-6, 2e-8 )
+
+ magnetic_field = kvector_t(0,0,0)
+
+ magParticle = MaterialManager.getHomogeneousMagneticMaterial("magParticle", 6e-4, 2e-8, magnetic_field )
+ # collection of particles
+ cylinder_ff = FormFactorCylinder(5*nanometer, 5*nanometer)
+ cylinder = Particle(magParticle, cylinder_ff)
+
+ particle_decoration = ParticleDecoration()
+ particle_decoration.addParticle(cylinder, 0.0, 1.0)
+ interference = InterferenceFunctionNone()
+ particle_decoration.addInterferenceFunction(interference)
+
+ # air layer with particles and substrate form multi layer
+ air_layer = Layer(mAmbience)
+ air_layer.setDecoration(particle_decoration)
+ substrate_layer = Layer(mSubstrate, 0)
+ multi_layer = MultiLayer()
+ multi_layer.addLayer(air_layer)
+ multi_layer.addLayer(substrate_layer)
+
+ # build and run experiment
+ simulation = Simulation()
+ simulation.setDetectorParameters(100,0*degree, 2.0*degree, 100, 0.0*degree, 2.0*degree, True)
+ simulation.setBeamParameters(1.0*angstrom, 0.2*degree, 0.0*degree)
+ simulation.setSample(multi_layer)
+ simulation.runSimulation()
+ ## intensity data
+ return simulation.getIntensityData().getArray()
+
+# ----------------------------------
+# read reference data from file
+# ----------------------------------
+def GetReferenceData():
+ path = os.path.split(__file__)[0]
+ if path: path +="/"
+ f = gzip.open(path+'../../ReferenceData/BornAgain/isgi_cylinder_DWBA.ima.gz', 'rb')
+ reference=numpy.fromstring(f.read(),numpy.float64,sep=' ')
+ f.close()
+ return reference
+
+
+# --------------------------------------------------------------
+# calculate numeric difference between result and reference data
+# --------------------------------------------------------------
+def GetDifference(data, reference):
+ reference = reference.reshape(data.shape)
+ # calculating relative average difference
+ data -= reference
+ diff=0.0
+ epsilon = sys.float_info.epsilon
+ for x, y in numpy.ndindex(data.shape):
+ v1 = data[x][y]
+ v2 = reference[x][y]
+ if v1 <= epsilon and v2 <= epsilon:
+ diff += 0.0
+ elif(v2 <= epsilon):
+ diff += abs(v1/epsilon)
+ else:
+ diff += abs(v1/v2)
+ return diff/data.size
+
+
+# --------------------------------------------------------------
+# run test and analyse test results
+# --------------------------------------------------------------
+def runTest():
+ result = RunSimulation()
+ reference = GetReferenceData()
+
+ diff = GetDifference(result, reference)
+ status = "OK"
+ if(diff > 2e-10 or numpy.isnan(diff)): status = "FAILED"
+ return "PolarizedDWBAZeroMag", "cylinder DWBA form factor with matrix calculation", status
+
+
+#-------------------------------------------------------------
+# main()
+#-------------------------------------------------------------
+if __name__ == '__main__':
+ name,description,status = runTest()
+ print name,description,status
+ if("FAILED" in status) : exit(1)
--
GitLab