Changeset bcb593 for tests/regression/Molecules/VerletIntegration/post

Timestamp:

Jul 11, 2013, 9:43:31 PM (12 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Action_Thermostats, Add_AtomRandomPerturbation, Add_FitFragmentPartialChargesAction, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Added_ParseSaveFragmentResults, AddingActions_SaveParseParticleParameters, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_ParticleName_to_Atom, Adding_StructOpt_integration_tests, AtomFragments, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.5.4, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, CombiningParticlePotentialParsing, Combining_Subpackages, Debian_Package_split, Debian_package_split_molecuildergui_only, Disabling_MemDebug, Docu_Python_wait, EmpiricalPotential_contain_HomologyGraph, EmpiricalPotential_contain_HomologyGraph_documentation, Enable_parallel_make_install, Enhance_userguide, Enhanced_StructuralOptimization, Enhanced_StructuralOptimization_continued, Example_ManyWaysToTranslateAtom, Exclude_Hydrogens_annealWithBondGraph, FitPartialCharges_GlobalError, Fix_BoundInBox_CenterInBox_MoleculeActions, Fix_ChargeSampling_PBC, Fix_ChronosMutex, Fix_FitPartialCharges, Fix_FitPotential_needs_atomicnumbers, Fix_ForceAnnealing, Fix_IndependentFragmentGrids, Fix_ParseParticles, Fix_ParseParticles_split_forward_backward_Actions, Fix_PopActions, Fix_QtFragmentList_sorted_selection, Fix_Restrictedkeyset_FragmentMolecule, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, Fix_fitting_potentials, Fixes, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, FragmentAction_writes_AtomFragments, FragmentMolecule_checks_bonddegrees, GeometryObjects, Gui_Fixes, Gui_displays_atomic_force_velocity, ImplicitCharges, IndependentFragmentGrids, IndependentFragmentGrids_IndividualZeroInstances, IndependentFragmentGrids_IntegrationTest, IndependentFragmentGrids_Sole_NN_Calculation, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, MoreRobust_FragmentAutomation, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PdbParser_setsAtomName, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, Rewrite_FitPartialCharges, RotateToPrincipalAxisSystem_UndoRedo, SaturateAtoms_findBestMatching, SaturateAtoms_singleDegree, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, Switchable_LogView, ThirdParty_MPQC_rebuilt_buildsystem, TrajectoryDependenant_MaxOrder, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, TremoloParser_setsAtomName, Ubuntu_1604_changes, stable

Children:

9c1324

Parents:

11f0fa

git-author:

Frederik Heber <heber@…> (06/18/13 14:35:36)

git-committer:

Frederik Heber <heber@…> (07/11/13 21:43:31)

Message:

FIX: VerletIntegrationAction now assumes forces have just been calculated.

• according to Wikipedia's Velocity_Verlet info, we first integrate position, then calculate forces, then integrate velocity. This assumes that forces based on next time step's position are already known. This is possible when parsed from file but not if they are calculated dynamically via fragmentation. Hence, we now integrate velocity from last time step to current, then integrate position from current time step to next. Then we are in the position to calculate forces and do this cycle again.
• for this, VelocityVerletUpdate was split into ..X() and ..U() for independent integration of position and velocity.
• VelocityVerletIntegration::operator() now first corrects forces, then integrates velocites, corrects them too and finally integrates positions according to the above new scheme.
• removed option MDSteps from VerletIntegrationAction.
• TESTFIX: Had to change regression rest on VerletIntegration since the cycle sequence has changed. It was nonsense before to have the forces already in some file respective to future position that actually first need to come out of the time integration.

File:

• tests/regression/Molecules/VerletIntegration/post/test.conf (modified) (1 diff)

tests/regression/Molecules/VerletIntegration/post/test.conf

@@ -71 +71 @@
 Ion_Type2       3       6       1.0     3       3       12.01100000000  Carbon  C
 #Ion_TypeNr._Nr.R[0]    R[1]    R[2]    MoveType (0 MoveIon, 1 FixedIon)
-Ion_Type1_1     9.782085945     2.645886050     2.645886050     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 0
-Ion_Type1_2     9.782085945     2.645886050     4.425886024     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 1
-Ion_Type1_3     10.672039608    3.904536878     3.535886037     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 2
-Ion_Type1_4     8.532785963     4.787886018     2.645886050     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 3
-Ion_Type1_5     8.532785963     4.787886018     4.425886024     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 4
-Ion_Type1_6     6.393632318     3.904536877     3.535886037     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 5
-Ion_Type1_7     7.283585982     2.645886050     2.645886050     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 6
-Ion_Type1_8     7.283585982     2.645886050     4.425886024     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 7
-Ion_Type2_1     9.782085945     3.275186040     3.535886037     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 8
-Ion_Type2_2     8.532785963     4.158586027     3.535886037     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 9
-Ion_Type2_3     7.283585982     3.275186040     3.535886037     0       1.000000e-01    0.000000e+00    0.000000e+00     # molecule nr 10
+Ion_Type1_1     9.782085945     2.645886050     2.645886050     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 0
+Ion_Type1_2     9.782085945     2.645886050     4.425886024     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 1
+Ion_Type1_3     10.672039608    3.904536878     3.535886037     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 2
+Ion_Type1_4     8.532785963     4.787886018     2.645886050     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 3
+Ion_Type1_5     8.532785963     4.787886018     4.425886024     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 4
+Ion_Type1_6     6.393632318     3.904536877     3.535886037     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 5
+Ion_Type1_7     7.283585982     2.645886050     2.645886050     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 6
+Ion_Type1_8     7.283585982     2.645886050     4.425886024     0       1.003242e-01    5.239287e-04    0.000000e+00     # molecule nr 7
+Ion_Type2_1     9.782085945     3.275186040     3.535886037     0       9.984420e-02    4.396970e-05    0.000000e+00     # molecule nr 8
+Ion_Type2_2     8.532785963     4.158586027     3.535886037     0       9.984420e-02    4.396970e-05    0.000000e+00     # molecule nr 9
+Ion_Type2_3     7.283585982     3.275186040     3.535886037     0       9.984420e-02    4.396970e-05    0.000000e+00     # molecule nr 10
 
-Ion_Type1_1     9.783090905     2.645891010     2.645886050     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 0
-Ion_Type1_2     9.783090905     2.645891010     4.425886024     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 1
-Ion_Type1_3     10.673044568    3.904541838     3.535886037     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 2
-Ion_Type1_4     8.533790923     4.787890978     2.645886050     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 3
-Ion_Type1_5     8.533790923     4.787890978     4.425886024     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 4
-Ion_Type1_6     6.394637278     3.904541837     3.535886037     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 5
-Ion_Type1_7     7.284590942     2.645891010     2.645886050     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 6
-Ion_Type1_8     7.284590942     2.645891010     4.425886024     0       1.002953e-01    4.950408e-04    0.000000e+00     # molecule nr 7
-Ion_Type2_1     9.783086361     3.275186456     3.535886037     0       9.984178e-02    4.154535e-05    0.000000e+00     # molecule nr 8
-Ion_Type2_2     8.533786379     4.158586443     3.535886037     0       9.984178e-02    4.154535e-05    0.000000e+00     # molecule nr 9
-Ion_Type2_3     7.284586398     3.275186456     3.535886037     0       9.984178e-02    4.154535e-05    0.000000e+00     # molecule nr 10
+Ion_Type1_1     9.783091811     2.645893914     2.645886050     0 # molecule nr 0
+Ion_Type1_2     9.783091811     2.645893914     4.425886024     0 # molecule nr 1
+Ion_Type1_3     10.673045474    3.904544742     3.535886037     0 # molecule nr 2
+Ion_Type1_4     8.533791829     4.787893882     2.645886050     0 # molecule nr 3
+Ion_Type1_5     8.533791829     4.787893882     4.425886024     0 # molecule nr 4
+Ion_Type1_6     6.394638184     3.904544741     3.535886037     0 # molecule nr 5
+Ion_Type1_7     7.284591848     2.645893914     2.645886050     0 # molecule nr 6
+Ion_Type1_8     7.284591848     2.645893914     4.425886024     0 # molecule nr 7
+Ion_Type2_1     9.783084607     3.275186700     3.535886037     0 # molecule nr 8
+Ion_Type2_2     8.533784625     4.158586687     3.535886037     0 # molecule nr 9
+Ion_Type2_3     7.284584644     3.275186700     3.535886037     0 # molecule nr 10