source: ThirdParty/vmg/src/discretization/discretization_poisson_fd.cpp@ 79b089

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.0 Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError 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 JobMarket_unresolvable_hostname_fix ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks RotateToPrincipalAxisSystem_UndoRedo StoppableMakroAction Subpackage_vmg TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps Ubuntu_1604_changes stable
Last change on this file since 79b089 was 7faa5c, checked in by Frederik Heber <heber@…>, 8 years ago

Merge commit 'de061d9d851257a04e924d4472df4523d33bb08b' as 'ThirdParty/vmg'
  • Property mode set to 100644
File size: 3.0 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 discretization_poisson_fd.cpp
21 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
22 * @date Mon Apr 18 13:03:47 2011
23 *
24 * @brief Discretization of the poisson equation
25 * using the standard 7-point stencil or
26 * the Collatz Mehrstellen Ansatz.
27 * Discretization error: O(h^2) or O(h^4)
28 *
29 */
30
31#ifdef HAVE_CONFIG_H
32#include <libvmg_config.h>
33#endif
34
35#include "discretization/discretization_poisson_fd.hpp"
36#include "grid/grid.hpp"
37#include "mg.hpp"
38
39using namespace VMG;
40
41void DiscretizationPoissonFD::InitDiscretizationPoissonFD()
42{
43 switch (order)
44 {
45 case 2:
46 stencil.SetDiag(6.0);
47 stencil.push_back(-1, 0, 0, -1.0);
48 stencil.push_back( 1, 0, 0, -1.0);
49 stencil.push_back( 0, -1, 0, -1.0);
50 stencil.push_back( 0, 1, 0, -1.0);
51 stencil.push_back( 0, 0, -1, -1.0);
52 stencil.push_back( 0, 0, 1, -1.0);
53 break;
54 case 4:
55 stencil.SetDiag(24.0/6.0);
56 stencil.push_back(-1, 0, 0, -2.0/6.0);
57 stencil.push_back( 1, 0, 0, -2.0/6.0);
58 stencil.push_back( 0, -1, 0, -2.0/6.0);
59 stencil.push_back( 0, 1, 0, -2.0/6.0);
60 stencil.push_back( 0, 0, -1, -2.0/6.0);
61 stencil.push_back( 0, 0, 1, -2.0/6.0);
62 stencil.push_back(-1, -1, 0, -1.0/6.0);
63 stencil.push_back(-1, 1, 0, -1.0/6.0);
64 stencil.push_back( 1, -1, 0, -1.0/6.0);
65 stencil.push_back( 1, 1, 0, -1.0/6.0);
66 stencil.push_back(-1, 0, -1, -1.0/6.0);
67 stencil.push_back(-1, 0, 1, -1.0/6.0);
68 stencil.push_back( 1, 0, -1, -1.0/6.0);
69 stencil.push_back( 1, 0, 1, -1.0/6.0);
70 stencil.push_back( 0, -1, -1, -1.0/6.0);
71 stencil.push_back( 0, -1, 1, -1.0/6.0);
72 stencil.push_back( 0, 1, -1, -1.0/6.0);
73 stencil.push_back( 0, 1, 1, -1.0/6.0);
74 break;
75 default:
76 assert(0 != "vmg choose discretization order 2 or 4");
77 break;
78 }
79}
80
81void DiscretizationPoissonFD::ModifyRightHandSide()
82{
83 if (order == 4) {
84
85 Stencil stencil(6.0/12.0);
86 stencil.push_back(-1, 0, 0, 1.0/12.0);
87 stencil.push_back( 1, 0, 0, 1.0/12.0);
88 stencil.push_back( 0, -1, 0, 1.0/12.0);
89 stencil.push_back( 0, 1, 0, 1.0/12.0);
90 stencil.push_back( 0, 0, -1, 1.0/12.0);
91 stencil.push_back( 0, 0, 1, 1.0/12.0);
92
93 stencil.Apply(MG::GetRhsMaxLevel());
94
95 }
96}
Note: See TracBrowser for help on using the repository browser.