source: ThirdParty/vmg/src/level/stencils.hpp@ 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: 3.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 stencils.hpp
21 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
22 * @date Mon Apr 18 13:06:11 2011
23 *
24 * @brief Some examples for stencils.
25 *
26 */
27
28#ifndef STENCILS_HPP_
29#define STENCILS_HPP_
30
31#include "base/stencil.hpp"
32
33namespace VMG
34{
35
36class Injection : public Stencil
37{
38public:
39 Injection() :
40 Stencil(1.0)
41 {}
42};
43
44class InterpolationTrilinear : public Stencil
45{
46public:
47 InterpolationTrilinear() :
48 Stencil(1.0)
49 {
50 this->push_back( 1, 0, 0, 0.5);
51 this->push_back(-1, 0, 0, 0.5);
52 this->push_back( 0, 1, 0, 0.5);
53 this->push_back( 0, -1, 0, 0.5);
54 this->push_back( 0, 0, 1, 0.5);
55 this->push_back( 0, 0, -1, 0.5);
56
57 this->push_back( 1, 1, 0, 0.25);
58 this->push_back( 1, -1, 0, 0.25);
59 this->push_back(-1, 1, 0, 0.25);
60 this->push_back(-1, -1, 0, 0.25);
61 this->push_back( 0, 1, 1, 0.25);
62 this->push_back( 0, 1, -1, 0.25);
63 this->push_back( 0, -1, 1, 0.25);
64 this->push_back( 0, -1, -1, 0.25);
65 this->push_back( 1, 0, 1, 0.25);
66 this->push_back( 1, 0, -1, 0.25);
67 this->push_back(-1, 0, 1, 0.25);
68 this->push_back(-1, 0, -1, 0.25);
69
70 this->push_back( 1, 1, 1, 0.125);
71 this->push_back( 1, 1, -1, 0.125);
72 this->push_back( 1, -1, 1, 0.125);
73 this->push_back(-1, 1, 1, 0.125);
74 this->push_back( 1, -1, -1, 0.125);
75 this->push_back(-1, 1, -1, 0.125);
76 this->push_back(-1, -1, 1, 0.125);
77 this->push_back(-1, -1, -1, 0.125);
78 }
79};
80
81class RestrictionFullWeight : public Stencil
82{
83public:
84 RestrictionFullWeight() :
85 Stencil(0.125)
86 {
87 this->push_back( 1, 0, 0, 0.0625);
88 this->push_back(-1, 0, 0, 0.0625);
89 this->push_back( 0, 1, 0, 0.0625);
90 this->push_back( 0, -1, 0, 0.0625);
91 this->push_back( 0, 0, 1, 0.0625);
92 this->push_back( 0, 0, -1, 0.0625);
93
94 this->push_back( 1, 1, 0, 0.03125);
95 this->push_back( 1, -1, 0, 0.03125);
96 this->push_back(-1, 1, 0, 0.03125);
97 this->push_back(-1, -1, 0, 0.03125);
98 this->push_back( 0, 1, 1, 0.03125);
99 this->push_back( 0, 1, -1, 0.03125);
100 this->push_back( 0, -1, 1, 0.03125);
101 this->push_back( 0, -1, -1, 0.03125);
102 this->push_back( 1, 0, 1, 0.03125);
103 this->push_back( 1, 0, -1, 0.03125);
104 this->push_back(-1, 0, 1, 0.03125);
105 this->push_back(-1, 0, -1, 0.03125);
106
107 this->push_back( 1, 1, 1, 0.015625);
108 this->push_back( 1, 1, -1, 0.015625);
109 this->push_back( 1, -1, 1, 0.015625);
110 this->push_back(-1, 1, 1, 0.015625);
111 this->push_back( 1, -1, -1, 0.015625);
112 this->push_back(-1, 1, -1, 0.015625);
113 this->push_back(-1, -1, 1, 0.015625);
114 this->push_back(-1, -1, -1, 0.015625);
115 }
116};
117
118namespace Stencils
119{
120 static const VMG::Stencil Injection = VMG::Injection();
121 static const VMG::Stencil InterpolationTrilinear = VMG::InterpolationTrilinear();
122 static const VMG::Stencil RestrictionFullWeight = VMG::RestrictionFullWeight();
123}
124
125}
126
127#endif /* STENCILS_HPP_ */
Note: See TracBrowser for help on using the repository browser.