Changeset 47d041 for src/Actions/MoleculeAction
- Timestamp:
- Nov 3, 2011, 7:44:01 PM (14 years ago)
- 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:
- 41a467
- Parents:
- 50e4e5
- git-author:
- Frederik Heber <heber@…> (10/27/11 11:53:58)
- git-committer:
- Frederik Heber <heber@…> (11/03/11 19:44:01)
- Location:
- src/Actions/MoleculeAction
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/MoleculeAction/FillVoidWithMoleculeAction.cpp
r50e4e5 r47d041 51 51 52 52 if (!boost::filesystem::exists(params.fillername)) { 53 DoeLog(1) && (eLog() << Verbose(1) << "File with filler molecule " << params.fillername << " does not exist!" << endl);53 ELOG(1, "File with filler molecule " << params.fillername << " does not exist!"); 54 54 return Action::failure; 55 55 } 56 56 57 DoLog(1) && (Log() << Verbose(1) << "Filling Box with water molecules, "57 LOG(1, "INFO: Filling Box with water molecules, " 58 58 << " minimum distance to molecules" << params.boundary 59 59 << ", random atom displacement " << params.RandAtomDisplacement … … 61 61 << ", distances between fillers (" << params.distances[0] << "," << params.distances[1] << "," << params.distances[2] 62 62 << "), MinDistance " << params.MinDistance 63 << ", DoRotate " << params.DoRotate << "." << endl);63 << ", DoRotate " << params.DoRotate << "."); 64 64 // construct water molecule 65 65 std::vector<molecule *> presentmolecules = World::getInstance().getAllMolecules(); 66 // DoLog(0) && (Log() << Verbose(0) << presentmolecules.size() << " molecules initially are present." << std::endl);66 // LOG(0, presentmolecules.size() << " molecules initially are present."); 67 67 std::string FilenameSuffix = params.fillername.string().substr(params.fillername.string().find_last_of('.')+1, params.fillername.string().length()); 68 68 ifstream input; … … 95 95 // (we can in general remove more quickly from a list than a vector) 96 96 std::vector<molecule *> fillermolecules = World::getInstance().getAllMolecules(); 97 // DoLog(0) && (Log() << Verbose(0) << fillermolecules.size() << " molecules are present." << std::endl);97 // LOG(0, fillermolecules.size() << " molecules are present."); 98 98 std::list<molecule *> fillermolecules_list; 99 99 std::copy( fillermolecules.begin(), fillermolecules.end(), std::back_inserter( fillermolecules_list )); 100 // DoLog(0) && (Log() << Verbose(0) << fillermolecules_list.size() << " molecules have been copied." << std::endl);100 // LOG(0, fillermolecules_list.size() << " molecules have been copied."); 101 101 for (std::vector<molecule *>::const_iterator iter = presentmolecules.begin(); 102 102 iter != presentmolecules.end(); … … 104 104 fillermolecules_list.remove(*iter); 105 105 } 106 // DoLog(0) && (Log() << Verbose(0) << fillermolecules_list.size() << " molecules left after removal." << std::endl);106 // LOG(0, fillermolecules_list.size() << " molecules left after removal."); 107 107 fillermolecules.clear(); 108 108 std::copy(fillermolecules_list.begin(), fillermolecules_list.end(), std::back_inserter( fillermolecules )); 109 109 110 // DoLog(0) && (Log() << Verbose(0) << fillermolecules.size() << " molecules have been inserted." << std::endl);110 // LOG(0, fillermolecules.size() << " molecules have been inserted."); 111 111 112 112 return Action::state_ptr(new MoleculeFillVoidWithMoleculeState(fillermolecules,params)); -
src/Actions/MoleculeAction/FillWithMoleculeAction.cpp
r50e4e5 r47d041 49 49 getParametersfromValueStorage(); 50 50 51 DoLog(1) && (Log() << Verbose(1) << "Filling Box with water molecules, "51 LOG(1, "INFO: Filling Box with water molecules, " 52 52 << " minimum distance to molecules" << params.boundary 53 53 << ", random atom displacement " << params.RandAtomDisplacement … … 55 55 << ", distances between fillers (" << params.distances[0] << "," << params.distances[1] << "," << params.distances[2] 56 56 << "), MinDistance " << params.MaxDistance 57 << ", DoRotate " << params.DoRotate << "." << endl);57 << ", DoRotate " << params.DoRotate << "."); 58 58 // construct water molecule 59 59 std::vector<molecule *> presentmolecules = World::getInstance().getAllMolecules(); 60 // DoLog(0) && (Log() << Verbose(0) << presentmolecules.size() << " molecules initially are present." << std::endl);60 // LOG(0, presentmolecules.size() << " molecules initially are present."); 61 61 std::string FilenameSuffix = params.fillername.string().substr(params.fillername.string().find_last_of('.')+1, params.fillername.string().length()); 62 62 ifstream input; … … 113 113 // (we can in general remove more quickly from a list than a vector) 114 114 std::vector<molecule *> fillermolecules = World::getInstance().getAllMolecules(); 115 // DoLog(0) && (Log() << Verbose(0) << fillermolecules.size() << " molecules are present." << std::endl);115 // LOG(0, fillermolecules.size() << " molecules are present."); 116 116 std::list<molecule *> fillermolecules_list; 117 117 std::copy( fillermolecules.begin(), fillermolecules.end(), std::back_inserter( fillermolecules_list )); 118 // DoLog(0) && (Log() << Verbose(0) << fillermolecules_list.size() << " molecules have been copied." << std::endl);118 // LOG(0, fillermolecules_list.size() << " molecules have been copied."); 119 119 for (std::vector<molecule *>::const_iterator iter = presentmolecules.begin(); 120 120 iter != presentmolecules.end(); … … 122 122 fillermolecules_list.remove(*iter); 123 123 } 124 // DoLog(0) && (Log() << Verbose(0) << fillermolecules_list.size() << " molecules left after removal." << std::endl);124 // LOG(0, fillermolecules_list.size() << " molecules left after removal."); 125 125 fillermolecules.clear(); 126 126 std::copy(fillermolecules_list.begin(), fillermolecules_list.end(), std::back_inserter( fillermolecules )); 127 127 128 // DoLog(0) && (Log() << Verbose(0) << fillermolecules.size() << " molecules have been inserted." << std::endl);128 // LOG(0, fillermolecules.size() << " molecules have been inserted."); 129 129 130 130 return Action::state_ptr(new MoleculeFillWithMoleculeState(fillermolecules,params)); -
src/Actions/MoleculeAction/LoadAction.cpp
r50e4e5 r47d041 48 48 // parsing file if present 49 49 if (!boost::filesystem::exists(params.filename)) { 50 DoLog(1) && (Log() << Verbose(1) << "Specified input file " << params.filename << " not found." << endl);50 LOG(1, "Specified input file " << params.filename << " not found."); 51 51 return Action::failure; 52 52 } else { 53 DoLog(1) && (Log() << Verbose(1) <<"Specified input file found, parsing ... ");53 LOG(1, "Specified input file found, parsing ... "); 54 54 55 55 // extract suffix … … 66 66 #endif 67 67 } else { 68 DoeLog(1) && (eLog() << Verbose(1) << "Input file does not have a suffix, cannot recognize format." << endl);68 ELOG(1, "Input file does not have a suffix, cannot recognize format."); 69 69 return Action::failure; 70 70 } … … 89 89 MoleculeList::const_reverse_iterator iter = World::getInstance().getMolecules()->ListOfMolecules.rbegin(); 90 90 (*iter)->SetNameFromFilename(FilenamePrefix.c_str()); 91 DoLog(0) && (Log() << Verbose(0) << "Chemical formula is " << (*iter)->getFormula() << std::endl);91 LOG(0, "Chemical formula is " << (*iter)->getFormula()); 92 92 93 93 return Action::state_ptr( … … 136 136 (*iter)->SetNameFromFilename(state->FilenamePrefix.c_str()); 137 137 (*iter)->setId(state->molId); 138 DoLog(0) && (Log() << Verbose(0) << "Chemical formula is " << (*iter)->getFormula() << std::endl);138 LOG(0, "Chemical formula is " << (*iter)->getFormula()); 139 139 140 140 return Action::state_ptr( -
src/Actions/MoleculeAction/RotateAroundSelfByAngleAction.cpp
r50e4e5 r47d041 59 59 // Creation Line that is the rotation axis 60 60 Vector *CenterOfGravity = mol->DetermineCenterOfGravity(); 61 DoLog(0) && (Log() << Verbose(0) << "Center of gravity is " << *CenterOfGravity << "." << std::endl);61 LOG(0, "Center of gravity is " << *CenterOfGravity << "."); 62 62 Line RotationAxis(*CenterOfGravity, params.Axis); 63 63 delete(CenterOfGravity); 64 DoLog(0) && (Log() << Verbose(0) << "Rotate " << mol->getName() << " around self by " << params.angle << " radian around axis " << RotationAxis << "." << endl);64 LOG(0, "Rotate " << mol->getName() << " around self by " << params.angle << " radian around axis " << RotationAxis << "."); 65 65 66 66 for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) { 67 67 (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), params.angle)); 68 68 } 69 DoLog(0) && (Log() << Verbose(0) << "done." << endl);69 LOG(0, "done."); 70 70 } 71 71 … … 78 78 BOOST_FOREACH(molecule *mol, state->selectedMolecules) { 79 79 Vector *CenterOfGravity = mol->DetermineCenterOfGravity(); 80 DoLog(0) && (Log() << Verbose(0) << "Center of gravity is " << *CenterOfGravity << "." << std::endl);80 LOG(0, "Center of gravity is " << *CenterOfGravity << "."); 81 81 Line RotationAxis(*CenterOfGravity, state->params.Axis); 82 82 delete(CenterOfGravity); 83 DoLog(0) && (Log() << Verbose(0) << "Rotate " << mol->getName() << " around self by " << -state->params.angle << " radian around axis " << RotationAxis << "." << endl);83 LOG(0, "Rotate " << mol->getName() << " around self by " << -state->params.angle << " radian around axis " << RotationAxis << "."); 84 84 85 85 for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) { … … 96 96 BOOST_FOREACH(molecule *mol, state->selectedMolecules) { 97 97 Vector *CenterOfGravity = mol->DetermineCenterOfGravity(); 98 DoLog(0) && (Log() << Verbose(0) << "Center of gravity is " << *CenterOfGravity << "." << std::endl);98 LOG(0, "Center of gravity is " << *CenterOfGravity << "."); 99 99 Line RotationAxis(*CenterOfGravity, state->params.Axis); 100 100 delete(CenterOfGravity); 101 DoLog(0) && (Log() << Verbose(0) << "Rotate " << mol->getName() << " around self by " << state->params.angle << " radian around axis " << RotationAxis << "." << endl);101 LOG(0, "Rotate " << mol->getName() << " around self by " << state->params.angle << " radian around axis " << RotationAxis << "."); 102 102 103 103 for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) { -
src/Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.cpp
r50e4e5 r47d041 48 48 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { 49 49 mol = iter->second; 50 DoLog(0) && (Log() << Verbose(0) << "Converting to prinicipal axis system." << endl);50 LOG(0, "Converting to prinicipal axis system."); 51 51 52 52 RealSpaceMatrix InertiaTensor = mol->getInertiaTensor(); -
src/Actions/MoleculeAction/SaveAdjacencyAction.cpp
r50e4e5 r47d041 47 47 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { 48 48 mol = iter->second; 49 DoLog(0) && (Log() << Verbose(0) << "Storing adjacency to path " << params.adjacencyfile << "." << endl);49 LOG(0, "Storing adjacency to path " << params.adjacencyfile << "."); 50 50 // TODO: sollte stream nicht filename benutzen, besser fuer unit test 51 51 #if BOOST_VERSION >= 104600 -
src/Actions/MoleculeAction/SaveBondsAction.cpp
r50e4e5 r47d041 47 47 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { 48 48 mol = iter->second; 49 DoLog(0) && (Log() << Verbose(0) << "Storing bonds to path " << params.bondsfile << "." << endl);49 LOG(0, "Storing bonds to path " << params.bondsfile << "."); 50 50 // TODO: sollte stream, nicht filenamen direkt nutzen, besser fuer unit tests 51 51 #if BOOST_VERSION >= 104600 -
src/Actions/MoleculeAction/SaveSelectedMoleculesAction.cpp
r50e4e5 r47d041 41 41 getParametersfromValueStorage(); 42 42 43 DoLog(1) && (Log() << Verbose(1) << "Storing world to file " << params.filename << "." << endl);43 LOG(1, "Storing world to file " << params.filename << "."); 44 44 45 45 // extract suffix … … 56 56 #endif 57 57 } else { 58 DoeLog(1) && (eLog() << Verbose(1) << "Output file does not have a suffix, cannot recognize format." << endl);58 ELOG(1, "Output file does not have a suffix, cannot recognize format."); 59 59 return Action::failure; 60 60 } … … 67 67 FormatParserStorage::getInstance().saveSelectedMolecules(output, FilenameSuffix); 68 68 } else { 69 DoeLog(1) && (eLog() << Verbose(1) << "Could not open file " << params.filename << "." << endl);69 ELOG(1, "Could not open file " << params.filename << "."); 70 70 } 71 71 output.close(); -
src/Actions/MoleculeAction/SaveTemperatureAction.cpp
r50e4e5 r47d041 45 45 getParametersfromValueStorage(); 46 46 47 DoLog(1) && (Log() << Verbose(1) << "Storing temperatures in " << params.temperaturefile << "." << endl);47 LOG(1, "Storing temperatures in " << params.temperaturefile << "."); 48 48 ofstream output; 49 49 output.open(params.temperaturefile.string().c_str(), ios::trunc); … … 52 52 OutputTemperature<std::vector<atom *> > writer(set); 53 53 if (output.fail() || !writer((ofstream * const) &output, 0, MDSteps)) 54 DoLog(2) && (Log() << Verbose(2) << "File could not be written." << endl);54 LOG(2, "File could not be written."); 55 55 else 56 DoLog(2) && (Log() << Verbose(2) << "File stored." << endl);56 LOG(2, "File stored."); 57 57 output.close(); 58 58 return Action::success; -
src/Actions/MoleculeAction/SuspendInWaterAction.cpp
r50e4e5 r47d041 46 46 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { 47 47 mol = iter->second; 48 DoLog(0) && (Log() << Verbose(0) <<"Evaluating necessary cell volume for a cluster suspended in water.");48 LOG(0, "Evaluating necessary cell volume for a cluster suspended in water."); 49 49 if (params.density < 1.0) { 50 DoeLog(1) && (eLog()<< Verbose(1) << "Density must be greater than 1.0g/cm^3!" << endl);50 ELOG(1, "Density must be greater than 1.0g/cm^3!"); 51 51 } else { 52 52 PrepareClustersinWater(World::getInstance().getConfig(), mol, volume, params.density); // if volume == 0, will calculate from ConvexEnvelope -
src/Actions/MoleculeAction/VerletIntegrationAction.cpp
r50e4e5 r47d041 46 46 47 47 48 DoLog(1) && (Log() << Verbose(1) << "Parsing forces file and Verlet integrating." << endl);48 LOG(1, "Parsing forces file and Verlet integrating."); 49 49 // TODO: sollte besser stream nutzen, nicht filename direkt (es sei denn, ist prefix), besser fuer unit test 50 50 char outputname[MAXSTRINGSIZE]; … … 54 54 VerletForceIntegration<std::vector<atom *> > Verlet(set, params.Deltat, step, false); 55 55 if (!Verlet(outputname, 1, 0, params.FixedCenterOfMass)) 56 DoLog(2) && (Log() << Verbose(2) << "File " << params.forcesfile << " not found." << endl);56 LOG(2, "File " << params.forcesfile << " not found."); 57 57 else 58 DoLog(2) && (Log() << Verbose(2) << "File " << params.forcesfile << " found and parsed." << endl);58 LOG(2, "File " << params.forcesfile << " found and parsed."); 59 59 } 60 60
Note:
See TracChangeset
for help on using the changeset viewer.