- Timestamp:
- Nov 23, 2009, 6:22:33 PM (15 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:
- e359a8
- Parents:
- f1ef60a
- Location:
- src
- Files:
-
- 22 edited
Legend:
- Unmodified
- Added
- Removed
-
src/analyzer.cpp
rf1ef60a r717e0c 94 94 if (!Hcorrection.ParseFragmentMatrix(argv[1], "", HCORRECTIONSUFFIX,0,0)) { 95 95 NoHCorrection = true; 96 Log() << Verbose(0) << "No HCorrection file found, skipping these." << endl;96 eLog() << Verbose(2) << "No HCorrection file found, skipping these." << endl; 97 97 } 98 98 … … 100 100 if (!Hessian.ParseFragmentMatrix(argv[1], dir, HessianSuffix,0,0)) { 101 101 NoHessian = true; 102 Log() << Verbose(0) << "No Hessian file found, skipping these." << endl;102 eLog() << Verbose(2) << "No Hessian file found, skipping these." << endl; 103 103 } 104 104 if (!Time.ParseFragmentMatrix(argv[1], dir, TimeSuffix, 10,1)) { 105 105 NoTime = true; 106 Log() << Verbose(0) << "No speed file found, skipping these." << endl;106 eLog() << Verbose(2) << "No speed file found, skipping these." << endl; 107 107 } 108 108 if (periode != NULL) { // also look for PAS values -
src/atom_bondedparticle.cpp
rf1ef60a r717e0c 44 44 void BondedParticle::OutputBondOfAtom() const 45 45 { 46 Log() << Verbose(4) << "Atom " << Name << "/" << nr << " with " << ListOfBonds.size() << " bonds: " ;46 Log() << Verbose(4) << "Atom " << Name << "/" << nr << " with " << ListOfBonds.size() << " bonds: " << endl; 47 47 int TotalDegree = 0; 48 48 for (BondList::const_iterator Runner = ListOfBonds.begin(); Runner != ListOfBonds.end(); ++Runner) { 49 Log() << Verbose( 0) << **Runner << "\t";49 Log() << Verbose(4) << **Runner << endl; 50 50 TotalDegree += (*Runner)->BondDegree; 51 51 } 52 Log() << Verbose( 0) << " -- TotalDegree: " << TotalDegree << endl;52 Log() << Verbose(4) << " -- TotalDegree: " << TotalDegree << endl; 53 53 }; 54 54 … … 75 75 status = true; 76 76 } else { 77 Log() << Verbose(1) << "ERROR: "<< *Binder << " does not contain " << *this << "." << endl;77 eLog() << Verbose(1) << *Binder << " does not contain " << *this << "." << endl; 78 78 } 79 79 } else { 80 Log() << Verbose(1) << "ERROR:Binder is " << Binder << "." << endl;80 eLog() << Verbose(1) << "Binder is " << Binder << "." << endl; 81 81 } 82 82 return status; … … 94 94 status = true; 95 95 } else { 96 Log() << Verbose(1) << "ERROR: "<< *Binder << " does not contain " << *this << "." << endl;96 eLog() << Verbose(1) << *Binder << " does not contain " << *this << "." << endl; 97 97 } 98 98 } else { 99 Log() << Verbose(1) << "ERROR:Binder is " << Binder << "." << endl;99 eLog() << Verbose(1) << "Binder is " << Binder << "." << endl; 100 100 } 101 101 return status; -
src/atom_graphnode.cpp
rf1ef60a r717e0c 29 29 Log() << Verbose(2) << "Atom " << Name << " is " << ((SeparationVertex) ? "a" : "not a") << " separation vertex, components are "; 30 30 OutputComponentNumber(); 31 Log() << Verbose( 0) << " with Lowpoint " << LowpointNr << " and Graph Nr. " << GraphNr << "." << endl;31 Log() << Verbose(3) << " with Lowpoint " << LowpointNr << " and Graph Nr. " << GraphNr << "." << endl; 32 32 }; 33 33 … … 40 40 if (ComponentNr != NULL) { 41 41 for (int i=0; ComponentNr[i] != -1; i++) 42 Log() << Verbose( 0) << ComponentNr[i] << " ";42 Log() << Verbose(2) << ComponentNr[i] << " "; 43 43 } 44 44 }; -
src/atom_trajectoryparticle.cpp
rf1ef60a r717e0c 203 203 U[d] = gsl_ran_gaussian (r, sigma); 204 204 } 205 Log() << Verbose( 0) << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl;205 Log() << Verbose(2) << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl; 206 206 } 207 207 for (int d=0; d<NDIM; d++) -
src/bond.cpp
rf1ef60a r717e0c 63 63 if(rightatom == Atom) 64 64 return leftatom; 65 eLog() << Verbose( 0) << "Bond " << *this << " does not contain atom " << *Atom << "!" << endl;65 eLog() << Verbose(1) << "Bond " << *this << " does not contain atom " << *Atom << "!" << endl; 66 66 return NULL; 67 67 }; … … 99 99 bool bond::MarkUsed(const enum Shading color) { 100 100 if (Used == black) { 101 eLog() << Verbose( 0) << "ERROR:Bond " << this << " was already marked black!." << endl;101 eLog() << Verbose(1) << "Bond " << this << " was already marked black!." << endl; 102 102 return false; 103 103 } else { -
src/bondgraph.cpp
rf1ef60a r717e0c 153 153 { 154 154 if (BondLengthMatrix == NULL) {// safety measure if no matrix has been parsed yet 155 eLog() << Verbose( 1) << "WARNING:BondLengthMatrixMinMaxDistance() called without having parsed the bond length matrix yet!" << endl;155 eLog() << Verbose(2) << "BondLengthMatrixMinMaxDistance() called without having parsed the bond length matrix yet!" << endl; 156 156 CovalentMinMaxDistance(Walker, OtherWalker, MinDistance, MaxDistance, IsAngstroem); 157 157 } else { -
src/boundary.cpp
rf1ef60a r717e0c 176 176 angle = 0.; // otherwise it's a vector in Axis Direction and unimportant for boundary issues 177 177 178 //Log() << Verbose( 0) << "Checking sign in quadrant : " << ProjectedVector.Projection(&AngleReferenceNormalVector) << "." << endl;178 //Log() << Verbose(2) << "Checking sign in quadrant : " << ProjectedVector.Projection(&AngleReferenceNormalVector) << "." << endl; 179 179 if (ProjectedVector.ScalarProduct(&AngleReferenceNormalVector) > 0) { 180 180 angle = 2. * M_PI - angle; … … 214 214 // for(Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++) { 215 215 // if (runner != BoundaryPoints[axis].begin()) 216 // Log() << Verbose( 0) << ", " << i << ": " << *runner->second.second;216 // Log() << Verbose(2) << ", " << i << ": " << *runner->second.second; 217 217 // else 218 // Log() << Verbose( 0) << i << ": " << *runner->second.second;218 // Log() << Verbose(2) << i << ": " << *runner->second.second; 219 219 // i++; 220 220 // } 221 // Log() << Verbose( 0) << endl;221 // Log() << Verbose(2) << endl; 222 222 // } 223 223 // 3c. throw out points whose distance is less than the mean of left and right neighbours … … 249 249 SideA.SubtractVector(MolCenter); 250 250 SideA.ProjectOntoPlane(&AxisVector); 251 // Log() << Verbose(0) << "SideA: "; 252 // SideA.Output(out); 253 // Log() << Verbose(0) << endl; 251 // Log() << Verbose(0) << "SideA: " << SideA << endl; 254 252 255 253 SideB.CopyVector(&right->second.second->x); 256 254 SideB.SubtractVector(MolCenter); 257 255 SideB.ProjectOntoPlane(&AxisVector); 258 // Log() << Verbose(0) << "SideB: "; 259 // SideB.Output(out); 260 // Log() << Verbose(0) << endl; 256 // Log() << Verbose(0) << "SideB: " << SideB << endl; 261 257 262 258 SideC.CopyVector(&left->second.second->x); 263 259 SideC.SubtractVector(&right->second.second->x); 264 260 SideC.ProjectOntoPlane(&AxisVector); 265 // Log() << Verbose(0) << "SideC: "; 266 // SideC.Output(out); 267 // Log() << Verbose(0) << endl; 261 // Log() << Verbose(0) << "SideC: " << SideC << endl; 268 262 269 263 SideH.CopyVector(&runner->second.second->x); 270 264 SideH.SubtractVector(MolCenter); 271 265 SideH.ProjectOntoPlane(&AxisVector); 272 // Log() << Verbose(0) << "SideH: "; 273 // SideH.Output(out); 274 // Log() << Verbose(0) << endl; 266 // Log() << Verbose(0) << "SideH: " << SideH << endl; 275 267 276 268 // calculate each length … … 335 327 for(Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++) { 336 328 if (runner != BoundaryPoints[axis].begin()) 337 Log() << Verbose( 0) << ", " << i << ": " << *runner->second.second;329 Log() << Verbose(2) << ", " << i << ": " << *runner->second.second; 338 330 else 339 Log() << Verbose( 0) << i << ": " << *runner->second.second;331 Log() << Verbose(2) << i << ": " << *runner->second.second; 340 332 i++; 341 333 } 342 Log() << Verbose( 0) << endl;334 Log() << Verbose(2) << endl; 343 335 } 344 336 … … 347 339 for (Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++) 348 340 if (!TesselStruct->AddBoundaryPoint(runner->second.second, 0)) 349 Log() << Verbose(3) << "WARNING:Point " << *(runner->second.second) << " is already present!" << endl;341 eLog() << Verbose(2) << "Point " << *(runner->second.second) << " is already present!" << endl; 350 342 351 343 Log() << Verbose(2) << "I found " << TesselStruct->PointsOnBoundaryCount << " points on the convex boundary." << endl; … … 355 347 // Log() << Verbose(1) << "Listing PointsOnBoundary:"; 356 348 // for(PointMap::iterator runner = PointsOnBoundary.begin(); runner != PointsOnBoundary.end(); runner++) { 357 // Log() << Verbose( 0) << " " << *runner->second;349 // Log() << Verbose(1) << " " << *runner->second; 358 350 // } 359 // Log() << Verbose( 0) << endl;351 // Log() << Verbose(1) << endl; 360 352 361 353 // 3a. guess starting triangle … … 406 398 // flip the line 407 399 if (TesselStruct->PickFarthestofTwoBaselines(line) == 0.) 408 Log() << Verbose(1) << "ERROR:Correction of concave baselines failed!" << endl;400 eLog() << Verbose(1) << "Correction of concave baselines failed!" << endl; 409 401 else { 410 402 TesselStruct->FlipBaseline(line); … … 462 454 463 455 if ((TesselStruct == NULL) || (TesselStruct->PointsOnBoundary.empty())) { 464 Log() << Verbose(2) << "ERROR:TesselStruct is empty." << endl;456 eLog() << Verbose(1) << "TesselStruct is empty." << endl; 465 457 return false; 466 458 } … … 470 462 Log() << Verbose(2) << "Remaining points are: "; 471 463 for (PointMap::iterator PointSprinter = TesselStruct->PointsOnBoundary.begin(); PointSprinter != TesselStruct->PointsOnBoundary.end(); PointSprinter++) 472 Log() << Verbose( 0) << *(PointSprinter->second) << "\t";473 Log() << Verbose( 0) << endl;464 Log() << Verbose(2) << *(PointSprinter->second) << "\t"; 465 Log() << Verbose(2) << endl; 474 466 475 467 PointRunner = TesselStruct->PointsOnBoundary.begin(); … … 528 520 // check whether there is something to work on 529 521 if (TesselStruct == NULL) { 530 Log() << Verbose(1) << "ERROR:TesselStruct is empty!" << endl;522 eLog() << Verbose(1) << "TesselStruct is empty!" << endl; 531 523 return volume; 532 524 } … … 756 748 Log() << Verbose(1) << "Minimum volume of the convex envelope contained in a rectangular box is " << minimumvolume << " atomicmassunit/" << (IsAngstroem ? "angstrom" : "atomiclength") << "^3." << endl; 757 749 if (minimumvolume > cellvolume) { 758 eLog() << Verbose( 0) << "ERROR:the containing box already has a greater volume than the envisaged cell volume!" << endl;750 eLog() << Verbose(1) << "the containing box already has a greater volume than the envisaged cell volume!" << endl; 759 751 Log() << Verbose(0) << "Setting Box dimensions to minimum possible, the greatest diameters." << endl; 760 752 for (int i = 0; i < NDIM; i++) … … 849 841 for(int i=0;i<NDIM;i++) { 850 842 N[i] = (int) ceil(1./FillerDistance.x[i]); 851 Log() << Verbose( 0) << N[i];843 Log() << Verbose(1) << N[i]; 852 844 if (i != NDIM-1) 853 Log() << Verbose( 0)<< ", ";845 Log() << Verbose(1)<< ", "; 854 846 else 855 Log() << Verbose( 0) << "." << endl;847 Log() << Verbose(1) << "." << endl; 856 848 } 857 849 … … 870 862 // get linked cell list 871 863 if (TesselStruct[i] == NULL) { 872 Log() << Verbose(1) << "ERROR:TesselStruct of " << (*ListRunner) << " is NULL. Didn't we pre-create it?" << endl;864 eLog() << Verbose(1) << "TesselStruct of " << (*ListRunner) << " is NULL. Didn't we pre-create it?" << endl; 873 865 FillIt = false; 874 866 } else { … … 1015 1007 //Log() << Verbose(1) << "Line " << *baseline->second << " has " << baseline->second->triangles.size() << " triangles adjacent" << endl; 1016 1008 if (baseline->second->triangles.size() != 2) 1017 eLog() << Verbose( 1) << "TESSELATION FINISHED WITH INVALID TRIANGLE COUNT!" << endl;1009 eLog() << Verbose(0) << "TESSELATION FINISHED WITH INVALID TRIANGLE COUNT!" << endl; 1018 1010 } 1019 1011 -
src/builder.cpp
rf1ef60a r717e0c 95 95 switch (choice) { 96 96 default: 97 Log() << Verbose(0) << "Not a valid choice." << endl;97 eLog() << Verbose(2) << "Not a valid choice." << endl; 98 98 break; 99 99 case 'a': // absolute coordinates of atom … … 109 109 valid = true; 110 110 do { 111 if (!valid) Log() << Verbose(0) << "Resulting position out of cell." << endl;111 if (!valid) eLog() << Verbose(2) << "Resulting position out of cell." << endl; 112 112 Log() << Verbose(0) << "Enter reference coordinates." << endl; 113 113 x.AskPosition(mol->cell_size, true); … … 125 125 valid = true; 126 126 do { 127 if (!valid) Log() << Verbose(0) << "Resulting position out of cell." << endl;127 if (!valid) eLog() << Verbose(2) << "Resulting position out of cell." << endl; 128 128 second = mol->AskAtom("Enter atom number: "); 129 129 Log() << Verbose(0) << "Enter relative coordinates." << endl; … … 142 142 do { 143 143 if (!valid) { 144 Log() << Verbose(0) << "Resulting coordinates out of cell - "; 145 first->x.Output(); 146 Log() << Verbose(0) << endl; 144 eLog() << Verbose(2) << "Resulting coordinates out of cell - " << first->x << endl; 147 145 } 148 146 Log() << Verbose(0) << "First, we need two atoms, the first atom is the central, while the second is the outer one." << endl; … … 671 669 Log() << Verbose(0) << "===============================================" << endl; 672 670 if (molecules->NumberOfActiveMolecules() > 1) 673 Log() << Verbose(0) << "WARNING:There is more than one molecule active! Atoms will be added to each." << endl;671 eLog() << Verbose(2) << "There is more than one molecule active! Atoms will be added to each." << endl; 674 672 Log() << Verbose(0) << "INPUT: "; 675 673 cin >> choice; … … 794 792 Log() << Verbose(0) << "===============================================" << endl; 795 793 if (molecules->NumberOfActiveMolecules() > 1) 796 Log() << Verbose(0) << "WARNING:There is more than one molecule active! Atoms will be added to each." << endl;794 eLog() << Verbose(2) << "There is more than one molecule active! Atoms will be added to each." << endl; 797 795 Log() << Verbose(0) << "INPUT: "; 798 796 cin >> choice; … … 827 825 } 828 826 if (count != j) 829 Log() << Verbose(0) << "ERROR:AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl;827 eLog() << Verbose(1) << "AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl; 830 828 x.Zero(); 831 829 y.Zero(); … … 1250 1248 1251 1249 if (!strcmp(configuration->configpath, configuration->GetDefaultPath())) { 1252 eLog() << Verbose( 0) << "WARNING:config is found under different path then stated in config file::defaultpath!" << endl;1250 eLog() << Verbose(2) << "config is found under different path then stated in config file::defaultpath!" << endl; 1253 1251 } 1254 1252 … … 1354 1352 1355 1353 if (!strcmp(configuration->configpath, configuration->GetDefaultPath())) { 1356 eLog() << Verbose( 0) << "WARNING:config is found under different path then stated in config file::defaultpath!" << endl;1354 eLog() << Verbose(2) << "config is found under different path then stated in config file::defaultpath!" << endl; 1357 1355 } 1358 1356 … … 1387 1385 molecule *mol = NULL; 1388 1386 string BondGraphFileName(""); 1387 int verbosity = 0; 1389 1388 strncpy(configuration.databasepath, LocalPath, MAXSTRINGSIZE-1); 1390 1389 … … 1432 1431 Log() << Verbose(0) << "\t-T x1 x2 x3\tTranslate periodically all atoms by this vector (x1,x2,x3)." << endl; 1433 1432 Log() << Verbose(0) << "\t-u rho\tsuspend in water solution and output necessary cell lengths, average density rho and repetition." << endl; 1434 Log() << Verbose(0) << "\t-v/-V\t\tGives version information." << endl; 1433 Log() << Verbose(0) << "\t-v\t\tsets verbosity (more is more)." << endl; 1434 Log() << Verbose(0) << "\t-V\t\tGives version information." << endl; 1435 1435 Log() << Verbose(0) << "Note: config files must not begin with '-' !" << endl; 1436 1436 return (1); 1437 1437 break; 1438 1438 case 'v': 1439 while (argv[argptr-1][verbosity+1] == 'v') { 1440 verbosity++; 1441 } 1442 setVerbosity(verbosity); 1443 Log() << Verbose(0) << "Setting verbosity to " << verbosity << "." << endl; 1444 break; 1439 1445 case 'V': 1440 1446 Log() << Verbose(0) << argv[0] << " " << VERSIONSTRING << endl; … … 1841 1847 } 1842 1848 } else { 1843 eLog() << Verbose( 0) << "Removal failed due to missing atoms on molecule or wrong id." << endl;1849 eLog() << Verbose(1) << "Removal failed due to missing atoms on molecule or wrong id." << endl; 1844 1850 } 1845 1851 argptr+=2; … … 2002 2008 j = atoi(argv[argptr++]); 2003 2009 if ((j<0) || (j>1)) { 2004 eLog() << Verbose(1) << " ERROR:Argument of '-m' should be either 0 for no-rotate or 1 for rotate." << endl;2010 eLog() << Verbose(1) << "Argument of '-m' should be either 0 for no-rotate or 1 for rotate." << endl; 2005 2011 j = 0; 2006 2012 } … … 2064 2070 // repetition[i] = atoi(argv[argptr++]); 2065 2071 // if (repetition[i] < 1) 2066 // eLog() << Verbose( 0) << "ERROR:repetition value must be greater 1!" << endl;2072 // eLog() << Verbose(1) << "repetition value must be greater 1!" << endl; 2067 2073 // repetition[i] = 1; 2068 2074 // } … … 2083 2089 Vector ** vectors; 2084 2090 if (faktor < 1) { 2085 eLog() << Verbose( 0) << "ERROR: Repetition faktor mus be greater than 1!" << endl;2091 eLog() << Verbose(1) << "Repetition factor mus be greater than 1!" << endl; 2086 2092 faktor = 1; 2087 2093 } … … 2100 2106 } 2101 2107 if (count != j) 2102 Log() << Verbose(0) << "ERROR:AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl;2108 eLog() << Verbose(1) << "AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl; 2103 2109 x.Zero(); 2104 2110 y.Zero(); … … 2165 2171 int j; 2166 2172 2167 setVerbosity( 2);2173 setVerbosity(0); 2168 2174 2169 2175 // =========================== PARSE COMMAND LINE OPTIONS ==================================== -
src/config.cpp
rf1ef60a r717e0c 72 72 file= new ifstream(filename); 73 73 if (file == NULL) { 74 eLog() << Verbose( 0) << "ERROR:config file " << filename << " missing!" << endl;74 eLog() << Verbose(1) << "config file " << filename << " missing!" << endl; 75 75 return; 76 76 } … … 87 87 // allocate buffer's 1st dimension 88 88 if (buffer != NULL) { 89 eLog() << Verbose( 0) << "ERROR:FileBuffer->buffer is not NULL!" << endl;89 eLog() << Verbose(1) << "FileBuffer->buffer is not NULL!" << endl; 90 90 return; 91 91 } else … … 653 653 { 654 654 if (FileBuffer != NULL) { 655 eLog() << Verbose( 1) << "WARNING:deleting present FileBuffer in PrepareFileBuffer()." << endl;655 eLog() << Verbose(2) << "deleting present FileBuffer in PrepareFileBuffer()." << endl; 656 656 delete(FileBuffer); 657 657 } … … 847 847 ifstream *file = new ifstream(filename); 848 848 if (file == NULL) { 849 eLog() << Verbose( 0) << "ERROR:config file " << filename << " missing!" << endl;849 eLog() << Verbose(1) << "config file " << filename << " missing!" << endl; 850 850 return; 851 851 } … … 1080 1080 ifstream *file = new ifstream(filename); 1081 1081 if (file == NULL) { 1082 eLog() << Verbose( 0) << "ERROR:config file " << filename << " missing!" << endl;1082 eLog() << Verbose(1) << "config file " << filename << " missing!" << endl; 1083 1083 return; 1084 1084 } -
src/ellipsoid.cpp
rf1ef60a r717e0c 241 241 x = new Vector[PointsToPick]; 242 242 } else { 243 eLog() << Verbose(2) << " WARNING:Given pointer to vector array seems already allocated." << endl;243 eLog() << Verbose(2) << "Given pointer to vector array seems already allocated." << endl; 244 244 } 245 245 … … 341 341 x = new Vector[PointsToPick]; 342 342 } else { 343 eLog() << Verbose(2) << " WARNING:Given pointer to vector array seems already allocated." << endl;343 eLog() << Verbose(2) << "Given pointer to vector array seems already allocated." << endl; 344 344 } 345 345 -
src/linkedcell.cpp
rf1ef60a r717e0c 47 47 Log() << Verbose(1) << "Begin of LinkedCell" << endl; 48 48 if (set->IsEmpty()) { 49 eLog() << Verbose( 0) << "ERROR:set contains no linked cell nodes!" << endl;49 eLog() << Verbose(1) << "set contains no linked cell nodes!" << endl; 50 50 return; 51 51 } … … 79 79 Log() << Verbose(2) << "Allocating cells ... "; 80 80 if (LC != NULL) { 81 Log() << Verbose(1) << "ERROR:Linked Cell list is already allocated, I do nothing." << endl;81 eLog() << Verbose(1) << "Linked Cell list is already allocated, I do nothing." << endl; 82 82 return; 83 83 } … … 122 122 Log() << Verbose(1) << "Begin of LinkedCell" << endl; 123 123 if (set->empty()) { 124 eLog() << Verbose( 0) << "ERROR:set contains no linked cell nodes!" << endl;124 eLog() << Verbose(1) << "set contains no linked cell nodes!" << endl; 125 125 return; 126 126 } … … 151 151 Log() << Verbose(2) << "Allocating cells ... "; 152 152 if (LC != NULL) { 153 Log() << Verbose(1) << "ERROR:Linked Cell list is already allocated, I do nothing." << endl;153 eLog() << Verbose(1) << "Linked Cell list is already allocated, I do nothing." << endl; 154 154 return; 155 155 } … … 199 199 status = status && ((n[i] >=0) && (n[i] < N[i])); 200 200 if (!status) 201 eLog() << Verbose( 0) << "ERROR:indices are out of bounds!" << endl;201 eLog() << Verbose(1) << "indices are out of bounds!" << endl; 202 202 return status; 203 203 }; … … 260 260 return status; 261 261 } else { 262 eLog() << Verbose(1) << " ERROR:Node at " << *Walker << " is out of bounds." << endl;262 eLog() << Verbose(1) << "Node at " << *Walker << " is out of bounds." << endl; 263 263 return false; 264 264 } -
src/memoryusageobserver.cpp
rf1ef60a r717e0c 88 88 89 89 if (current == memoryUsers.end()) { 90 Log() << Verbose(0) << "WARNING:There is non-tracked memory to be freed. Pointer "90 eLog() << Verbose(2) << "There is non-tracked memory to be freed. Pointer " 91 91 << pointer << " is not registered by MemoryUsageObserver: "; 92 92 if (msg != NULL) -
src/molecule.cpp
rf1ef60a r717e0c 192 192 BondRescale = TopOrigin->type->HBondDistance[TopBond->BondDegree-1]; 193 193 if (BondRescale == -1) { 194 eLog() << Verbose( 3) << "ERROR:There is no typical hydrogen bond distance in replacing bond (" << TopOrigin->Name << "<->" << TopReplacement->Name << ") of degree " << TopBond->BondDegree << "!" << endl;194 eLog() << Verbose(1) << "There is no typical hydrogen bond distance in replacing bond (" << TopOrigin->Name << "<->" << TopReplacement->Name << ") of degree " << TopBond->BondDegree << "!" << endl; 195 195 return false; 196 196 BondRescale = bondlength; … … 235 235 SecondOtherAtom = (*Runner)->GetOtherAtom(TopOrigin); 236 236 } else { 237 Log() << Verbose(3) << "WARNING:Detected more than four bonds for atom " << TopOrigin->Name;237 eLog() << Verbose(2) << "Detected more than four bonds for atom " << TopOrigin->Name; 238 238 } 239 239 } … … 272 272 bondangle = TopOrigin->type->HBondAngle[1]; 273 273 if (bondangle == -1) { 274 Log() << Verbose(3) << "ERROR:There is no typical hydrogen bond angle in replacing bond (" << TopOrigin->Name << "<->" << TopReplacement->Name << ") of degree " << TopBond->BondDegree << "!" << endl;274 eLog() << Verbose(1) << "There is no typical hydrogen bond angle in replacing bond (" << TopOrigin->Name << "<->" << TopReplacement->Name << ") of degree " << TopBond->BondDegree << "!" << endl; 275 275 return false; 276 276 bondangle = 0; … … 394 394 break; 395 395 default: 396 eLog() << Verbose( 0) << "ERROR:BondDegree does not state single, double or triple bond!" << endl;396 eLog() << Verbose(1) << "BondDegree does not state single, double or triple bond!" << endl; 397 397 AllWentWell = false; 398 398 break; … … 541 541 add(Binder, last); 542 542 } else { 543 eLog() << Verbose(1) << " ERROR:Could not add bond between " << atom1->Name << " and " << atom2->Name << " as one or both are not present in the molecule." << endl;543 eLog() << Verbose(1) << "Could not add bond between " << atom1->Name << " and " << atom2->Name << " as one or both are not present in the molecule." << endl; 544 544 } 545 545 return Binder; … … 619 619 AtomCount--; 620 620 } else 621 eLog() << Verbose( 0) << "ERROR:Atom " << pointer->Name << " is of element " << pointer->type->Z << " but the entry in the table of the molecule is 0!" << endl;621 eLog() << Verbose(1) << "Atom " << pointer->Name << " is of element " << pointer->type->Z << " but the entry in the table of the molecule is 0!" << endl; 622 622 if (ElementsInMolecule[pointer->type->Z] == 0) // was last atom of this element? 623 623 ElementCount--; … … 637 637 ElementsInMolecule[pointer->type->Z]--; // decrease number of atom of this element 638 638 else 639 eLog() << Verbose( 0) << "ERROR:Atom " << pointer->Name << " is of element " << pointer->type->Z << " but the entry in the table of the molecule is 0!" << endl;639 eLog() << Verbose(1) << "Atom " << pointer->Name << " is of element " << pointer->type->Z << " but the entry in the table of the molecule is 0!" << endl; 640 640 if (ElementsInMolecule[pointer->type->Z] == 0) // was last atom of this element? 641 641 ElementCount--; -
src/molecule_dynamics.cpp
rf1ef60a r717e0c 428 428 } 429 429 if (Potential > Params.PenaltyConstants[2]) { 430 eLog() << Verbose( 0) << "ERROR:The two-step permutation procedure did not maintain injectivity!" << endl;430 eLog() << Verbose(1) << "The two-step permutation procedure did not maintain injectivity!" << endl; 431 431 exit(255); 432 432 } 433 433 //Log() << Verbose(0) << endl; 434 434 } else { 435 eLog() << Verbose( 0) << "ERROR: "<< *Runner << " was not the owner of " << *Sprinter << "!" << endl;435 eLog() << Verbose(1) << *Runner << " was not the owner of " << *Sprinter << "!" << endl; 436 436 exit(255); 437 437 } -
src/molecule_fragmentation.cpp
rf1ef60a r717e0c 441 441 // transmorph graph keyset list into indexed KeySetList 442 442 if (GlobalKeySetList == NULL) { 443 Log() << Verbose(1) << "ERROR:Given global key set list (graph) is NULL!" << endl;443 eLog() << Verbose(1) << "Given global key set list (graph) is NULL!" << endl; 444 444 return false; 445 445 } … … 898 898 } 899 899 } else { 900 Log() << Verbose(0) << "ERROR:Son " << Runner->Name << " has father " << FatherOfRunner->Name << " but its entry in SonList is " << SonList[FatherOfRunner->nr] << "!" << endl;900 eLog() << Verbose(1) << "Son " << Runner->Name << " has father " << FatherOfRunner->Name << " but its entry in SonList is " << SonList[FatherOfRunner->nr] << "!" << endl; 901 901 } 902 902 if ((LonelyFlag) && (Leaf->AtomCount > 1)) { … … 1133 1133 Log() << Verbose(0) << endl; 1134 1134 //if (!CheckForConnectedSubgraph(FragmentSearch->FragmentSet)) 1135 // Log() << Verbose(0) << "ERROR:The found fragment is not a connected subgraph!" << endl;1135 //eLog() << Verbose(1) << "The found fragment is not a connected subgraph!" << endl; 1136 1136 InsertFragmentIntoGraph(FragmentSearch); 1137 1137 } -
src/moleculelist.cpp
rf1ef60a r717e0c 311 311 Tesselation *TesselStruct = NULL; 312 312 if ((srcmol == NULL) || (mol == NULL)) { 313 Log() << Verbose(1) << "ERROR:Either fixed or variable molecule is given as NULL." << endl;313 eLog() << Verbose(1) << "Either fixed or variable molecule is given as NULL." << endl; 314 314 return false; 315 315 } … … 319 319 FindNonConvexBorder(mol, TesselStruct, (const LinkedCell *&)LCList, 4., NULL); 320 320 if (TesselStruct == NULL) { 321 Log() << Verbose(1) << "ERROR:Could not tesselate the fixed molecule." << endl;321 eLog() << Verbose(1) << "Could not tesselate the fixed molecule." << endl; 322 322 return false; 323 323 } -
src/parser.cpp
rf1ef60a r717e0c 226 226 } 227 227 } else { 228 eLog() << Verbose( 0) << "ERROR:Matrix nr. " << MatrixNr << " has column and row count of (" << ColumnCounter[MatrixNr] << "," << RowCounter[MatrixNr] << "), could not allocate nor parse!" << endl;228 eLog() << Verbose(1) << "Matrix nr. " << MatrixNr << " has column and row count of (" << ColumnCounter[MatrixNr] << "," << RowCounter[MatrixNr] << "), could not allocate nor parse!" << endl; 229 229 } 230 230 input.close(); -
src/periodentafel.cpp
rf1ef60a r717e0c 313 313 314 314 if (!otherstatus) 315 eLog() << Verbose( 0) << "WARNING:Something went wrong while parsing the other databases!" << endl;315 eLog() << Verbose(2) << "Something went wrong while parsing the other databases!" << endl; 316 316 317 317 return status; -
src/stackclass.hpp
rf1ef60a r717e0c 72 72 return true; 73 73 } else { 74 eLog() << Verbose( 0) << "ERROR:Stack is full, " << "Stack: CurrentLastEntry " << CurrentLastEntry<< "\tCurrentFirstEntry " << CurrentFirstEntry << "\tNextFreeField " << NextFreeField << "\tEntryCount " << EntryCount << "!" << endl;74 eLog() << Verbose(1) << "Stack is full, " << "Stack: CurrentLastEntry " << CurrentLastEntry<< "\tCurrentFirstEntry " << CurrentFirstEntry << "\tNextFreeField " << NextFreeField << "\tEntryCount " << EntryCount << "!" << endl; 75 75 return false; 76 76 } … … 87 87 Walker = StackList[CurrentFirstEntry]; 88 88 if (Walker == NULL) 89 eLog() << Verbose( 0) << "ERROR:Stack's field is empty!" << endl;89 eLog() << Verbose(1) << "Stack's field is empty!" << endl; 90 90 StackList[CurrentFirstEntry] = NULL; 91 91 if (CurrentFirstEntry != CurrentLastEntry) { // hasn't last item been popped as well? … … 96 96 } 97 97 } else 98 eLog() << Verbose( 0) << "ERROR:Stack is empty!" << endl;98 eLog() << Verbose(1) << "Stack is empty!" << endl; 99 99 return Walker; 100 100 }; … … 111 111 StackList[CurrentLastEntry] = NULL; 112 112 if (Walker == NULL) 113 eLog() << Verbose( 0) << "ERROR:Stack's field is empty!" << endl;113 eLog() << Verbose(1) << "Stack's field is empty!" << endl; 114 114 NextFreeField = CurrentLastEntry; 115 115 if (CurrentLastEntry != CurrentFirstEntry) // has there been more than one item on stack 116 116 CurrentLastEntry = (CurrentLastEntry + (EntryCount-1)) % EntryCount; // step back from current free field to last (modulo does not work in -1, thus go EntryCount-1 instead) 117 117 } else { 118 eLog() << Verbose( 0) << "ERROR:Stack is empty!" << endl;118 eLog() << Verbose(1) << "Stack is empty!" << endl; 119 119 } 120 120 return Walker; … … 151 151 } while (i!=NextFreeField); 152 152 else 153 eLog() << Verbose( 0) << "ERROR:Stack is already empty!" << endl;153 eLog() << Verbose(1) << "Stack is already empty!" << endl; 154 154 if (found) { 155 155 NextFreeField = CurrentLastEntry; -
src/tesselation.cpp
rf1ef60a r717e0c 48 48 //Log() << Verbose(5) << "Erasing point nr. " << Nr << "." << endl; 49 49 if (!lines.empty()) 50 eLog() << Verbose( 0) << "WARNING:Memory Leak! I " << *this << " am still connected to some lines." << endl;50 eLog() << Verbose(2) << "Memory Leak! I " << *this << " am still connected to some lines." << endl; 51 51 node = NULL; 52 52 }; … … 153 153 } 154 154 if (!triangles.empty()) 155 eLog() << Verbose( 0) << "WARNING:Memory Leak! I " << *this << " am still connected to some triangles." << endl;155 eLog() << Verbose(2) << "Memory Leak! I " << *this << " am still connected to some triangles." << endl; 156 156 }; 157 157 … … 188 188 // get the two triangles 189 189 if (triangles.size() != 2) { 190 Log() << Verbose(1) << "ERROR:Baseline " << *this << " is connected to less than two triangles, Tesselation incomplete!" << endl;190 eLog() << Verbose(1) << "Baseline " << *this << " is connected to less than two triangles, Tesselation incomplete!" << endl; 191 191 return true; 192 192 } … … 226 226 i++; 227 227 } else { 228 // Log() << Verbose(2) << "ERROR:I cannot find third node in triangle, something's wrong." << endl;228 //eLog() << Verbose(1) << "I cannot find third node in triangle, something's wrong." << endl; 229 229 return true; 230 230 } … … 409 409 } while (CrossPoint.NormSquared() < MYEPSILON); 410 410 if (i==3) { 411 Log() << Verbose(1) << "ERROR:Could not find any cross points, something's utterly wrong here!" << endl;411 eLog() << Verbose(1) << "Could not find any cross points, something's utterly wrong here!" << endl; 412 412 exit(255); 413 413 } … … 634 634 runner->second = NULL; 635 635 } else 636 eLog() << Verbose(1) << " ERROR:The triangle " << runner->first << " has already been free'd." << endl;636 eLog() << Verbose(1) << "The triangle " << runner->first << " has already been free'd." << endl; 637 637 } 638 638 Log() << Verbose(1) << "This envelope was written to file " << TriangleFilesWritten << " times(s)." << endl; … … 1174 1174 if (NewLines[j]->IsConnectedTo(BLS[0])) { 1175 1175 if (n>2) { 1176 Log() << Verbose(1) << "ERROR: " <<BLS[0] << " connects to all of the new lines?!" << endl;1176 Log() << Verbose(1) << BLS[0] << " connects to all of the new lines?!" << endl; 1177 1177 return false; 1178 1178 } else … … 1191 1191 } 1192 1192 } else { // something is wrong with FindClosestTriangleToPoint! 1193 Log() << Verbose(1) << "ERROR:The closest triangle did not produce an intersection!" << endl;1193 eLog() << Verbose(1) << "The closest triangle did not produce an intersection!" << endl; 1194 1194 return false; 1195 1195 } … … 1379 1379 triangle->lines[i] = NULL; // free'd or not: disconnect 1380 1380 } else 1381 eLog() << Verbose( 0) << "ERROR:This line " << i << " has already been free'd." << endl;1381 eLog() << Verbose(1) << "This line " << i << " has already been free'd." << endl; 1382 1382 } 1383 1383 … … 1432 1432 line->endpoints[i] = NULL; // free'd or not: disconnect 1433 1433 } else 1434 eLog() << Verbose( 0) << "ERROR:Endpoint " << i << " has already been free'd." << endl;1434 eLog() << Verbose(1) << "Endpoint " << i << " has already been free'd." << endl; 1435 1435 } 1436 1436 if (!line->triangles.empty()) 1437 eLog() << Verbose( 0) << "WARNING:Memory Leak! I " << *line << " am still connected to some triangles." << endl;1437 eLog() << Verbose(2) << "Memory Leak! I " << *line << " am still connected to some triangles." << endl; 1438 1438 1439 1439 if (LinesOnBoundary.erase(line->Nr)) … … 1604 1604 } 1605 1605 } else { 1606 eLog() << Verbose( 0) << "ERROR:The current cell " << LC->n[0] << "," << LC->n[1] << "," << LC->n[2] << " is invalid!" << endl;1606 eLog() << Verbose(1) << "The current cell " << LC->n[0] << "," << LC->n[1] << "," << LC->n[2] << " is invalid!" << endl; 1607 1607 } 1608 1608 } … … 1786 1786 if (fabs(OldSphereCenter.ScalarProduct(&SearchDirection)) > HULLEPSILON) { 1787 1787 // rotated the wrong way! 1788 eLog() << Verbose( 0) << "ERROR:SearchDirection and RelativeOldSphereCenter are still not orthogonal!" << endl;1788 eLog() << Verbose(1) << "SearchDirection and RelativeOldSphereCenter are still not orthogonal!" << endl; 1789 1789 } 1790 1790 … … 1924 1924 if (fabs(OldSphereCenter.ScalarProduct(&SearchDirection)) > HULLEPSILON) { 1925 1925 // rotated the wrong way! 1926 eLog() << Verbose( 0) << "ERROR:SearchDirection and RelativeOldSphereCenter are still not orthogonal!" << endl;1926 eLog() << Verbose(1) << "SearchDirection and RelativeOldSphereCenter are still not orthogonal!" << endl; 1927 1927 } 1928 1928 … … 1935 1935 1936 1936 if (OptCandidates->begin() == OptCandidates->end()) { 1937 eLog() << Verbose( 0) << "WARNING:Could not find a suitable candidate." << endl;1937 eLog() << Verbose(2) << "Could not find a suitable candidate." << endl; 1938 1938 return false; 1939 1939 } … … 1983 1983 //Log() << Verbose(1) << "We have "<< TrianglesOnBoundaryCount << " for line " << *BaseRay << "." << endl; 1984 1984 } else { 1985 Log() << Verbose(1) << "WARNING:This triangle consisting of ";1985 eLog() << Verbose(2) << "This triangle consisting of "; 1986 1986 Log() << Verbose(0) << *(*it)->point << ", "; 1987 1987 Log() << Verbose(0) << *BaseRay->endpoints[0]->node << " and "; … … 2009 2009 (*it)->OtherOptCenter.Scale(-1.); 2010 2010 2011 Log() << Verbose(0) << "--> WARNING: Special new triangle with " << *BTS << " and normal vector " << BTS->NormalVector 2012 << " for this triangle ... " << endl; 2011 eLog() << Verbose(2) << "--> WARNING: Special new triangle with " << *BTS << " and normal vector " << BTS->NormalVector << " for this triangle ... " << endl; 2013 2012 Log() << Verbose(1) << "We have "<< BaseRay->triangles.size() << " for line " << BaseRay << "." << endl; 2014 2013 } else { 2015 Log() << Verbose(1) << "WARNING: This triangle consisting of "; 2016 Log() << Verbose(0) << *(*it)->point << ", "; 2017 Log() << Verbose(0) << *BaseRay->endpoints[0]->node << " and "; 2018 Log() << Verbose(0) << *BaseRay->endpoints[1]->node << " "; 2019 Log() << Verbose(0) << "exists and is not added, as it does not seem helpful!" << endl; 2014 eLog() << Verbose(2) << "This triangle consisting of " << *(*it)->point << ", " << *BaseRay->endpoints[0]->node << " and " << *BaseRay->endpoints[1]->node << " " << "exists and is not added, as it does not seem helpful!" << endl; 2020 2015 result = false; 2021 2016 } … … 2033 2028 BaseRay = BLS[0]; 2034 2029 if ((BTS != NULL) && (BTS->NormalVector.NormSquared() < MYEPSILON)) { 2035 Log() << Verbose(1) << "CRITICAL:Triangle " << *BTS << " has zero normal vector!" << endl;2030 eLog() << Verbose(1) << "Triangle " << *BTS << " has zero normal vector!" << endl; 2036 2031 exit(255); 2037 2032 } … … 2174 2169 BaseLineNormal.Zero(); 2175 2170 if (Base->triangles.size() < 2) { 2176 Log() << Verbose(2) << "ERROR:Less than two triangles are attached to this baseline!" << endl;2171 eLog() << Verbose(1) << "Less than two triangles are attached to this baseline!" << endl; 2177 2172 return 0.; 2178 2173 } … … 2214 2209 BaseLineNormal.Zero(); 2215 2210 if (Base->triangles.size() < 2) { 2216 Log() << Verbose(2) << "ERROR:Less than two triangles are attached to this baseline!" << endl;2211 eLog() << Verbose(1) << "Less than two triangles are attached to this baseline!" << endl; 2217 2212 return NULL; 2218 2213 } … … 2250 2245 // check whether everything is in place to create new lines and triangles 2251 2246 if (i<4) { 2252 Log() << Verbose(1) << "ERROR:We have not gathered enough baselines!" << endl;2247 eLog() << Verbose(1) << "We have not gathered enough baselines!" << endl; 2253 2248 return NULL; 2254 2249 } 2255 2250 for (int j=0;j<4;j++) 2256 2251 if (OldLines[j] == NULL) { 2257 Log() << Verbose(1) << "ERROR:We have not gathered enough baselines!" << endl;2252 eLog() << Verbose(1) << "We have not gathered enough baselines!" << endl; 2258 2253 return NULL; 2259 2254 } 2260 2255 for (int j=0;j<2;j++) 2261 2256 if (OldPoints[j] == NULL) { 2262 Log() << Verbose(1) << "ERROR:We have not gathered enough endpoints!" << endl;2257 eLog() << Verbose(1) << "We have not gathered enough endpoints!" << endl; 2263 2258 return NULL; 2264 2259 } … … 2336 2331 N[i] = LC->n[i]; 2337 2332 } else { 2338 eLog() << Verbose( 0) << "ERROR:Point " << *a << " is not found in cell " << LC->index << "." << endl;2333 eLog() << Verbose(1) << "Point " << *a << " is not found in cell " << LC->index << "." << endl; 2339 2334 return; 2340 2335 } … … 2475 2470 // test whether old center is on the band's plane 2476 2471 if (fabs(OldSphereCenter.ScalarProduct(&CirclePlaneNormal)) > HULLEPSILON) { 2477 eLog() << Verbose( 0) << "ERROR:Something's very wrong here: OldSphereCenter is not on the band's plane as desired by " << fabs(OldSphereCenter.ScalarProduct(&CirclePlaneNormal)) << "!" << endl;2472 eLog() << Verbose(1) << "Something's very wrong here: OldSphereCenter is not on the band's plane as desired by " << fabs(OldSphereCenter.ScalarProduct(&CirclePlaneNormal)) << "!" << endl; 2478 2473 OldSphereCenter.ProjectOntoPlane(&CirclePlaneNormal); 2479 2474 } … … 2485 2480 //Log() << Verbose(2) << "INFO: SearchDirection is " << SearchDirection << "." << endl; 2486 2481 if (fabs(OldSphereCenter.ScalarProduct(&SearchDirection)) > HULLEPSILON) { // rotated the wrong way! 2487 eLog() << Verbose( 0) << "ERROR:SearchDirection and RelativeOldSphereCenter are not orthogonal!" << endl;2482 eLog() << Verbose(1) << "SearchDirection and RelativeOldSphereCenter are not orthogonal!" << endl; 2488 2483 } 2489 2484 … … 2494 2489 //Log() << Verbose(2) << "INFO: Center cell is " << N[0] << ", " << N[1] << ", " << N[2] << " with No. " << LC->index << "." << endl; 2495 2490 } else { 2496 eLog() << Verbose( 0) << "ERROR:Vector " << CircleCenter << " is outside of LinkedCell's bounding box." << endl;2491 eLog() << Verbose(1) << "Vector " << CircleCenter << " is outside of LinkedCell's bounding box." << endl; 2497 2492 return; 2498 2493 } … … 2599 2594 } 2600 2595 } else { 2601 eLog() << Verbose( 2) << "ERROR:The projected center of the old sphere has radius " << radius << " instead of " << CircleRadius << "." << endl;2596 eLog() << Verbose(1) << "The projected center of the old sphere has radius " << radius << " instead of " << CircleRadius << "." << endl; 2602 2597 } 2603 2598 } else { … … 2690 2685 triangles->unique(); 2691 2686 } else { 2692 Log() << Verbose(1) << "ERROR:I cannot find a boundary point to the tessel point " << *trianglePoints[0] << "." << endl;2687 eLog() << Verbose(1) << "I cannot find a boundary point to the tessel point " << *trianglePoints[0] << "." << endl; 2693 2688 return NULL; 2694 2689 } … … 2701 2696 for (int i=0;i<3;i++) { 2702 2697 if (trianglePoints[i] == NULL) { 2703 Log() << Verbose(1) << "ERROR:IsInnerPoint encounters serious error, point " << i << " not found." << endl;2698 eLog() << Verbose(1) << "IsInnerPoint encounters serious error, point " << i << " not found." << endl; 2704 2699 } 2705 2700 //Log() << Verbose(2) << "List of triangle points:" << endl; … … 2720 2715 2721 2716 if ((triangles == NULL) || (triangles->empty())) { 2722 Log() << Verbose(0) << "ERROR:There is no nearest triangle. Please check the tesselation structure.";2717 eLog() << Verbose(1) << "There is no nearest triangle. Please check the tesselation structure."; 2723 2718 delete(triangles); 2724 2719 return NULL; … … 2754 2749 Log() << Verbose(2) << "Normal Vector of this back side is " << result->NormalVector << "." << endl; 2755 2750 if (Center.ScalarProduct(&result->NormalVector) < 0) { 2756 Log() << Verbose(1) << "ERROR:Front and back side yield NormalVector in wrong direction!" << endl;2751 eLog() << Verbose(1) << "Front and back side yield NormalVector in wrong direction!" << endl; 2757 2752 } 2758 2753 } … … 2856 2851 2857 2852 if (connectedPoints->size() == 0) { // if have not found any points 2858 Log() << Verbose(1) << "ERROR:We have not found any connected points to " << *Point<< "." << endl;2853 eLog() << Verbose(1) << "We have not found any connected points to " << *Point<< "." << endl; 2859 2854 return NULL; 2860 2855 } … … 2979 2974 ReferencePoint = PointRunner->second; 2980 2975 } else { 2981 Log() << Verbose(2) << "ERROR:GetPathOfConnectedPoints() could not find the BoundaryPoint belonging to " << *Point << "." << endl;2976 eLog() << Verbose(1) << "GetPathOfConnectedPoints() could not find the BoundaryPoint belonging to " << *Point << "." << endl; 2982 2977 return NULL; 2983 2978 } … … 2996 2991 LineRunner = TouchedLine.find(runner->second); 2997 2992 if (LineRunner == TouchedLine.end()) { 2998 Log() << Verbose(2) << "ERROR:I could not find " << *runner->second << " in the touched list." << endl;2993 eLog() << Verbose(1) << "I could not find " << *runner->second << " in the touched list." << endl; 2999 2994 } else if (!LineRunner->second) { 3000 2995 LineRunner->second = true; … … 3026 3021 } 3027 3022 } else { 3028 Log() << Verbose(2) << "ERROR:I could not find " << *triangle << " in the touched list." << endl;3023 eLog() << Verbose(1) << "I could not find " << *triangle << " in the touched list." << endl; 3029 3024 triangle = NULL; 3030 3025 } … … 3043 3038 LineRunner = TouchedLine.find(CurrentLine); 3044 3039 if (LineRunner == TouchedLine.end()) 3045 Log() << Verbose(2) << "ERROR:I could not find " << *CurrentLine << " in the touched list." << endl;3040 eLog() << Verbose(1) << "I could not find " << *CurrentLine << " in the touched list." << endl; 3046 3041 else 3047 3042 LineRunner->second = true; … … 3061 3056 } 3062 3057 } else { 3063 Log() << Verbose(1) << "ERROR:There are no lines attached to " << *ReferencePoint << "." << endl;3058 eLog() << Verbose(1) << "There are no lines attached to " << *ReferencePoint << "." << endl; 3064 3059 } 3065 3060 … … 3139 3134 3140 3135 if (Point == NULL) { 3141 Log() << Verbose(1) << "ERROR:Point given is NULL." << endl;3136 eLog() << Verbose(1) << "Point given is NULL." << endl; 3142 3137 } else { 3143 3138 // go through its lines and insert all triangles … … 3171 3166 3172 3167 if (point == NULL) { 3173 Log() << Verbose(1) << "ERROR:Cannot remove the point " << point << ", it's NULL!" << endl;3168 eLog() << Verbose(1) << "Cannot remove the point " << point << ", it's NULL!" << endl; 3174 3169 return 0.; 3175 3170 } else … … 3181 3176 // get list of connected points 3182 3177 if (point->lines.empty()) { 3183 Log() << Verbose(1) << "ERROR:Cannot remove the point " << *point << ", it's connected to no lines!" << endl;3178 eLog() << Verbose(1) << "Cannot remove the point " << *point << ", it's connected to no lines!" << endl; 3184 3179 return 0.; 3185 3180 } … … 3281 3276 triangle = GetPresentTriangle(TriangleCandidates); 3282 3277 if (triangle != NULL) { 3283 Log() << Verbose(1) << "WARNING:New triangle already present, skipping!" << endl;3278 eLog() << Verbose(2) << "New triangle already present, skipping!" << endl; 3284 3279 StartNode++; 3285 3280 MiddleNode++; … … 3619 3614 NearestBoundaryPoint = PointRunner->second; 3620 3615 } else { 3621 Log() << Verbose(1) << "ERROR:I cannot find the boundary point." << endl;3616 eLog() << Verbose(1) << "I cannot find the boundary point." << endl; 3622 3617 return; 3623 3618 } -
src/tesselationhelpers.cpp
rf1ef60a r717e0c 77 77 78 78 if (fabs(m11) < MYEPSILON) 79 eLog() << Verbose( 0) << "ERROR:three points are colinear." << endl;79 eLog() << Verbose(1) << "three points are colinear." << endl; 80 80 81 81 center->x[0] = 0.5 * m12/ m11; … … 84 84 85 85 if (fabs(a.Distance(center) - RADIUS) > MYEPSILON) 86 eLog() << Verbose( 0) << "ERROR:The given center is further way by " << fabs(a.Distance(center) - RADIUS) << " from a than RADIUS." << endl;86 eLog() << Verbose(1) << "The given center is further way by " << fabs(a.Distance(center) - RADIUS) << " from a than RADIUS." << endl; 87 87 88 88 gsl_matrix_free(A); … … 224 224 // test whether new center is on the parameter circle's plane 225 225 if (fabs(helper.ScalarProduct(&CirclePlaneNormal)) > HULLEPSILON) { 226 eLog() << Verbose( 0) << "ERROR:Something's very wrong here: NewSphereCenter is not on the band's plane as desired by " <<fabs(helper.ScalarProduct(&CirclePlaneNormal)) << "!" << endl;226 eLog() << Verbose(1) << "Something's very wrong here: NewSphereCenter is not on the band's plane as desired by " <<fabs(helper.ScalarProduct(&CirclePlaneNormal)) << "!" << endl; 227 227 helper.ProjectOntoPlane(&CirclePlaneNormal); 228 228 } … … 230 230 // test whether the new center vector has length of CircleRadius 231 231 if (fabs(radius - CircleRadius) > HULLEPSILON) 232 eLog() << Verbose(1) << " ERROR:The projected center of the new sphere has radius " << radius << " instead of " << CircleRadius << "." << endl;232 eLog() << Verbose(1) << "The projected center of the new sphere has radius " << radius << " instead of " << CircleRadius << "." << endl; 233 233 alpha = helper.Angle(&OldSphereCenter); 234 234 // make the angle unique by checking the halfplanes/search direction … … 494 494 Vector BaseLineVector, OrthogonalVector, helper; 495 495 if (candidate1->BaseLine != candidate2->BaseLine) { // sanity check 496 Log() << Verbose(0) << "ERROR:sortCandidates was called for two different baselines: " << candidate1->BaseLine << " and " << candidate2->BaseLine << "." << endl;496 eLog() << Verbose(1) << "sortCandidates was called for two different baselines: " << candidate1->BaseLine << " and " << candidate2->BaseLine << "." << endl; 497 497 //return false; 498 498 exit(1); … … 577 577 } 578 578 } else { 579 eLog() << Verbose( 0) << "ERROR:The current cell " << LC->n[0] << "," << LC->n[1] << ","579 eLog() << Verbose(1) << "The current cell " << LC->n[0] << "," << LC->n[1] << "," 580 580 << LC->n[2] << " is invalid!" << endl; 581 581 } … … 633 633 } 634 634 } else { 635 eLog() << Verbose( 0) << "ERROR:The current cell " << LC->n[0] << "," << LC->n[1] << ","635 eLog() << Verbose(1) << "The current cell " << LC->n[0] << "," << LC->n[1] << "," 636 636 << LC->n[2] << " is invalid!" << endl; 637 637 } … … 741 741 } 742 742 } else { 743 eLog() << Verbose( 0) << "ERROR:Given vrmlfile is " << vrmlfile << "." << endl;743 eLog() << Verbose(1) << "Given vrmlfile is " << vrmlfile << "." << endl; 744 744 } 745 745 delete(center); … … 811 811 *rasterfile << "9\n# terminating special property\n"; 812 812 } else { 813 eLog() << Verbose( 0) << "ERROR:Given rasterfile is " << rasterfile << "." << endl;813 eLog() << Verbose(1) << "Given rasterfile is " << rasterfile << "." << endl; 814 814 } 815 815 IncludeSphereinRaster3D(rasterfile, Tess, cloud); -
src/vector.cpp
rf1ef60a r717e0c 229 229 factor = Direction.ScalarProduct(PlaneNormal); 230 230 if (factor < MYEPSILON) { // Uniqueness: line parallel to plane? 231 Log() << Verbose(2) << "WARNING:Line is parallel to plane, no intersection." << endl;231 eLog() << Verbose(2) << "Line is parallel to plane, no intersection." << endl; 232 232 return false; 233 233 } … … 253 253 return true; 254 254 } else { 255 Log() << Verbose(2) << "WARNING:Intersection point " << *this << " is not on plane." << endl;255 eLog() << Verbose(2) << "Intersection point " << *this << " is not on plane." << endl; 256 256 return false; 257 257 } … … 747 747 x[i] = C.x[i]; 748 748 } else { 749 eLog() << Verbose( 0) << "ERROR:inverse of matrix does not exists: det A = " << detA << "." << endl;749 eLog() << Verbose(1) << "inverse of matrix does not exists: det A = " << detA << "." << endl; 750 750 } 751 751 }; … … 799 799 x2.SubtractVector(y2); 800 800 if ((fabs(x1.Norm()) < MYEPSILON) || (fabs(x2.Norm()) < MYEPSILON) || (fabs(x1.Angle(&x2)) < MYEPSILON)) { 801 Log() << Verbose(4) << "WARNING:Given vectors are linear dependent." << endl;801 eLog() << Verbose(2) << "Given vectors are linear dependent." << endl; 802 802 return false; 803 803 } … … 833 833 Zero(); 834 834 if ((fabs(x1.Norm()) < MYEPSILON) || (fabs(x2.Norm()) < MYEPSILON) || (fabs(x1.Angle(&x2)) < MYEPSILON)) { 835 Log() << Verbose(4) << "WARNING:Given vectors are linear dependent." << endl;835 eLog() << Verbose(2) << "Given vectors are linear dependent." << endl; 836 836 return false; 837 837 }
Note:
See TracChangeset
for help on using the changeset viewer.