|
105 | 105 | else: |
106 | 106 | cellml_file = "input/fixedlidvelocity.cellml" |
107 | 107 |
|
| 108 | +worldRegion = iron.Region() |
| 109 | +iron.Context.WorldRegionGet(worldRegion) |
| 110 | + |
108 | 111 | # Diagnostics |
109 | 112 | #iron.DiagnosticsSetOn(iron.DiagnosticTypes.ALL,[1,2,3,4,5],"Diagnostics",[""]) |
110 | 113 | #iron.ErrorHandlingModeSet(iron.ErrorHandlingModes.TRAP_ERROR) |
111 | 114 | iron.OutputSetOn("Testing") |
112 | 115 |
|
113 | 116 | # Get the computational nodes info |
114 | 117 | computationEnvironment = iron.ComputationEnvironment() |
| 118 | +iron.Context.ComputationEnvironmentGet(computationEnvironment) |
115 | 119 | numberOfComputationalNodes = computationEnvironment.NumberOfWorldNodesGet() |
116 | 120 | computationalNodeNumber = computationEnvironment.WorldNodeNumberGet() |
117 | 121 |
|
|
173 | 177 |
|
174 | 178 | # Create a RC coordinate system for the fluid region |
175 | 179 | fluidCoordinateSystem = iron.CoordinateSystem() |
176 | | -fluidCoordinateSystem.CreateStart(fluidCoordinateSystemUserNumber) |
| 180 | +fluidCoordinateSystem.CreateStart(fluidCoordinateSystemUserNumber,iron.Context) |
177 | 181 | fluidCoordinateSystem.DimensionSet(2) |
178 | 182 | fluidCoordinateSystem.CreateFinish() |
179 | 183 | if (progressDiagnostics): |
|
188 | 192 |
|
189 | 193 | # Create a fluid region |
190 | 194 | fluidRegion = iron.Region() |
191 | | -fluidRegion.CreateStart(fluidRegionUserNumber,iron.WorldRegion) |
| 195 | +fluidRegion.CreateStart(fluidRegionUserNumber,worldRegion) |
192 | 196 | fluidRegion.label = 'FluidRegion' |
193 | 197 | fluidRegion.coordinateSystem = fluidCoordinateSystem |
194 | 198 | fluidRegion.CreateFinish() |
|
204 | 208 | print('Basis functions ...') |
205 | 209 |
|
206 | 210 | linearBasis = iron.Basis() |
207 | | -linearBasis.CreateStart(linearBasisUserNumber) |
| 211 | +linearBasis.CreateStart(linearBasisUserNumber,iron.Context) |
208 | 212 | linearBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP |
209 | 213 | linearBasis.numberOfXi = 2 |
210 | 214 | linearBasis.interpolationXi = [iron.BasisInterpolationSpecifications.LINEAR_LAGRANGE]*2 |
|
213 | 217 | if (useHermite): |
214 | 218 | numberOfNodesXi = 2 |
215 | 219 | hermiteBasis = iron.Basis() |
216 | | - hermiteBasis.CreateStart(hermiteBasisUserNumber) |
| 220 | + hermiteBasis.CreateStart(hermiteBasisUserNumber,iron.Context) |
217 | 221 | hermiteBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP |
218 | 222 | hermiteBasis.numberOfXi = 2 |
219 | 223 | hermiteBasis.interpolationXi = [iron.BasisInterpolationSpecifications.CUBIC_HERMITE]*2 |
|
222 | 226 | else: |
223 | 227 | numberOfNodesXi = 3 |
224 | 228 | quadraticBasis = iron.Basis() |
225 | | - quadraticBasis.CreateStart(quadraticBasisUserNumber) |
| 229 | + quadraticBasis.CreateStart(quadraticBasisUserNumber,iron.Context) |
226 | 230 | quadraticBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP |
227 | 231 | quadraticBasis.numberOfXi = 2 |
228 | 232 | quadraticBasis.interpolationXi = [iron.BasisInterpolationSpecifications.QUADRATIC_LAGRANGE]*2 |
|
604 | 608 | fluidProblemSpecification = [iron.ProblemClasses.FLUID_MECHANICS, |
605 | 609 | iron.ProblemTypes.NAVIER_STOKES_EQUATION, |
606 | 610 | iron.ProblemSubtypes.TRANSIENT_NAVIER_STOKES] |
607 | | -fluidProblem.CreateStart(fluidProblemUserNumber,fluidProblemSpecification) |
| 611 | +fluidProblem.CreateStart(fluidProblemUserNumber,iron.Context,fluidProblemSpecification) |
608 | 612 | fluidProblem.CreateFinish() |
609 | 613 |
|
610 | 614 | if (progressDiagnostics): |
|
0 commit comments