Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Contexts #5

Open
wants to merge 3 commits into
base: develop
Choose a base branch
from
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 18 additions & 12 deletions src/fortran/membrane_extension_2d_space.F90
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
!> Main program
PROGRAM MEMBRANEEXTENSION2DSPACE
PROGRAM MembraneExtension2DSpace

USE OpenCMISS
USE OpenCMISS_Iron
Expand Down Expand Up @@ -69,7 +69,9 @@ PROGRAM MEMBRANEEXTENSION2DSPACE

TYPE(cmfe_BasisType) :: Basis
TYPE(cmfe_BoundaryConditionsType) :: BoundaryConditions
TYPE(cmfe_CoordinateSystemType) :: CoordinateSystem, WorldCoordinateSystem
TYPE(cmfe_ComputationEnvironmentType) :: computationEnvironment
TYPE(cmfe_ContextType) :: context
TYPE(cmfe_CoordinateSystemType) :: CoordinateSystem
TYPE(cmfe_MeshType) :: Mesh
TYPE(cmfe_DecompositionType) :: Decomposition
TYPE(cmfe_EquationsType) :: Equations
Expand Down Expand Up @@ -106,22 +108,26 @@ PROGRAM MEMBRANEEXTENSION2DSPACE
#endif

!Intialise cmiss
CALL cmfe_Initialise(WorldCoordinateSystem,WorldRegion,Err)

CALL cmfe_ErrorHandlingModeSet(CMFE_ERRORS_TRAP_ERROR,Err)
CALL cmfe_Context_Initialise(context,err)
CALL cmfe_Initialise(context,err)
CALL cmfe_ErrorHandlingModeSet(CMFE_ERRORS_TRAP_ERROR,err)
CALL cmfe_Region_Initialise(worldRegion,err)
CALL cmfe_Context_WorldRegionGet(context,worldRegion,err)

WRITE(*,'(A)') "Program starting."

!Set all diganostic levels on for testing
CALL cmfe_DiagnosticsSetOn(CMFE_FROM_DIAG_TYPE,[1,2,3,4,5],"Diagnostics",["PROBLEM_FINITE_ELEMENT_CALCULATE"],Err)

!Get the number of computational nodes and this computational node number
CALL cmfe_ComputationalNumberOfNodesGet(NumberOfComputationalNodes,Err)
CALL cmfe_ComputationalNodeNumberGet(ComputationalNodeNumber,Err)
CALL cmfe_ComputationEnvironment_Initialise(computationEnvironment,err)
CALL cmfe_Context_ComputationEnvironmentGet(context,computationEnvironment,err)
CALL cmfe_ComputationEnvironment_NumberOfWorldNodesGet(computationEnvironment,numberOfComputationalNodes,err)
CALL cmfe_ComputationEnvironment_WorldNodeNumberGet(computationEnvironment,computationalNodeNumber,err)

!Create a CS - default is 3D rectangular cartesian CS with 0,0,0 as origin
CALL cmfe_CoordinateSystem_Initialise(CoordinateSystem,Err)
CALL cmfe_CoordinateSystem_CreateStart(CoordinateSystemUserNumber,CoordinateSystem,Err)
CALL cmfe_CoordinateSystem_CreateStart(CoordinateSystemUserNumber,context,CoordinateSystem,Err)
CALL cmfe_CoordinateSystem_TypeSet(CoordinateSystem,CMFE_COORDINATE_RECTANGULAR_CARTESIAN_TYPE,Err)
CALL cmfe_CoordinateSystem_DimensionSet(CoordinateSystem,NumberOfSpatialCoordinates,Err)
CALL cmfe_CoordinateSystem_OriginSet(CoordinateSystem,[0.0_CMISSRP,0.0_CMISSRP,0.0_CMISSRP],Err)
Expand All @@ -135,7 +141,7 @@ PROGRAM MEMBRANEEXTENSION2DSPACE

!Define basis function - tri-linear Lagrange
CALL cmfe_Basis_Initialise(Basis,Err)
CALL cmfe_Basis_CreateStart(BasisUserNumber,Basis,Err)
CALL cmfe_Basis_CreateStart(BasisUserNumber,context,Basis,Err)
CALL cmfe_Basis_TypeSet(Basis,CMFE_BASIS_LAGRANGE_HERMITE_TP_TYPE,Err)
CALL cmfe_Basis_NumberOfXiSet(Basis,NumberOfXiCoordinates,Err)
CALL cmfe_Basis_InterpolationXiSet(Basis,[CMFE_BASIS_LINEAR_LAGRANGE_INTERPOLATION, &
Expand Down Expand Up @@ -276,7 +282,7 @@ PROGRAM MEMBRANEEXTENSION2DSPACE

!Define the problem
CALL cmfe_Problem_Initialise(Problem,Err)
CALL cmfe_Problem_CreateStart(ProblemUserNumber,[CMFE_PROBLEM_ELASTICITY_CLASS,CMFE_PROBLEM_FINITE_ELASTICITY_TYPE, &
CALL cmfe_Problem_CreateStart(ProblemUserNumber,context,[CMFE_PROBLEM_ELASTICITY_CLASS,CMFE_PROBLEM_FINITE_ELASTICITY_TYPE, &
& CMFE_PROBLEM_NO_SUBTYPE],Problem,Err)
CALL cmfe_Problem_CreateFinish(Problem,Err)

Expand Down Expand Up @@ -337,11 +343,11 @@ PROGRAM MEMBRANEEXTENSION2DSPACE
CALL cmfe_Fields_ElementsExport(Fields,"./results/MembraneExtension2DSpace","FORTRAN",Err)
CALL cmfe_Fields_Finalise(Fields,Err)

CALL cmfe_Finalise(Err)
CALL cmfe_Finalise(context,Err)

WRITE(*,'(A)') "Program successfully completed."

STOP

END PROGRAM MEMBRANEEXTENSION2DSPACE
END PROGRAM MembraneExtension2DSpace