Skip to content
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
29 changes: 18 additions & 11 deletions src/fortran/cellml_split_reaction_diffusion_equation.F90
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION
PROGRAM CellMLSplitReactionDiffusionEquation

USE OpenCMISS
USE OpenCMISS_Iron
Expand Down Expand Up @@ -58,7 +58,9 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION
!CMISS variables

TYPE(cmfe_BasisType) :: Basis
TYPE(cmfe_CoordinateSystemType) :: CoordinateSystem,WorldCoordinateSystem
TYPE(cmfe_ComputationEnvironmentType) :: computationEnvironment
TYPE(cmfe_ContextType) :: context
TYPE(cmfe_CoordinateSystemType) :: CoordinateSystem
TYPE(cmfe_DecompositionType) :: Decomposition
TYPE(cmfe_EquationsType) :: Equations
TYPE(cmfe_EquationsSetType) :: EquationsSet
Expand Down Expand Up @@ -107,11 +109,16 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION
!-----------------------------------------------------------------------------------------------------------

!Intialise OpenCMISS
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)
!Get the computational nodes information
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)

NUMBER_GLOBAL_X_ELEMENTS=10
NUMBER_OF_DOMAINS=NumberOfComputationalNodes
Expand All @@ -127,7 +134,7 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION

!Start the creation of a new RC coordinate system
CALL cmfe_CoordinateSystem_Initialise(CoordinateSystem,Err)
CALL cmfe_CoordinateSystem_CreateStart(CoordinateSystemUserNumber,CoordinateSystem,Err)
CALL cmfe_CoordinateSystem_CreateStart(CoordinateSystemUserNumber,context,CoordinateSystem,Err)
!Set the coordinate system to be 1D
CALL cmfe_CoordinateSystem_DimensionSet(CoordinateSystem,1,Err)
!Finish the creation of the coordinate system
Expand All @@ -152,7 +159,7 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION

!Start the creation of a basis (default is linear lagrange)
CALL cmfe_Basis_Initialise(Basis,Err)
CALL cmfe_Basis_CreateStart(BasisUserNumber,Basis,Err)
CALL cmfe_Basis_CreateStart(BasisUserNumber,context,Basis,Err)
!Set the basis to be a linear Lagrange basis
CALL cmfe_Basis_NumberOfXiSet(Basis,1,Err)
!Finish the creation of the basis
Expand Down Expand Up @@ -387,7 +394,7 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION

!Create the problem
CALL cmfe_Problem_Initialise(Problem,Err)
CALL cmfe_Problem_CreateStart(ProblemUserNumber,[CMFE_PROBLEM_CLASSICAL_FIELD_CLASS, &
CALL cmfe_Problem_CreateStart(ProblemUserNumber,context,[CMFE_PROBLEM_CLASSICAL_FIELD_CLASS, &
& CMFE_PROBLEM_REACTION_DIFFUSION_EQUATION_TYPE,CMFE_PROBLEM_CELLML_REAC_INTEG_REAC_DIFF_STRANG_SPLIT_SUBTYPE],Problem,Err)
!Finish the creation of a problem.
CALL cmfe_Problem_CreateFinish(Problem,Err)
Expand Down Expand Up @@ -534,10 +541,10 @@ PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION
CALL cmfe_Fields_Finalise(Fields,Err)
ENDIF

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

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

STOP

END PROGRAM CELLML_SPLIT_REACTION_DIFFUSION_EQUATION
END PROGRAM CellMLSplitReactionDiffusionEquation