Fixes for contexts.

Author: chrispbradley

src/python/fsi_tube_flow.py

@@ -211,8 +211,12 @@
 #  Initialise OpenCMISS
 #================================================================================================================================
 
+worldRegion = iron.Region()
+iron.Context.WorldRegionGet(worldRegion)
+
 # Get the computational nodes info
 computationEnvironment = iron.ComputationEnvironment()
+iron.Context.ComputationEnvironmentGet(computationEnvironment)
 numberOfComputationalNodes = computationEnvironment.NumberOfWorldNodesGet()
 computationalNodeNumber = computationEnvironment.WorldNodeNumberGet()
 
@@ -228,17 +232,17 @@
 
 # Create a RC coordinate system for the fluid region
 fluidCoordinateSystem = iron.CoordinateSystem()
-fluidCoordinateSystem.CreateStart(fluidCoordinateSystemUserNumber)
+fluidCoordinateSystem.CreateStart(fluidCoordinateSystemUserNumber,iron.Context)
 fluidCoordinateSystem.DimensionSet(3)
 fluidCoordinateSystem.CreateFinish()
 # Create a RC coordinate system for the solid region
 solidCoordinateSystem = iron.CoordinateSystem()
-solidCoordinateSystem.CreateStart(solidCoordinateSystemUserNumber)
+solidCoordinateSystem.CreateStart(solidCoordinateSystemUserNumber,iron.Context)
 solidCoordinateSystem.DimensionSet(3)
 solidCoordinateSystem.CreateFinish()
 # Create a RC coordinate system for the interface region
 interfaceCoordinateSystem = iron.CoordinateSystem()
-interfaceCoordinateSystem.CreateStart(interfaceCoordinateSystemUserNumber)
+interfaceCoordinateSystem.CreateStart(interfaceCoordinateSystemUserNumber,iron.Context)
 interfaceCoordinateSystem.DimensionSet(3)
 interfaceCoordinateSystem.CreateFinish()
 
@@ -254,13 +258,13 @@
 
 # Create a fluid region
 fluidRegion = iron.Region()
-fluidRegion.CreateStart(fluidRegionUserNumber,iron.WorldRegion)
+fluidRegion.CreateStart(fluidRegionUserNumber,worldRegion)
 fluidRegion.label = 'FluidRegion'
 fluidRegion.coordinateSystem = fluidCoordinateSystem
 fluidRegion.CreateFinish()
 # Create a solid region
 solidRegion = iron.Region()
-solidRegion.CreateStart(solidRegionUserNumber,iron.WorldRegion)
+solidRegion.CreateStart(solidRegionUserNumber,worldRegion)
 solidRegion.label = 'SolidRegion'
 solidRegion.coordinateSystem = solidCoordinateSystem
 solidRegion.CreateFinish()
@@ -279,7 +283,7 @@
 numberOfGaussXi = uInterpolation+1
 
 uBasis = iron.Basis()
-uBasis.CreateStart(uBasisUserNumber)
+uBasis.CreateStart(uBasisUserNumber,iron.Context)
 uBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP
 uBasis.numberOfXi = 3
 if (uInterpolation == LINEAR):
@@ -293,15 +297,15 @@
 uBasis.CreateFinish()
 
 pBasis = iron.Basis()
-pBasis.CreateStart(pBasisUserNumber)
+pBasis.CreateStart(pBasisUserNumber,iron.Context)
 pBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP
 pBasis.numberOfXi = 3
 pBasis.interpolationXi = [iron.BasisInterpolationSpecifications.LINEAR_LAGRANGE]*3
 pBasis.quadratureNumberOfGaussXi = [numberOfGaussXi]*3
 pBasis.CreateFinish()
 
 interfaceBasis = iron.Basis()
-interfaceBasis.CreateStart(interfaceBasisUserNumber)
+interfaceBasis.CreateStart(interfaceBasisUserNumber,iron.Context)
 interfaceBasis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP
 interfaceBasis.numberOfXi = 2
 if (uInterpolation == LINEAR):
@@ -800,7 +804,7 @@
 
 # Create an interface between the two meshes
 interface = iron.Interface()
-interface.CreateStart(interfaceUserNumber,iron.WorldRegion)
+interface.CreateStart(interfaceUserNumber,worldRegion)
 interface.LabelSet('Interface')
 # Add in the two meshes
 solidMeshIndex = interface.MeshAdd(solidMesh)
@@ -1667,7 +1671,7 @@
     fsiProblemSpecification = [iron.ProblemClasses.MULTI_PHYSICS,
                                iron.ProblemTypes.FINITE_ELASTICITY_NAVIER_STOKES,
                                iron.ProblemSubtypes.FINITE_ELASTICITY_NAVIER_STOKES_ALE]
-fsiProblem.CreateStart(fsiProblemUserNumber,fsiProblemSpecification)
+fsiProblem.CreateStart(fsiProblemUserNumber,iron.Context,fsiProblemSpecification)
 fsiProblem.CreateFinish()
 
 if (progressDiagnostics):