1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
// This file is part of ff3d - http://www.freefem.org/ff3d
// Copyright (C) 2001, 2002, 2003 Stphane Del Pino
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2, or (at your option)
// any later version.
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software Foundation,
// Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
// $Id: FunctionExpressionConvection.cpp,v 1.5 2007/04/22 21:46:10 delpinux Exp $
#include <FunctionExpressionConvection.hpp>
#include <MeshOfHexahedra.hpp>
#include <Information.hpp>
#include <Convection.hpp>
std::ostream& FunctionExpressionConvection::
put(std::ostream& os) const
{
if (__scalarFunction != 0) {
os << *__scalarFunction;
} else {
os << "convect(" << *__field << ','
<< *__timeStep << ',' << *__convectedFunction << ')';
}
return os;
}
void FunctionExpressionConvection::
execute()
{
__convectedFunction->execute();
__field->execute();
__timeStep->execute();
if (Information::instance().usesMesh()) {
const Mesh& mesh = *Information::instance().getMesh();
if (__field->numberOfComponents() != 3) {
throw ErrorHandler(__FILE__,__LINE__,
"convection needs a 3 component field:\n"
+stringify(*__field)+" has "+stringify(__field->numberOfComponents())+" components",
ErrorHandler::normal);
}
switch (mesh.type()) {
case Mesh::cartesianHexahedraMesh: {
__scalarFunction
= new Convection<Structured3DMesh>(*__convectedFunction->function(),
*__field->field(),
__timeStep->realValue(),
static_cast<const Structured3DMesh&>(mesh));
break;
}
case Mesh::tetrahedraMesh: {
__scalarFunction
= new Convection<MeshOfTetrahedra>(*__convectedFunction->function(),
*__field->field(),
__timeStep->realValue(),
static_cast<const MeshOfTetrahedra&>(mesh));
break;
}
case Mesh::hexahedraMesh: {
__scalarFunction
= new Convection<MeshOfHexahedra>(*__convectedFunction->function(),
*__field->field(),
__timeStep->realValue(),
static_cast<const MeshOfHexahedra&>(mesh));
break;
}
default: {
throw ErrorHandler(__FILE__,__LINE__,
"convection is not implemented for this kind of mesh",
ErrorHandler::unexpected);
}
}
__isToEvaluate = false;
} else {
__isToEvaluate = true;
}
}
FunctionExpressionConvection::
FunctionExpressionConvection(ReferenceCounting<FieldExpression> field,
ReferenceCounting<RealExpression> dt,
ReferenceCounting<FunctionExpression> phi)
: FunctionExpression(FunctionExpression::convection),
__isToEvaluate(true),
__field(field),
__timeStep(dt),
__convectedFunction(phi)
{
;
}
FunctionExpressionConvection::
FunctionExpressionConvection(const FunctionExpressionConvection& f)
: FunctionExpression(f),
__isToEvaluate(f.__isToEvaluate),
__field(f.__field),
__timeStep(f.__timeStep),
__convectedFunction(f.__convectedFunction)
{
;
}
FunctionExpressionConvection::~FunctionExpressionConvection()
{
;
}
|
