Tests for polarization functions.

Part of #2795

Change-Id: I30d0c67bac3e16412ca2afaa12ea7e28e8c01f6e

diff --git a/src/gromacs/listed_forces/tests/bonded.cpp b/src/gromacs/listed_forces/tests/bonded.cpp
index 24723bd574705b76aef9f9e3cc777593f635108c..fc8a0528d7900bb9fbc95b1274cfdbc50052000d 100644 (file)
--- a/src/gromacs/listed_forces/tests/bonded.cpp
+++ b/src/gromacs/listed_forces/tests/bonded.cpp
@@ -54,6 +54,7 @@
 #include "gromacs/math/units.h"
 #include "gromacs/math/vec.h"
 #include "gromacs/math/vectypes.h"
+#include "gromacs/mdtypes/mdatom.h"
 #include "gromacs/pbcutil/ishift.h"
 #include "gromacs/pbcutil/pbc.h"
 #include "gromacs/topology/idef.h"
@@ -407,6 +408,85 @@ struct iListInput
         {
             return setRbDihedrals(rbc, rbc);
         }
+        /*! \brief Set parameters for Polarization
+         *
+         * \param[in] alpha Polarizability
+         * \return The structure itself.
+         */
+        iListInput setPolarization(real alpha)
+        {
+            ftype                  = F_POLARIZATION;
+            fep                    = false;
+            iparams.polarize.alpha = alpha;
+            return *this;
+        }
+        /*! \brief Set parameters for Anharmonic Polarization
+         *
+         * \param[in] alpha Polarizability (nm^3)
+         * \param[in] drcut The cut-off distance (nm) after which the
+         *                  fourth power kicks in
+         * \param[in] khyp  The force constant for the fourth power
+         * \return The structure itself.
+         */
+        iListInput setAnharmPolarization(real alpha,
+                                         real drcut,
+                                         real khyp)
+        {
+            ftype                         = F_ANHARM_POL;
+            fep                           = false;
+            iparams.anharm_polarize.alpha = alpha;
+            iparams.anharm_polarize.drcut = drcut;
+            iparams.anharm_polarize.khyp  = khyp;
+            return *this;
+        }
+        /*! \brief Set parameters for Thole Polarization
+         *
+         * \param[in] a      Thole factor
+         * \param[in] alpha1 Polarizability 1 (nm^3)
+         * \param[in] alpha2 Polarizability 2 (nm^3)
+         * \param[in] rfac   Distance factor
+         * \return The structure itself.
+         */
+        iListInput setTholePolarization(real a,
+                                        real alpha1,
+                                        real alpha2,
+                                        real rfac)
+        {
+            ftype                = F_THOLE_POL;
+            fep                  = false;
+            iparams.thole.a      = a;
+            iparams.thole.alpha1 = alpha1;
+            iparams.thole.alpha2 = alpha2;
+            iparams.thole.rfac   = rfac;
+            return *this;
+        }
+        /*! \brief Set parameters for Water Polarization
+         *
+         * \param[in] alpha_x Polarizability X (nm^3)
+         * \param[in] alpha_y Polarizability Y (nm^3)
+         * \param[in] alpha_z Polarizability Z (nm^3)
+         * \param[in] rOH     Oxygen-Hydrogen distance
+         * \param[in] rHH     Hydrogen-Hydrogen distance
+         * \param[in] rOD     Oxygen-Dummy distance
+         * \return The structure itself.
+         */
+        iListInput setWaterPolarization(real alpha_x,
+                                        real alpha_y,
+                                        real alpha_z,
+                                        real rOH,
+                                        real rHH,
+                                        real rOD)
+        {
+            ftype             = F_WATER_POL;
+            fep               = false;
+            iparams.wpol.al_x = alpha_x;
+            iparams.wpol.al_y = alpha_y;
+            iparams.wpol.al_z = alpha_z;
+            iparams.wpol.rOH  = rOH;
+            iparams.wpol.rHH  = rHH;
+            iparams.wpol.rOD  = rOD;
+            return *this;
+        }
 };
 
 /*! \brief Utility to fill iatoms struct
@@ -419,7 +499,8 @@ void fillIatoms(int ftype, std::vector<t_iatom> *iatoms)
     std::unordered_map<int, std::vector<int> > ia =
     { { 2, { 0, 0, 1, 0, 1, 2, 0, 2, 3 } },
       { 3, { 0, 0, 1, 2, 0, 1, 2, 3 } },
-      { 4, { 0, 0, 1, 2, 3 } }};
+      { 4, { 0, 0, 1, 2, 3 } },
+      { 5, { 0, 0, 1, 2, 3, 0 } } };
     EXPECT_TRUE(ftype >= 0 && ftype < F_NRE);
     int nral = interaction_function[ftype].nratoms;
     for (auto &i : ia[nral])
@@ -460,8 +541,11 @@ class ListedForcesTest : public ::testing::TestWithParam<std::tuple<iListInput,
                           const real                  lambda)
         {
             SCOPED_TRACE(std::string("Testing PBC ") + epbc_names[epbc_]);
-            std::vector<int> ddgatindex = { 0, 1, 2, 3 };
-            OutputQuantities output;
+            std::vector<int>  ddgatindex = { 0, 1, 2, 3 };
+            std::vector<real> chargeA    = { 1.5, -2.0, 1.5, -1.0 };
+            t_mdatoms         mdatoms    = {0};
+            mdatoms.chargeA              = chargeA.data();
+            OutputQuantities  output;
             output.energy = bondedFunction(input_.ftype) (iatoms.size(),
                                                           iatoms.data(),
                                                           &input_.iparams,
@@ -470,7 +554,7 @@ class ListedForcesTest : public ::testing::TestWithParam<std::tuple<iListInput,
                                                           &pbc_,
                                                           /* const struct t_graph *g */ nullptr,
                                                           lambda, &output.dvdlambda,
