Nektar++
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
•
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
solvers
ADRSolver
EquationSystems
Laplace.cpp
Go to the documentation of this file.
1
///////////////////////////////////////////////////////////////////////////////
2
//
3
// File Laplace.cpp
4
//
5
// For more information, please see: http://www.nektar.info
6
//
7
// The MIT License
8
//
9
// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10
// Department of Aeronautics, Imperial College London (UK), and Scientific
11
// Computing and Imaging Institute, University of Utah (USA).
12
//
13
// License for the specific language governing rights and limitations under
14
// Permission is hereby granted, free of charge, to any person obtaining a
15
// copy of this software and associated documentation files (the "Software"),
16
// to deal in the Software without restriction, including without limitation
17
// the rights to use, copy, modify, merge, publish, distribute, sublicense,
18
// and/or sell copies of the Software, and to permit persons to whom the
19
// Software is furnished to do so, subject to the following conditions:
20
//
21
// The above copyright notice and this permission notice shall be included
22
// in all copies or substantial portions of the Software.
23
//
24
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
25
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
26
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
27
// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
28
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
29
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
30
// DEALINGS IN THE SOFTWARE.
31
//
32
// Description: Laplace solve routines
33
//
34
///////////////////////////////////////////////////////////////////////////////
35
36
#include <
ADRSolver/EquationSystems/Laplace.h
>
37
38
namespace
Nektar
39
{
40
string
Laplace::className
=
GetEquationSystemFactory
().
41
RegisterCreatorFunction(
"Laplace"
,
Laplace::create
);
42
43
Laplace::Laplace
(
44
const
LibUtilities::SessionReaderSharedPtr
& pSession)
45
:
EquationSystem
(pSession),
46
m_factors()
47
{
48
m_factors
[
StdRegions::eFactorLambda
] = 0.0;
49
m_factors
[
StdRegions::eFactorTau
] = 1.0;
50
}
51
52
void
Laplace::v_InitObject
()
53
{
54
EquationSystem::v_InitObject
();
55
}
56
57
Laplace::~Laplace
()
58
{
59
60
}
61
62
void
Laplace::v_GenerateSummary
(
SolverUtils::SummaryList
& s)
63
{
64
EquationSystem::SessionSummary
(s);
65
SolverUtils::AddSummaryItem
(s,
"Lambda"
,
66
m_factors
[
StdRegions::eFactorLambda
]);
67
}
68
69
void
Laplace::v_DoSolve
()
70
{
71
for
(
int
i = 0; i <
m_fields
.num_elements(); ++i)
72
{
73
// Zero field so initial conditions are zero
74
Vmath::Zero
(
m_fields
[i]->
GetNcoeffs
(),
75
m_fields
[i]->UpdateCoeffs(), 1);
76
m_fields
[i]->HelmSolve(
m_fields
[i]->GetPhys(),
77
m_fields
[i]->UpdateCoeffs(),
78
NullFlagList
,
79
m_factors
);
80
m_fields
[i]->SetPhysState(
false
);
81
}
82
}
83
84
Array<OneD, bool>
Laplace::v_GetSystemSingularChecks
()
85
{
86
return
Array<OneD, bool>(
m_session
->GetVariables().size(),
true
);
87
}
88
}
Generated on Wed Mar 4 2015 16:33:15 for Nektar++ by
1.8.1.2