- Timestamp:
- May 1, 2013, 7:53:26 AM (12 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:
- 10f0cb
- Parents:
- 6ff62c
- git-author:
- Frederik Heber <heber@…> (03/27/13 12:07:15)
- git-committer:
- Frederik Heber <heber@…> (05/01/13 07:53:26)
- Location:
- src/Jobs
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Jobs/InterfaceVMGJob.cpp
r6ff62c rb6b21a 71 71 vmg_float _box_end, 72 72 const int& near_field_cells, 73 const bool _DoPrintDebug, 73 74 int coarseningSteps, 74 75 double alpha) : … … 78 79 sampled_input(_sampled_input), 79 80 returndata(_returndata), 80 level(levelMax) 81 level(levelMax), 82 DoPrintDebug(_DoPrintDebug) 81 83 { 82 84 for (size_t i=0;i<3;++i) { … … 144 146 comm.CommFromGhosts(particle_grid); 145 147 146 // print nuclei grid to vtk 147 comm.PrintGrid(particle_grid, "Sampled Nuclei Density"); 148 if (DoPrintDebug) { 149 // print nuclei grid to vtk 150 comm.PrintGrid(particle_grid, "Sampled Nuclei Density"); 151 } 148 152 149 153 // add sampled electron charge density onto grid … … 153 157 -1.); 154 158 155 // print electron grid to vtk 156 comm.PrintGrid(grid, "Sampled Electron Density"); 159 if (DoPrintDebug) { 160 // print electron grid to vtk 161 comm.PrintGrid(grid, "Sampled Electron Density"); 162 } 157 163 158 164 // add particle_grid onto grid … … 180 186 comm.PrintStringOnce("Grid charge integral: %e", charge_sum/(4.0 * VMG::Math::pi)); 181 187 182 // print total grid to vtk 183 comm.PrintGrid(grid, "Total Charge Density"); 188 if (DoPrintDebug) { 189 // print total grid to vtk 190 comm.PrintGrid(grid, "Total Charge Density"); 191 } 184 192 185 193 // delete temp_grid; … … 197 205 Index i; 198 206 199 // print output grid to vtk 200 comm.PrintGrid(grid, "Potential Solution"); 207 if (DoPrintDebug) { 208 // print output grid to vtk 209 comm.PrintGrid(grid, "Potential Solution"); 210 } 201 211 202 212 // obtain sampled potential from grid -
src/Jobs/InterfaceVMGJob.hpp
r6ff62c rb6b21a 46 46 vmg_float _box_end, 47 47 const int& near_field_cells, 48 const bool _DoPrintDebug=false, 48 49 int coarseningSteps=0, 49 double alpha=1.6); 50 double alpha=1.6 51 ); 50 52 51 53 ~InterfaceVMGJob() {} … … 71 73 //!> contains end of grid intervals 72 74 double box_end[3]; 75 76 //!> whether we do print grid for debug visualization or not 77 const bool DoPrintDebug; 73 78 }; 74 79 -
src/Jobs/VMGJob.cpp
r6ff62c rb6b21a 82 82 const std::vector< double > &_particle_charges, 83 83 const size_t _near_field_cells, 84 const size_t _interpolation_degree) : 84 const size_t _interpolation_degree, 85 const bool _DoPrintDebug) : 85 86 FragmentJob(_JobId), 86 87 density_grid(_density_grid), … … 90 91 interpolation_degree(_interpolation_degree), 91 92 returndata(static_cast<const SamplingGridProperties &>(_density_grid)), 92 particles() 93 particles(), 94 DoPrintDebug(_DoPrintDebug) 93 95 {} 94 96 … … 97 99 near_field_cells(3), 98 100 interpolation_degree(3), 99 particles() 101 particles(), 102 DoPrintDebug(false) 100 103 {} 101 104 … … 216 219 Vector(density_grid.begin), 217 220 density_grid.end[0]-density_grid.begin[0], 218 near_field_cells); 221 near_field_cells, 222 DoPrintDebug); 219 223 new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 220 224 new Givens<SolverSingular>(); -
src/Jobs/VMGJob.hpp
r6ff62c rb6b21a 42 42 * @param _interpolation_degree degree of interpolation polynomial for getting nuclei 43 43 * potential from grid 44 * @param _DoPrintDebug whether we do print grid for debug visualization or not 44 45 */ 45 46 VMGJob(const JobId_t _JobId, … … 48 49 const std::vector< double > &_particle_charges, 49 50 const size_t _near_field_cells, 50 const size_t _interpolation_degree 51 const size_t _interpolation_degree, 52 const bool _DoPrintDebug=false 51 53 ); 52 54 virtual ~VMGJob(); … … 70 72 //!> interpolation degree used in sampling the potential of the nuclei 71 73 const size_t interpolation_degree; 74 //!> whether we do print grid for debug visualization or not 75 const bool DoPrintDebug; 72 76 73 77 private: … … 120 124 ar & const_cast< size_t &>(near_field_cells); 121 125 ar & const_cast< size_t &>(interpolation_degree); 126 ar & const_cast< bool &>(DoPrintDebug); 122 127 ar & returndata; 123 128 }
Note:
See TracChangeset
for help on using the changeset viewer.