-                                                          /* const struct t_mdatoms *md */ nullptr,
+                                                          &mdatoms,
                                                           /* struct t_fcdata *fcd */ nullptr,
                                                           ddgatindex.data());
             // Internal consistency test of both test input
@@ -559,6 +643,15 @@ std::vector<iListInput> c_InputDihs =
     { iListInput(3e-4, 1e-8).setRbDihedrals(rbc) }
 };
 
+//! Function types for testing polarization. Add new terms at the end.
+std::vector<iListInput> c_InputPols =
+{
+    { iListInput(2e-5, 1e-8).setPolarization(0.12) },
+    { iListInput(1.7e-3, 1e-8).setAnharmPolarization(0.0013, 0.02, 1235.6) },
+    { iListInput(1.4e-3, 1e-8).setTholePolarization(0.26, 0.07, 0.09, 1.6) },
+    { iListInput().setWaterPolarization(0.001, 0.0012, 0.0016, 0.095, 0.15, 0.02) },
+};
+
 //! Coordinates for testing
 std::vector<std::vector<gmx::RVec> > c_coordinatesForTests =
 {
@@ -576,6 +669,8 @@ INSTANTIATE_TEST_CASE_P(Angle, ListedForcesTest, ::testing::Combine(::testing::V
 
 INSTANTIATE_TEST_CASE_P(Dihedral, ListedForcesTest, ::testing::Combine(::testing::ValuesIn(c_InputDihs), ::testing::ValuesIn(c_coordinatesForTests), ::testing::ValuesIn(c_pbcForTests)));
 
+INSTANTIATE_TEST_CASE_P(Polarize, ListedForcesTest, ::testing::Combine(::testing::ValuesIn(c_InputPols), ::testing::ValuesIn(c_coordinatesForTests), ::testing::ValuesIn(c_pbcForTests)));
+
 }  // namespace
 
 }  // namespace gmx

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_0.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_0.xml
new file mode 100644 (file)
index 0000000..1a412cb
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_0.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">117.28019586833446</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">926.23638559389656</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">13.025199172414171</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-1186.74036904218</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-19.971972064368394</Real>
+          <Real Name="Y">115.77954819923707</Real>
+          <Real Name="Z">144.72443524904634</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">6.9467728919542244</Real>
+          <Real Name="Y">-115.77954819923707</Real>
+          <Real Name="Z">115.77954819923707</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_1.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_1.xml
new file mode 100644 (file)
index 0000000..1a412cb
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_1.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">117.28019586833446</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">926.23638559389656</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">13.025199172414171</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-1186.74036904218</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-19.971972064368394</Real>
+          <Real Name="Y">115.77954819923707</Real>
+          <Real Name="Z">144.72443524904634</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">6.9467728919542244</Real>
+          <Real Name="Y">-115.77954819923707</Real>
+          <Real Name="Z">115.77954819923707</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_10.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_10.xml
new file mode 100644 (file)
index 0000000..1a52d5b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_10.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">10827.481942307291</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">85527.567362390473</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1202.4530379358848</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-109576.62812110817</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1844.0698251282447</Real>
+          <Real Name="Y">10693.613119872665</Real>
+          <Real Name="Z">13355.447638845031</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">641.61678719235988</Real>
+          <Real Name="Y">-10693.613119872665</Real>
+          <Real Name="Z">10693.613119872665</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_11.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_11.xml
new file mode 100644 (file)
index 0000000..1a52d5b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_11.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">10827.481942307291</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">85527.567362390473</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1202.4530379358848</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-109576.62812110817</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1844.0698251282447</Real>
+          <Real Name="Y">10693.613119872665</Real>
+          <Real Name="Z">13355.447638845031</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">641.61678719235988</Real>
+          <Real Name="Y">-10693.613119872665</Real>
+          <Real Name="Z">10693.613119872665</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_12.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_12.xml
new file mode 100644 (file)
index 0000000..f9e3c86
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_12.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">6634.6501936621835</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">64134.915916992832</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">16834.287828835641</Real>
+          <Real Name="Z">-47300.628088157187</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-5074.6076174669543</Real>
+          <Real Name="Z">-16834.287828835641</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-11759.680211368686</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_13.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_13.xml
new file mode 100644 (file)
index 0000000..f9e3c86
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_13.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">6634.6501936621835</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">64134.915916992832</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">16834.287828835641</Real>
+          <Real Name="Z">-47300.628088157187</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-5074.6076174669543</Real>
+          <Real Name="Z">-16834.287828835641</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-11759.680211368686</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_14.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_14.xml
new file mode 100644 (file)
index 0000000..f9e3c86
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_14.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">6634.6501936621835</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">64134.915916992832</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">16834.287828835641</Real>
+          <Real Name="Z">-47300.628088157187</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-5074.6076174669543</Real>
+          <Real Name="Z">-16834.287828835641</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-11759.680211368686</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_15.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_15.xml
new file mode 100644 (file)
index 0000000..57c320b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_15.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">8420.0362853425831</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">48868.881042349596</Real>
+          <Real Name="Y">30612.527407106132</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-20367.569274272631</Real>
+          <Real Name="Y">-44370.122623426825</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-43273.123232299702</Real>
+          <Real Name="Y">31269.684749268814</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">14771.811464222737</Real>
+          <Real Name="Y">-17512.089532948117</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_16.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_16.xml
new file mode 100644 (file)
index 0000000..57c320b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_16.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">8420.0362853425831</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">48868.881042349596</Real>
+          <Real Name="Y">30612.527407106132</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-20367.569274272631</Real>
+          <Real Name="Y">-44370.122623426825</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-43273.123232299702</Real>
+          <Real Name="Y">31269.684749268814</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">14771.811464222737</Real>
+          <Real Name="Y">-17512.089532948117</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_17.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_17.xml
new file mode 100644 (file)
index 0000000..57c320b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_17.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">8420.0362853425831</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">48868.881042349596</Real>
+          <Real Name="Y">30612.527407106132</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-20367.569274272631</Real>
+          <Real Name="Y">-44370.122623426825</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-43273.123232299702</Real>
+          <Real Name="Y">31269.684749268814</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">14771.811464222737</Real>
+          <Real Name="Y">-17512.089532948117</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_18.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_18.xml
new file mode 100644 (file)
index 0000000..f3d9429
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_18.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.189967472820328</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-0.058802567113213383</Real>
+          <Real Name="Y">0.98008709084019618</Real>
+          <Real Name="Z">-0.98003392409622847</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.058270584850865621</Real>
+          <Real Name="Y">-0.92688886460541997</Real>
+          <Real Name="Z">0.87374380511461147</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">1.9600678481924569</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.00053198226234776186</Real>
+          <Real Name="Y">-0.053198226234776214</Real>
+          <Real Name="Z">-1.8537777292108399</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_19.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_19.xml
new file mode 100644 (file)
index 0000000..f3d9429
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_19.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.189967472820328</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-0.058802567113213383</Real>
+          <Real Name="Y">0.98008709084019618</Real>
+          <Real Name="Z">-0.98003392409622847</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.058270584850865621</Real>
+          <Real Name="Y">-0.92688886460541997</Real>
+          <Real Name="Z">0.87374380511461147</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">1.9600678481924569</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.00053198226234776186</Real>
+          <Real Name="Y">-0.053198226234776214</Real>
+          <Real Name="Z">-1.8537777292108399</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_2.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_2.xml
new file mode 100644 (file)
index 0000000..1a412cb
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_2.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">117.28019586833446</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">926.23638559389656</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">13.025199172414171</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-1186.74036904218</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-19.971972064368394</Real>
+          <Real Name="Y">115.77954819923707</Real>
+          <Real Name="Z">144.72443524904634</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">6.9467728919542244</Real>
+          <Real Name="Y">-115.77954819923707</Real>
+          <Real Name="Z">115.77954819923707</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_20.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_20.xml
new file mode 100644 (file)
index 0000000..f3d9429
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_20.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.189967472820328</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-0.058802567113213383</Real>
+          <Real Name="Y">0.98008709084019618</Real>
+          <Real Name="Z">-0.98003392409622847</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.058270584850865621</Real>
+          <Real Name="Y">-0.92688886460541997</Real>
+          <Real Name="Z">0.87374380511461147</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">1.9600678481924569</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0.00053198226234776186</Real>
+          <Real Name="Y">-0.053198226234776214</Real>
+          <Real Name="Z">-1.8537777292108399</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_21.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_21.xml
new file mode 100644 (file)
index 0000000..49465db
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_21.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.0062300725190453932</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.97676043227863252</Real>
+          <Real Name="Z">-0.048263050343893665</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-1.0051085931258372</Real>
+          <Real Name="Z">0.02857552486453363</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-0.039670992865678656</Real>
+          <Real Name="Z">1.3461222560797761</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.068019153712883362</Real>
+          <Real Name="Z">-1.3264347306004161</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_22.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_22.xml
new file mode 100644 (file)
index 0000000..49465db
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_22.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.0062300725190453932</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.97676043227863252</Real>
+          <Real Name="Z">-0.048263050343893665</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-1.0051085931258372</Real>
+          <Real Name="Z">0.02857552486453363</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-0.039670992865678656</Real>
+          <Real Name="Z">1.3461222560797761</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.068019153712883362</Real>
+          <Real Name="Z">-1.3264347306004161</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_23.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_23.xml
new file mode 100644 (file)
index 0000000..49465db
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_23.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.0062300725190453932</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.97676043227863252</Real>
+          <Real Name="Z">-0.048263050343893665</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-1.0051085931258372</Real>
+          <Real Name="Z">0.02857552486453363</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-0.039670992865678656</Real>
+          <Real Name="Z">1.3461222560797761</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0.068019153712883362</Real>
+          <Real Name="Z">-1.3264347306004161</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_24.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_24.xml
new file mode 100644 (file)
index 0000000..e6bf0d6
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_24.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.036552277896660712</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-1.2610253745363624</Real>
+          <Real Name="Y">1.5333660304021515</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.3353509548743419</Real>
+          <Real Name="Y">-1.5916299364299373</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.0035104610033714</Real>
+          <Real Name="Y">0.68555613834201967</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1.0778360413413508</Real>
+          <Real Name="Y">-0.62729223231423381</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_25.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_25.xml
new file mode 100644 (file)
index 0000000..e6bf0d6
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_25.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.036552277896660712</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-1.2610253745363624</Real>
+          <Real Name="Y">1.5333660304021515</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.3353509548743419</Real>
+          <Real Name="Y">-1.5916299364299373</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.0035104610033714</Real>
+          <Real Name="Y">0.68555613834201967</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1.0778360413413508</Real>
+          <Real Name="Y">-0.62729223231423381</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_26.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_26.xml
new file mode 100644 (file)
index 0000000..e6bf0d6
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_26.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="THOLE">
+    <FEP Name="No">
+      <Real Name="Epot ">-0.036552277896660712</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-1.2610253745363624</Real>
+          <Real Name="Y">1.5333660304021515</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.3353509548743419</Real>
+          <Real Name="Y">-1.5916299364299373</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1.0035104610033714</Real>
+          <Real Name="Y">0.68555613834201967</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1.0778360413413508</Real>
+          <Real Name="Y">-0.62729223231423381</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_27.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_27.xml
new file mode 100644 (file)
index 0000000..aa00c5a
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_27.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">976.98762692485559</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-195.3779875862125</Real>
+          <Real Name="Y">19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">195.3779875862125</Real>
+          <Real Name="Y">-19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_28.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_28.xml
new file mode 100644 (file)
index 0000000..aa00c5a
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_28.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">976.98762692485559</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-195.3779875862125</Real>
+          <Real Name="Y">19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">195.3779875862125</Real>
+          <Real Name="Y">-19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_29.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_29.xml
new file mode 100644 (file)
index 0000000..aa00c5a
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_29.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">976.98762692485559</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">-195.3779875862125</Real>
+          <Real Name="Y">19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">195.3779875862125</Real>
+          <Real Name="Y">-19537.798758621251</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_3.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_3.xml
new file mode 100644 (file)
index 0000000..02d659f
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_3.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">71.870154544676424</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">694.67728919542242</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">182.35278841379841</Real>
+          <Real Name="Z">-512.32450078162401</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-54.995285394637648</Real>
+          <Real Name="Z">-182.35278841379841</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-127.35750301916076</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_30.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_30.xml
new file mode 100644 (file)
index 0000000..4087b90
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_30.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">7893.2706984829874</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">35168.037765518253</Real>
+          <Real Name="Z">42983.157268966759</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-35168.037765518253</Real>
+          <Real Name="Z">-42983.157268966759</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_31.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_31.xml
new file mode 100644 (file)
index 0000000..4087b90
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_31.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">7893.2706984829874</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">35168.037765518253</Real>
+          <Real Name="Z">42983.157268966759</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-35168.037765518253</Real>
+          <Real Name="Z">-42983.157268966759</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_32.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_32.xml
new file mode 100644 (file)
index 0000000..4087b90
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_32.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">7893.2706984829874</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">35168.037765518253</Real>
+          <Real Name="Z">42983.157268966759</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-35168.037765518253</Real>
+          <Real Name="Z">-42983.157268966759</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_33.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_33.xml
new file mode 100644 (file)
index 0000000..0f1a4f1
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_33.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">3973.1129741191753</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">18521.833223172944</Real>
+          <Real Name="Y">34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-18521.833223172944</Real>
+          <Real Name="Y">-34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_34.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_34.xml
new file mode 100644 (file)
index 0000000..0f1a4f1
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_34.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">3973.1129741191753</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">18521.833223172944</Real>
+          <Real Name="Y">34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-18521.833223172944</Real>
+          <Real Name="Y">-34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_35.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_35.xml
new file mode 100644 (file)
index 0000000..0f1a4f1
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_35.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="WATERPOL">
+    <FEP Name="No">
+      <Real Name="Epot ">3973.1129741191753</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">18521.833223172944</Real>
+          <Real Name="Y">34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-18521.833223172944</Real>
+          <Real Name="Y">-34777.281790345827</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_4.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_4.xml
new file mode 100644 (file)
index 0000000..02d659f
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_4.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">71.870154544676424</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">694.67728919542242</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">182.35278841379841</Real>
+          <Real Name="Z">-512.32450078162401</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-54.995285394637648</Real>
+          <Real Name="Z">-182.35278841379841</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-127.35750301916076</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_5.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_5.xml
new file mode 100644 (file)
index 0000000..02d659f
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_5.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">71.870154544676424</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">694.67728919542242</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">182.35278841379841</Real>
+          <Real Name="Z">-512.32450078162401</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-54.995285394637648</Real>
+          <Real Name="Z">-182.35278841379841</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">-127.35750301916076</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_6.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_6.xml
new file mode 100644 (file)
index 0000000..22bc607
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_6.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">91.193681540689695</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">529.34409436691192</Real>
+          <Real Name="Y">331.59262604261494</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-220.64687398069611</Real>
+          <Real Name="Y">-480.60090457503304</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-468.4729969011629</Real>
+          <Real Name="Y">338.42361938636998</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">159.77577651494713</Real>
+          <Real Name="Y">-189.41534085395185</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_7.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_7.xml
new file mode 100644 (file)
index 0000000..22bc607
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_7.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">91.193681540689695</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">529.34409436691192</Real>
+          <Real Name="Y">331.59262604261494</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-220.64687398069611</Real>
+          <Real Name="Y">-480.60090457503304</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-468.4729969011629</Real>
+          <Real Name="Y">338.42361938636998</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">159.77577651494713</Real>
+          <Real Name="Y">-189.41534085395185</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_8.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_8.xml
new file mode 100644 (file)
index 0000000..22bc607
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_8.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="POLARIZATION">
+    <FEP Name="No">
+      <Real Name="Epot ">91.193681540689695</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">529.34409436691192</Real>
+          <Real Name="Y">331.59262604261494</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-220.64687398069611</Real>
+          <Real Name="Y">-480.60090457503304</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-468.4729969011629</Real>
+          <Real Name="Y">338.42361938636998</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">159.77577651494713</Real>
+          <Real Name="Y">-189.41534085395185</Real>
+          <Real Name="Z">0</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>

diff --git a/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_9.xml b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_9.xml
new file mode 100644 (file)
index 0000000..1a52d5b
--- /dev/null
+++ b/src/gromacs/listed_forces/tests/refdata/Polarize_ListedForcesTest_Ifunc_9.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <FunctionType Name="ANHARM_POL">
+    <FEP Name="No">
+      <Real Name="Epot ">10827.481942307291</Real>
+      <Real Name="dVdlambda ">0</Real>
+      <Vector Name="Central shift forces">
+        <Real Name="X">0</Real>
+        <Real Name="Y">0</Real>
+        <Real Name="Z">0</Real>
+      </Vector>
+      <Sequence Name="Forces">
+        <Int Name="Length">4</Int>
+        <Vector>
+          <Real Name="X">0</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">85527.567362390473</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">1202.4530379358848</Real>
+          <Real Name="Y">0</Real>
+          <Real Name="Z">-109576.62812110817</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">-1844.0698251282447</Real>
+          <Real Name="Y">10693.613119872665</Real>
+          <Real Name="Z">13355.447638845031</Real>
+        </Vector>
+        <Vector>
+          <Real Name="X">641.61678719235988</Real>
+          <Real Name="Y">-10693.613119872665</Real>
+          <Real Name="Z">10693.613119872665</Real>
+        </Vector>
+      </Sequence>
+    </FEP>
+  </FunctionType>
+</ReferenceData>