@@ -99,18 +99,18 @@ void Biharmonic::JR::initialize()
99
99
if (_biharmonic ._verbose )
100
100
libMesh ::out << ">>> Initializing Biharmonic::JR\n" ;
101
101
102
- Parameters parameters ;
103
- parameters .set < Point > ("center" ) = _biharmonic ._initialCenter ;
104
- parameters .set < Real > ("width" ) = _biharmonic ._initialWidth ;
102
+ Parameters params ;
103
+ params .set < Point > ("center" ) = _biharmonic ._initialCenter ;
104
+ params .set < Real > ("width" ) = _biharmonic ._initialWidth ;
105
105
106
106
if (_biharmonic ._initialState == Biharmonic ::BALL )
107
- project_solution (Biharmonic ::JR ::InitialDensityBall , Biharmonic ::JR ::InitialGradientZero , parameters );
107
+ project_solution (Biharmonic ::JR ::InitialDensityBall , Biharmonic ::JR ::InitialGradientZero , params );
108
108
109
109
if (_biharmonic ._initialState == Biharmonic ::ROD )
110
- project_solution (Biharmonic ::JR ::InitialDensityRod , Biharmonic ::JR ::InitialGradientZero , parameters );
110
+ project_solution (Biharmonic ::JR ::InitialDensityRod , Biharmonic ::JR ::InitialGradientZero , params );
111
111
112
112
if (_biharmonic ._initialState == Biharmonic ::STRIP )
113
- project_solution (Biharmonic ::JR ::InitialDensityStrip , Biharmonic ::JR ::InitialGradientZero , parameters );
113
+ project_solution (Biharmonic ::JR ::InitialDensityStrip , Biharmonic ::JR ::InitialGradientZero , params );
114
114
115
115
// both states are equal
116
116
* (old_local_solution ) = * (current_local_solution );
@@ -125,13 +125,13 @@ void Biharmonic::JR::initialize()
125
125
126
126
127
127
Number Biharmonic ::JR ::InitialDensityBall (const Point & p ,
128
- const Parameters & parameters ,
128
+ const Parameters & params ,
129
129
const std ::string & ,
130
130
const std ::string & )
131
131
{
132
132
// Initialize with a ball in the middle, which is a segment in 1D, a disk in 2D and a ball in 3D.
133
- Point center = parameters .get < Point > ("center" );
134
- Real width = parameters .get < Real > ("width" );
133
+ Point center = params .get < Point > ("center" );
134
+ Real width = params .get < Real > ("width" );
135
135
Point pc = p - center ;
136
136
Real r = pc .norm ();
137
137
return (r < width ) ? 1.0 : -0.5 ;
@@ -141,13 +141,13 @@ Number Biharmonic::JR::InitialDensityBall(const Point & p,
141
141
142
142
143
143
Number Biharmonic ::JR ::InitialDensityRod (const Point & p ,
144
- const Parameters & parameters ,
144
+ const Parameters & params ,
145
145
const std ::string & ,
146
146
const std ::string & )
147
147
{
148
148
// Initialize with a rod in the middle so that we have a z-homogeneous system to model the 2D disk.
149
- Point center = parameters .get < Point > ("center" );
150
- Real width = parameters .get < Real > ("width" );
149
+ Point center = params .get < Point > ("center" );
150
+ Real width = params .get < Real > ("width" );
151
151
Point pc = p - center ;
152
152
#if LIBMESH_DIM > 2
153
153
pc (2 ) = 0 ;
@@ -161,13 +161,13 @@ Number Biharmonic::JR::InitialDensityRod(const Point & p,
161
161
162
162
163
163
Number Biharmonic ::JR ::InitialDensityStrip (const Point & p ,
164
- const Parameters & parameters ,
164
+ const Parameters & params ,
165
165
const std ::string & ,
166
166
const std ::string & )
167
167
{
168
168
// Initialize with a wide strip in the middle so that we have a yz-homogeneous system to model the 1D.
169
- Point center = parameters .get < Point > ("center" );
170
- Real width = parameters .get < Real > ("width" );
169
+ Point center = params .get < Point > ("center" );
170
+ Real width = params .get < Real > ("width" );
171
171
Real r = sqrt ((p (0 )- center (0 ))* (p (0 )- center (0 )));
172
172
return (r < width ) ? 1.0 : -0.5 ;
173
173
}
0 commit comments