Changeset 47d041 for src/Fragmentation/ForceMatrix.cpp
- 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, Candidate_v1.7.0, 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)
- File:
-
- 1 edited
-
src/Fragmentation/ForceMatrix.cpp (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
src/Fragmentation/ForceMatrix.cpp
r50e4e5 r47d041 43 43 stringstream line; 44 44 45 DoLog(0) && (Log() << Verbose(0) << "Parsing force indices for " << MatrixCounter << " matrices." << endl);45 LOG(0, "Parsing force indices for " << MatrixCounter << " matrices."); 46 46 Indices.resize(MatrixCounter + 1); 47 47 line << name << FRAGMENTPREFIX << FORCESFILE; 48 48 input.open(line.str().c_str(), ios::in); 49 //L og() << Verbose(0) << "Opening " << line.str() << " ... " << input << endl;49 //LOG(0, "Opening " << line.str() << " ... " << input); 50 50 if (input.fail()) { 51 DoLog(0) && (Log() << Verbose(0) << endl << "ForceMatrix::ParseIndices: Unable to open " << line.str() << ", is the directory correct?" << endl);51 LOG(0, endl << "ForceMatrix::ParseIndices: Unable to open " << line.str() << ", is the directory correct?"); 52 52 return false; 53 53 } … … 59 59 Indices[i].resize(RowCounter[i]); 60 60 FragmentNumber = FixedDigitNumber(MatrixCounter, i); 61 //Log() << Verbose(0) << FRAGMENTPREFIX << FragmentNumber << "[" << RowCounter[i] << "]:"; 61 //std::stringstream output; 62 //output << FRAGMENTPREFIX << FragmentNumber << "[" << RowCounter[i] << "]:"; 62 63 delete[](FragmentNumber); 63 64 for(int j=0;(j<RowCounter[i]) && (!line.eof());j++) { 64 65 line >> Indices[i][j]; 65 // Log() << Verbose(0)<< " " << Indices[i][j];66 //output << " " << Indices[i][j]; 66 67 } 67 //L og() << Verbose(0) << endl;68 //LOG(0, output.str()); 68 69 } 69 70 Indices[MatrixCounter].resize(RowCounter[MatrixCounter]); … … 92 93 int j = Indices[ FragmentNr ][l]; 93 94 if (j > RowCounter[MatrixCounter]) { 94 DoeLog(0) && (eLog()<< Verbose(0) << "Current force index " << j << " is greater than " << RowCounter[MatrixCounter] << "!" << endl);95 ELOG(0, "Current force index " << j << " is greater than " << RowCounter[MatrixCounter] << "!"); 95 96 performCriticalExit(); 96 97 return false; 97 98 } 98 99 if (j != -1) { 99 //if (j == 0) L og() << Verbose(0) << "Summing onto ion 0, type 0 from fragment " << FragmentNr << ", ion " << l << "." << endl;100 //if (j == 0) LOG(0, "Summing onto ion 0, type 0 from fragment " << FragmentNr << ", ion " << l << "."); 100 101 for(int k=2;k<ColumnCounter[MatrixCounter];k++) 101 102 Matrix[MatrixCounter][j][k] += sign*Fragments.Matrix[ FragmentNr ][l][k]; … … 128 129 input.open(file.str().c_str(), ios::in); 129 130 if (input.fail()) { 130 DoLog(0) && (Log() << Verbose(0) << endl << "ForceMatrix::ParseFragmentMatrix: Unable to open " << file.str() << ", is the directory correct?" << endl);131 LOG(0, endl << "ForceMatrix::ParseFragmentMatrix: Unable to open " << file.str() << ", is the directory correct?"); 131 132 return false; 132 133 } … … 137 138 while (!zeile.eof()) { 138 139 zeile >> nr; 139 //L og() << Verbose(0) << "Current index: " << getNr() << "." << endl;140 //LOG(0, "Current index: " << getNr() << "."); 140 141 if (nr > RowCounter[MatrixCounter]) 141 142 RowCounter[MatrixCounter] = nr;
Note:
See TracChangeset
for help on using the changeset viewer.
