source: ThirdParty/vmg/src/solver/solver.cpp@ 33a694

Action_Thermostats Add_AtomRandomPerturbation Add_RotateAroundBondAction Add_SelectAtomByNameAction Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_StructOpt_integration_tests AutomationFragmentation_failures Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator Docu_Python_wait EmpiricalPotential_contain_HomologyGraph_documentation Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph Fix_ChronosMutex Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion GeometryObjects Gui_displays_atomic_force_velocity IndependentFragmentGrids_IntegrationTest JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks RotateToPrincipalAxisSystem_UndoRedo StoppableMakroAction TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps Ubuntu_1604_changes
Last change on this file since 33a694 was 7faa5c, checked in by Frederik Heber <heber@…>, 8 years ago

Merge commit 'de061d9d851257a04e924d4472df4523d33bb08b' as 'ThirdParty/vmg'

  • Property mode set to 100644
File size: 1.7 KB
Line 
1/*
2 * vmg - a versatile multigrid solver
3 * Copyright (C) 2012 Institute for Numerical Simulation, University of Bonn
4 *
5 * vmg is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * (at your option) any later version.
9 *
10 * vmg is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19/**
20 * @file solver.cpp
21 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
22 * @date Mon Apr 18 13:10:55 2011
23 *
24 * @brief VMG::Solver
25 *
26 */
27
28#ifdef HAVE_CONFIG_H
29#include <libvmg_config.h>
30#endif
31
32#include "solver/solver.hpp"
33
34using namespace VMG;
35
36void Solver::Run(Grid& sol, Grid& rhs)
37{
38#ifdef DEBUG_MATRIX_CHECKS
39 sol.IsConsistent();
40 rhs.IsConsistent();
41#endif
42
43#ifdef DEBUG
44 sol.IsCompatible(rhs);
45#endif
46
47 if (rhs.Global().LocalSize().Product() > 0) {
48 this->Realloc(rhs);
49 this->AssembleMatrix(rhs);
50 this->Compute();
51 this->ExportSol(sol, rhs);
52 }
53}
54
55void Solver::Realloc(int n)
56{
57 //Reallocate memory if necessary
58 this->size = n;
59
60 if (static_cast<int>(A.size()) < n*n)
61 A.resize(n*n);
62
63 if (static_cast<int>(b.size()) < n)
64 b.resize(n);
65
66 if (static_cast<int>(x.size()) < n)
67 x.resize(n);
68}
69
70void Solver::Realloc(Grid& sol)
71{
72 this->Realloc(sol.Global().GlobalSize().Product());
73}
Note: See TracBrowser for help on using the repository browser.