types.f90

Go to the documentation of this file.

!#### Index: ne
!###  Description:
!###    Index label for a element.
!#### Index: ng
!###  Description:
!###    Index label for a gauss point.
!#### Index: ni
!###  Description:
!###    Index label for a xi direction.
!#### Index: nk
!###  Description:
!###    Index label for a derivative with respect to the global directions.
!#### Index: nn
!###  Description:
!###    Index for a local node within an element.
!#### Index: np
!###  Description:
!###    Index for a node.
!#### Index: ns
!###  Description:
!###    Index for a element parameter within an element.
!#### Index: nu
!###  Description:
!###    Index for a partial derivative.

MODULE types

  USE cmisspetsctypes, ONLY : petsciscoloringtype,petscksptype,petscmattype,petscmatcoloringtype,petscmatfdcoloringtype, &
    & petscpctype,petscsnestype,petscsneslinesearchtype,petscvectype
  USE constants
  USE kinds
  USE iso_c_binding
  USE iso_varying_string
  USE trees
  use linkedlist_routines

  IMPLICIT NONE

  !
  !================================================================================================================================
  !
  ! Base types
  !

  TYPE real_dp_ptr_type
    REAL(DP), POINTER :: ptr(:)
  END TYPE real_dp_ptr_type
  
  TYPE integer_intg_ptr_type
    INTEGER(INTG), POINTER :: ptr(:)
  END TYPE integer_intg_ptr_type

  TYPE integer_cint_alloc_type
    INTEGER(C_INT), ALLOCATABLE :: array(:)
  END TYPE integer_cint_alloc_type
  
  !
  !================================================================================================================================
  !
  ! List types
  !

  TYPE list_ptr_type
    TYPE(list_type), POINTER :: ptr
  END TYPE list_ptr_type

  TYPE list_type
    LOGICAL :: mutable
    LOGICAL :: list_finished
    INTEGER(INTG) :: number_in_list
    INTEGER(INTG) :: data_dimension
    INTEGER(INTG) :: initial_size
    INTEGER(INTG) :: size
    INTEGER(INTG) :: data_type
    INTEGER(INTG) :: key_dimension
    INTEGER(INTG) :: sort_order
    INTEGER(INTG) :: sort_method
    INTEGER(INTG), ALLOCATABLE :: list_intg(:)
    INTEGER(INTG), ALLOCATABLE :: list_intg2(:,:)
    REAL(SP), ALLOCATABLE :: list_sp(:)
    REAL(SP), ALLOCATABLE :: list_sp2(:,:)
    REAL(DP), ALLOCATABLE :: list_dp(:)
    REAL(DP), ALLOCATABLE :: list_dp2(:,:)
    INTEGER(C_INT), ALLOCATABLE :: list_c_int(:)
    INTEGER(C_INT), ALLOCATABLE :: list_c_int2(:,:)
  END TYPE list_type
    
  !
  !================================================================================================================================
  !
  ! Quadrature types
  !

  TYPE quadrature_scheme_type
    INTEGER(INTG) :: global_number
    TYPE(quadrature_type), POINTER :: quadrature
    INTEGER(INTG) :: number_of_gauss
    REAL(DP), ALLOCATABLE :: gauss_positions(:,:)
    REAL(DP), ALLOCATABLE :: gauss_weights(:)
    REAL(DP), ALLOCATABLE :: gauss_basis_fns(:,:,:)
    !Quadrature information at faces
    INTEGER(INTG), ALLOCATABLE :: number_of_face_gauss(:)
    REAL(DP), ALLOCATABLE :: face_gauss_basis_fns(:,:,:,:)
    REAL(DP), ALLOCATABLE :: face_gauss_positions(:,:,:)
    REAL(DP), ALLOCATABLE :: face_gauss_weights(:,:)
  END TYPE quadrature_scheme_type

  TYPE quadrature_scheme_ptr_type
    TYPE(quadrature_scheme_type), POINTER :: ptr
  END TYPE quadrature_scheme_ptr_type

  TYPE quadrature_type
    INTEGER(INTG) :: TYPE
    TYPE(basis_type), POINTER :: basis
    INTEGER(INTG), ALLOCATABLE :: number_of_gauss_xi(:)
    INTEGER(INTG) :: gauss_order
    TYPE(quadrature_scheme_ptr_type), ALLOCATABLE :: quadrature_scheme_map(:)
    INTEGER(INTG) :: number_of_schemes
    TYPE(quadrature_scheme_ptr_type), POINTER :: schemes(:)
    LOGICAL :: evaluate_face_gauss=.false. !! \todo should this be here??
  END TYPE quadrature_type

  !
  !================================================================================================================================
  !
  ! Basis types
  !

  TYPE basis_ptr_type
    TYPE(basis_type), POINTER :: ptr
  END TYPE basis_ptr_type

  TYPE basis_type
    !\todo Add in different sub types for the different types of bases???
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: family_number
    LOGICAL :: basis_finished
    LOGICAL :: hermite
    INTEGER(INTG) :: TYPE
    INTEGER(INTG) :: number_of_xi
    INTEGER(INTG) :: number_of_xi_coordinates
    INTEGER(INTG), ALLOCATABLE :: interpolation_xi(:)
    INTEGER(INTG), ALLOCATABLE :: interpolation_type(:)
    INTEGER(INTG), ALLOCATABLE :: interpolation_order(:)
    !Degenerate information
    LOGICAL :: degenerate
    INTEGER(INTG), ALLOCATABLE :: collapsed_xi(:)
    INTEGER(INTG) :: number_of_collapsed_xi
    LOGICAL, ALLOCATABLE :: node_at_collapse(:)
    !Quadrature
    TYPE(quadrature_type) :: quadrature
    INTEGER(INTG) :: number_of_partial_derivatives
    INTEGER(INTG) :: number_of_nodes
    !\todo
    INTEGER(INTG), ALLOCATABLE :: number_of_nodes_xic(:)
    INTEGER(INTG) :: number_of_element_parameters
    INTEGER(INTG) :: maximum_number_of_derivatives
    INTEGER(INTG), ALLOCATABLE :: number_of_derivatives(:)
    INTEGER(INTG), ALLOCATABLE :: node_position_index(:,:)
    INTEGER(INTG), ALLOCATABLE :: node_position_index_inv(:,:,:,:)
    INTEGER(INTG), ALLOCATABLE :: derivative_order_index(:,:,:)
    INTEGER(INTG), ALLOCATABLE :: derivative_order_index_inv(:,:,:,:)
    INTEGER(INTG), ALLOCATABLE :: partial_derivative_index(:,:)
    INTEGER(INTG), ALLOCATABLE :: element_parameter_index(:,:)
    INTEGER(INTG), ALLOCATABLE :: element_parameter_index_inv(:,:)
    !Line information
    INTEGER(INTG) :: number_of_local_lines
    INTEGER(INTG), ALLOCATABLE :: local_line_xi_direction(:)
    INTEGER(INTG), ALLOCATABLE :: number_of_nodes_in_local_line(:)
    INTEGER(INTG), ALLOCATABLE :: node_numbers_in_local_line(:,:)
    INTEGER(INTG), ALLOCATABLE :: derivative_numbers_in_local_line(:,:)
    INTEGER(INTG), ALLOCATABLE :: element_parameters_in_local_line(:,:)
    !Face information
    INTEGER(INTG) :: number_of_local_faces
    INTEGER(INTG), ALLOCATABLE :: local_face_xi_direction(:)
    INTEGER(INTG), ALLOCATABLE :: number_of_nodes_in_local_face(:)
    INTEGER(INTG), ALLOCATABLE :: node_numbers_in_local_face(:,:)
    INTEGER(INTG), ALLOCATABLE :: derivative_numbers_in_local_face(:,:,:)
    INTEGER(INTG), ALLOCATABLE :: element_parameters_in_local_face(:,:)
    !\todo What is the difference between LOCAL_XI_NORMAL and LOCAL_FACE_XI_DIRECTION ? They're the same
    INTEGER(INTG), ALLOCATABLE :: local_xi_normal(:)
    !Sub-basis information
    TYPE(basis_ptr_type), POINTER :: line_bases(:)
    TYPE(basis_ptr_type), POINTER :: face_bases(:)
    INTEGER(INTG) :: number_of_sub_bases
    TYPE(basis_ptr_type), POINTER :: sub_bases(:)
    TYPE(basis_type), POINTER :: parent_basis
  END TYPE basis_type

  TYPE basis_functions_type
   INTEGER(INTG) :: number_basis_functions
    TYPE(basis_ptr_type), POINTER :: bases(:)
  END TYPE basis_functions_type
  
  !
  !================================================================================================================================
  !
  ! Coordinate system types
  !

  TYPE, bind(c) :: coordinate_system_type
    INTEGER(INTG) :: user_number
    LOGICAL :: coordinate_system_finished
    INTEGER(INTG) :: TYPE
    INTEGER(INTG) :: radial_interpolation_type
    INTEGER(INTG) :: number_of_dimensions
    REAL(DP) :: focus
    REAL(DP) :: origin(3)
    REAL(DP) :: orientation(3,3)
  END TYPE coordinate_system_type

  TYPE coordinate_system_ptr_type
    TYPE(coordinate_system_type), POINTER :: ptr
  END TYPE coordinate_system_ptr_type

  !
  !================================================================================================================================
  !
  ! Data projection types
  !
  
  TYPE data_projection_result_type
    INTEGER(INTG) :: user_number
    REAL(DP) :: distance
    INTEGER(INTG) :: element_number
    INTEGER(INTG) :: element_face_number
    INTEGER(INTG) :: element_line_number
    INTEGER(INTG) :: exit_tag
    REAL(DP), ALLOCATABLE :: xi(:)
    REAL(DP), ALLOCATABLE :: projectionvector(:)
  END TYPE data_projection_result_type

  TYPE data_projection_type
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(varying_string) :: label
    LOGICAL :: data_projection_finished
    TYPE(data_points_type), POINTER :: data_points
    TYPE(field_type), POINTER :: projection_field
    INTEGER(INTG) :: coordinate_system_dimensions
    REAL(DP) :: maximum_iteration_update
    INTEGER(INTG) :: maximum_number_of_iterations
    TYPE(mesh_type), POINTER :: mesh
    INTEGER(INTG) :: number_of_closest_elements
    INTEGER(INTG) :: number_of_xi
    INTEGER(INTG) :: projection_type
    REAL(DP), ALLOCATABLE :: starting_xi(:)
    REAL(DP) :: absolute_tolerance
    REAL(DP) :: relative_tolerance
    INTEGER(INTG), ALLOCATABLE :: candidateelementnumbers(:)
    INTEGER(INTG), ALLOCATABLE :: localfacelinenumbers(:)
    LOGICAL :: data_projection_projected
    TYPE(data_projection_result_type), ALLOCATABLE :: data_projection_results(:)
  END TYPE data_projection_type

  TYPE data_projection_ptr_type
    TYPE(data_projection_type), POINTER :: ptr
  END TYPE data_projection_ptr_type

  !
  !================================================================================================================================
  !
  ! Data point types
  !

  TYPE data_point_type
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(varying_string) :: label
    REAL(DP), ALLOCATABLE :: position(:) !Values of the data point specifying the spatial position in the region, has the size of region dimension the data point belongs to.
    REAL(DP), ALLOCATABLE :: weights(:)
  END TYPE data_point_type

  TYPE data_points_type
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    LOGICAL :: data_points_finished
    INTEGER(INTG) :: number_of_data_points
    TYPE(data_point_type), ALLOCATABLE :: data_points(:)
    TYPE(tree_type), POINTER :: data_points_tree
    INTEGER(INTG) :: number_of_data_projections
    TYPE(data_projection_ptr_type), ALLOCATABLE :: data_projections(:)
    TYPE(tree_type), POINTER :: data_projections_tree
  END TYPE data_points_type

  !
  !================================================================================================================================
  !
  ! Node types
  !

  TYPE node_type
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(varying_string) :: label
  END TYPE node_type

  TYPE nodes_type
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    LOGICAL :: nodes_finished
    INTEGER(INTG) :: number_of_nodes
    TYPE(node_type), ALLOCATABLE :: nodes(:)
    INTEGER(INTG), ALLOCATABLE :: coupled_nodes(:,:)
    TYPE(tree_type), POINTER :: nodes_tree
  END TYPE nodes_type

  !
  !================================================================================================================================
  !
  ! Mesh types
  !

  TYPE meshdofstype
    TYPE(meshcomponenttopologytype), POINTER :: meshcomponenttopology
    INTEGER(INTG) :: numberofdofs
  END TYPE meshdofstype

  TYPE mesh_adjacent_element_type
    INTEGER(INTG) :: number_of_adjacent_elements
    INTEGER(INTG), ALLOCATABLE :: adjacent_elements(:)
  END TYPE mesh_adjacent_element_type
  
  TYPE mesh_element_type
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(basis_type), POINTER :: basis
    INTEGER(INTG), ALLOCATABLE :: mesh_element_nodes(:)
    INTEGER(INTG), ALLOCATABLE :: global_element_nodes(:)
    INTEGER(INTG), ALLOCATABLE :: user_element_node_versions(:,:)
    INTEGER(INTG), ALLOCATABLE :: user_element_nodes(:)
    TYPE(mesh_adjacent_element_type), ALLOCATABLE :: adjacent_elements(:)
    !INTEGER(INTG), ALLOCATABLE :: NUMBER_OF_ADJACENT_ELEMENTS(:) !<NUMBER_OF_ADJACENT_ELEMENTS(-ni:ni). The number of elements adjacent to this element in the ni'th xi direction. Note that -ni gives the adjacent element before the element in the ni'th direction and +ni gives the adjacent element after the element in the ni'th direction. The ni=0 index should be 1 for the current element. Old CMISS name NXI(-ni:ni,0:nei,ne).
    !INTEGER(INTG), ALLOCATABLE :: ADJACENT_ELEMENTS(:,:) !<ADJACENT_ELEMENTS(nei,-ni:ni). The local element numbers of the elements adjacent to this element in the ni'th xi direction. Note that -ni gives the adjacent elements before the element in the ni'th direction and +ni gives the adjacent elements after the element in the ni'th direction. The ni=0 index should give the current element number. Old CMISS name NXI(-ni:ni,0:nei,ne)
    LOGICAL :: boundary_element
  END TYPE mesh_element_type

  TYPE meshelementstype
    TYPE(meshcomponenttopologytype), POINTER :: meshcomponenttopology
    INTEGER(INTG) :: number_of_elements
    LOGICAL :: elements_finished
    TYPE(mesh_element_type), POINTER :: elements(:)
    TYPE(tree_type), POINTER :: elements_tree
  END TYPE meshelementstype

  TYPE meshnodederivativetype
    INTEGER(INTG) :: numberofversions !The number of global versions at the node for the mesh.
    INTEGER(INTG), ALLOCATABLE :: userversionnumbers(:) !userVersionNumbers(versionIdx). The user version numbers for the versionIdx'th version for the node.
    INTEGER(INTG), ALLOCATABLE :: dofindex(:) !The global dof version index (nv) in the domain of the nk'th global derivative for the node.
    INTEGER(INTG) :: globalderivativeindex !The global derivative index of the nk'th global derivative for the node.
    INTEGER(INTG) :: partialderivativeindex !The partial derivative index (nu) of the nk'th global derivative for the node. Old CMISS name NUNK(nk,nj,np).
  END TYPE meshnodederivativetype

  TYPE meshnodetype
    INTEGER(INTG) :: meshnumber
    INTEGER(INTG) :: globalnumber
    INTEGER(INTG) :: usernumber
    INTEGER(INTG) :: numberofderivatives
    TYPE(meshnodederivativetype), ALLOCATABLE :: derivatives(:)
    INTEGER(INTG) :: numberofsurroundingelements
    INTEGER(INTG), POINTER :: surroundingelements(:)
    LOGICAL :: boundarynode
  END TYPE meshnodetype

  TYPE meshnodestype
    TYPE(meshcomponenttopologytype), POINTER :: meshcomponenttopology
    INTEGER(INTG) :: numberofnodes
    TYPE(meshnodetype), ALLOCATABLE :: nodes(:)
    TYPE(tree_type), POINTER :: nodestree
  END TYPE meshnodestype
  
  TYPE meshelementdatapointtype
    INTEGER(INTG) :: usernumber
    INTEGER(INTG) :: globalnumber
  END TYPE meshelementdatapointtype
  
  TYPE meshelementdatapointstype
    INTEGER(INTG) :: numberofprojecteddata
    INTEGER(INTG) :: elementnumber
    TYPE(meshelementdatapointtype), ALLOCATABLE :: dataindices(:)
  END TYPE meshelementdatapointstype
  
  TYPE meshdatapointtype 
    INTEGER(INTG) :: usernumber
    INTEGER(INTG) :: globalnumber
    INTEGER(INTG) :: elementnumber
  END TYPE meshdatapointtype 
  
  
  TYPE meshdatapointstype
    TYPE(meshcomponenttopologytype), POINTER :: meshcomponenttopology
    INTEGER(INTG) :: totalnumberofprojecteddata
    TYPE(meshdatapointtype), ALLOCATABLE :: datapoints(:)
    TYPE(meshelementdatapointstype), ALLOCATABLE :: elementdatapoint(:)
  END TYPE meshdatapointstype

  TYPE meshcomponenttopologytype
    TYPE(mesh_type), POINTER :: mesh
    INTEGER(INTG) :: meshcomponentnumber
    TYPE(meshnodestype), POINTER :: nodes
    TYPE(meshelementstype), POINTER :: elements
    TYPE(meshdofstype), POINTER :: dofs
    TYPE(meshdatapointstype), POINTER :: datapoints
  END TYPE meshcomponenttopologytype

  TYPE meshcomponenttopologyptrtype
    TYPE(meshcomponenttopologytype), POINTER :: ptr
  END TYPE meshcomponenttopologyptrtype

  TYPE embedding_xi_type
    INTEGER(INTG) :: number_of_nodes
    INTEGER(INTG), ALLOCATABLE :: node_numbers(:)
    REAL(DP), ALLOCATABLE :: xi_coords(:,:)
  END TYPE embedding_xi_type

  TYPE embedding_gausspoint_type
    INTEGER(INTG) :: element_number
    REAL(DP), ALLOCATABLE :: child_xi_coord(:)
    REAL(DP), ALLOCATABLE :: parent_xi_coord(:)
  END TYPE embedding_gausspoint_type

  TYPE mesh_embedding_type
    TYPE(mesh_type), POINTER :: child_mesh, parent_mesh
    TYPE(embedding_xi_type), ALLOCATABLE :: child_node_xi_position(:)
    TYPE(embedding_gausspoint_type), ALLOCATABLE :: gauss_point_xi_position(:,:)
  END TYPE mesh_embedding_type


  TYPE mesh_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    LOGICAL :: mesh_finished
    TYPE(meshes_type), POINTER :: meshes
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    TYPE(generated_mesh_type), POINTER :: generated_mesh
    INTEGER(INTG) :: number_of_dimensions
    INTEGER(INTG) :: number_of_components
    LOGICAL :: mesh_embedded
    TYPE(mesh_type), POINTER :: embedding_mesh
    INTEGER(INTG) :: number_of_embedded_meshes
    TYPE(mesh_ptr_type), POINTER :: embedded_meshes(:)
    INTEGER(INTG) :: number_of_elements
    TYPE(meshcomponenttopologyptrtype), POINTER :: topology(:)
    TYPE(decompositions_type), POINTER :: decompositions
    LOGICAL :: surrounding_elements_calculate
  END TYPE mesh_type

  TYPE mesh_ptr_type
    TYPE(mesh_type), POINTER :: ptr
  END TYPE mesh_ptr_type

  TYPE meshes_type
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: number_of_meshes
    TYPE(mesh_ptr_type), POINTER :: meshes(:)
  END TYPE meshes_type

  !
  !================================================================================================================================
  !
  ! Generated Mesh types
  !

  TYPE generated_mesh_regular_type
    TYPE(generated_mesh_type), POINTER :: generated_mesh
    TYPE(basis_ptr_type), ALLOCATABLE :: bases(:)
    INTEGER(INTG) :: coordinate_dimension
    INTEGER(INTG) :: mesh_dimension
    REAL(DP), ALLOCATABLE :: origin(:)
    REAL(DP), ALLOCATABLE :: maximum_extent(:)
    REAL(DP), ALLOCATABLE :: base_vectors(:,:)
    INTEGER(INTG), ALLOCATABLE :: number_of_elements_xi(:)
  END TYPE generated_mesh_regular_type

  TYPE generated_mesh_cylinder_type
    TYPE(generated_mesh_type), POINTER :: generated_mesh
    REAL(DP), ALLOCATABLE :: origin(:)
    REAL(DP), ALLOCATABLE :: cylinder_extent(:)
    INTEGER(INTG) :: mesh_dimension
    INTEGER(INTG), ALLOCATABLE :: number_of_elements_xi(:)
    TYPE(basis_ptr_type), ALLOCATABLE :: bases(:)
    LOGICAL :: append_linear_component=.false. 
 END TYPE generated_mesh_cylinder_type
 
 TYPE generated_mesh_ellipsoid_type
    TYPE(generated_mesh_type), POINTER :: generated_mesh
    REAL(DP), ALLOCATABLE :: origin(:)
    REAL(DP), ALLOCATABLE :: ellipsoid_extent(:)
    INTEGER(INTG) :: mesh_dimension
    INTEGER(INTG), ALLOCATABLE :: number_of_elements_xi(:)
    TYPE(basis_ptr_type), ALLOCATABLE :: bases(:)
    LOGICAL :: append_linear_component=.false. 
END TYPE generated_mesh_ellipsoid_type

  TYPE generated_mesh_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    TYPE(generated_meshes_type), POINTER :: generated_meshes
    LOGICAL :: generated_mesh_finished
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: generated_type
    TYPE(generated_mesh_regular_type), POINTER :: regular_mesh
    TYPE(generated_mesh_cylinder_type), POINTER :: cylinder_mesh
    TYPE(generated_mesh_ellipsoid_type), POINTER :: ellipsoid_mesh
    TYPE(mesh_type), POINTER :: mesh
  END TYPE generated_mesh_type
  
  TYPE generated_mesh_ptr_type
    TYPE(generated_mesh_type), POINTER :: ptr
  END TYPE generated_mesh_ptr_type
       
  TYPE generated_meshes_type
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: number_of_generated_meshes
    TYPE(generated_mesh_ptr_type), POINTER :: generated_meshes(:)
   END TYPE generated_meshes_type
  
  !
  !================================================================================================================================
  !
  ! Domain types
  !
  
  TYPE domain_dofs_type
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: number_of_dofs
    INTEGER(INTG) :: total_number_of_dofs
    INTEGER(INTG) :: number_of_global_dofs
    INTEGER(INTG), ALLOCATABLE :: dof_index(:,:)
  END TYPE domain_dofs_type

  TYPE domain_line_type
    INTEGER(INTG) :: number
    TYPE(basis_type), POINTER :: basis
    INTEGER(INTG), ALLOCATABLE :: nodes_in_line(:)
    INTEGER(INTG), ALLOCATABLE :: derivatives_in_line(:,:,:)
    LOGICAL :: boundary_line
    INTEGER(INTG) :: element_number
  END TYPE domain_line_type

  TYPE domain_line_ptr_type
    TYPE(domain_line_type), POINTER :: ptr
  END TYPE domain_line_ptr_type

  TYPE domain_lines_type
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: number_of_lines
    TYPE(domain_line_type), ALLOCATABLE :: lines(:)
  END TYPE domain_lines_type

  TYPE domain_face_type
    INTEGER(INTG) :: number
    INTEGER(INTG) :: xi_direction1
    INTEGER(INTG) :: xi_direction2
    TYPE(basis_type), POINTER :: basis
    INTEGER(INTG), ALLOCATABLE :: nodes_in_face(:)
    INTEGER(INTG), ALLOCATABLE :: derivatives_in_face(:,:,:)
    LOGICAL :: boundary_face
    INTEGER(INTG) :: element_number
  END TYPE domain_face_type

  TYPE domain_face_ptr_type
    TYPE(domain_face_type), POINTER :: ptr
  END TYPE domain_face_ptr_type

  TYPE domain_faces_type
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: number_of_faces
    TYPE(domain_face_type), ALLOCATABLE :: faces(:)
  END TYPE domain_faces_type

  TYPE domain_element_type
    INTEGER(INTG) :: number
    TYPE(basis_type), POINTER :: basis
    INTEGER(INTG), ALLOCATABLE :: element_nodes(:)
    INTEGER(INTG), ALLOCATABLE :: element_derivatives(:,:)
    INTEGER(INTG), ALLOCATABLE :: elementversions(:,:)
  END TYPE domain_element_type
  
  TYPE domain_elements_type
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: number_of_elements
    INTEGER(INTG) :: total_number_of_elements
    INTEGER(INTG) :: number_of_global_elements
    TYPE(domain_element_type), POINTER :: elements(:)
    INTEGER(INTG) :: maximum_number_of_element_parameters
  END TYPE domain_elements_type

  TYPE domain_node_derivative_type
    INTEGER(INTG) :: numberofversions !The number of global versions at the node for the mesh.
    INTEGER(INTG), ALLOCATABLE :: userversionnumbers(:) !The user version index of the nk'th global derivative for the node.
    INTEGER(INTG), ALLOCATABLE :: dof_index(:) !The local dof derivative version index in the domain of the nk'th global derivative for the node.
    INTEGER(INTG) :: global_derivative_index !The global derivative index of the nk'th global derivative for the node.
    INTEGER(INTG) :: partial_derivative_index !The partial derivative index (nu) of the nk'th global derivative for the node. Old CMISS name NUNK(nk,nj,np).
  END TYPE domain_node_derivative_type

  TYPE domain_node_type
    INTEGER(INTG) :: local_number
    INTEGER(INTG) :: mesh_number
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: number_of_derivatives
    TYPE(domain_node_derivative_type), ALLOCATABLE :: derivatives(:)
    INTEGER(INTG) :: number_of_surrounding_elements
    INTEGER(INTG), POINTER :: surrounding_elements(:)
    INTEGER(INTG) :: number_of_node_lines
    INTEGER(INTG), ALLOCATABLE :: node_lines(:)
    INTEGER(INTG) :: number_of_node_faces
    INTEGER(INTG), ALLOCATABLE :: node_faces(:)
    LOGICAL :: boundary_node
  END TYPE domain_node_type

  TYPE domain_nodes_type
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: number_of_nodes
    INTEGER(INTG) :: total_number_of_nodes
    INTEGER(INTG) :: number_of_global_nodes
    INTEGER(INTG) :: maximum_number_of_derivatives
    TYPE(domain_node_type), POINTER :: nodes(:)
    TYPE(tree_type), POINTER :: nodes_tree
  END TYPE domain_nodes_type

  TYPE domain_topology_type
    TYPE(domain_type), POINTER :: domain
    TYPE(domain_nodes_type), POINTER :: nodes
    TYPE(domain_dofs_type), POINTER :: dofs
    TYPE(domain_elements_type), POINTER :: elements
    TYPE(domain_faces_type), POINTER :: faces
    TYPE(domain_lines_type), POINTER :: lines
  END TYPE domain_topology_type

  !
  !================================================================================================================================
  !
  ! Distributed matrix vector types
  !

  TYPE distributed_vector_transfer_type
    TYPE(distributed_vector_cmiss_type), POINTER :: cmiss_vector
    INTEGER(INTG) :: data_type
    INTEGER(INTG) :: send_buffer_size
    INTEGER(INTG) :: receive_buffer_size
    INTEGER(INTG) :: send_tag_number
    INTEGER(INTG) :: receive_tag_number
    INTEGER(INTG) :: mpi_send_request
    INTEGER(INTG) :: mpi_receive_request
    INTEGER(INTG), ALLOCATABLE :: send_buffer_intg(:)
    REAL(DP), ALLOCATABLE :: send_buffer_dp(:)
    REAL(SP), ALLOCATABLE :: send_buffer_sp(:)
    LOGICAL, ALLOCATABLE :: send_buffer_l(:)
    INTEGER(INTG), ALLOCATABLE :: receive_buffer_intg(:)
    REAL(DP), ALLOCATABLE :: receive_buffer_dp(:)
    REAL(SP), ALLOCATABLE :: receive_buffer_sp(:)
    LOGICAL, ALLOCATABLE :: receive_buffer_l(:)
  END TYPE distributed_vector_transfer_type

  TYPE distributed_vector_cmiss_type
    TYPE(distributed_vector_type), POINTER :: distributed_vector
    INTEGER(INTG) :: base_tag_number
    INTEGER(INTG) :: n
    INTEGER(INTG) :: data_size
    INTEGER(INTG), ALLOCATABLE :: data_intg(:)
    REAL(DP), ALLOCATABLE :: data_dp(:)
    REAL(SP), ALLOCATABLE :: data_sp(:)
    LOGICAL, ALLOCATABLE :: data_l(:)
    TYPE(distributed_vector_transfer_type), ALLOCATABLE :: transfers(:)
  END TYPE distributed_vector_cmiss_type

  TYPE distributed_vector_petsc_type
    TYPE(distributed_vector_type), POINTER :: distributed_vector
    INTEGER(INTG) :: n
    INTEGER(INTG) :: global_n
    INTEGER(INTG) :: data_size
    INTEGER(INTG), ALLOCATABLE :: global_numbers(:)
    LOGICAL :: use_override_vector
    TYPE(petscvectype) :: vector
    TYPE(petscvectype) :: override_vector
  END TYPE distributed_vector_petsc_type
  
  TYPE distributed_vector_type
    LOGICAL :: vector_finished
    INTEGER(INTG) :: library_type
    INTEGER(INTG) :: ghosting_type
    TYPE(domain_mapping_type), POINTER :: domain_mapping
    INTEGER(INTG) :: data_type
    TYPE(distributed_vector_cmiss_type), POINTER :: cmiss
    TYPE(distributed_vector_petsc_type), POINTER :: petsc
  END TYPE distributed_vector_type

  TYPE distributed_matrix_cmiss_type
    TYPE(distributed_matrix_type), POINTER :: distributed_matrix
    INTEGER(INTG) :: base_tag_number
    TYPE(matrix_type), POINTER :: matrix
!   TYPE(LINKEDLIST),POINTER :: LIST(:) !< \todo Comment
  END TYPE distributed_matrix_cmiss_type

  TYPE distributed_matrix_petsc_type
    TYPE(distributed_matrix_type), POINTER :: distributed_matrix
    INTEGER(INTG) :: m
    INTEGER(INTG) :: n
    INTEGER(INTG) :: global_m
    INTEGER(INTG) :: global_n
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: number_non_zeros
    INTEGER(INTG) :: data_size
    INTEGER(INTG) :: maximum_column_indices_per_row
    INTEGER(INTG), ALLOCATABLE :: diagonal_number_non_zeros(:)
    INTEGER(INTG), ALLOCATABLE :: offdiagonal_number_non_zeros(:)
    INTEGER(INTG), ALLOCATABLE :: row_indices(:)
    INTEGER(INTG), ALLOCATABLE :: column_indices(:)
    TYPE(linkedlist),POINTER :: list(:)
    INTEGER(INTG), ALLOCATABLE :: global_row_numbers(:)
    REAL(DP), POINTER :: data_dp(:)
    LOGICAL :: use_override_matrix
    TYPE(petscmattype) :: matrix
    TYPE(petscmattype) :: override_matrix
  END TYPE distributed_matrix_petsc_type
  
  TYPE distributed_matrix_type
    LOGICAL :: matrix_finished
    INTEGER(INTG) :: library_type
    INTEGER(INTG) :: ghosting_type
    TYPE(domain_mapping_type), POINTER :: row_domain_mapping
    TYPE(domain_mapping_type), POINTER :: column_domain_mapping
    INTEGER(INTG) :: data_type
    TYPE(distributed_matrix_cmiss_type), POINTER :: cmiss
    TYPE(distributed_matrix_petsc_type), POINTER :: petsc
  END TYPE distributed_matrix_type

  !
  !================================================================================================================================
  !
  ! Matrix vector types
  !

  TYPE vector_type
    INTEGER(INTG) :: id
    LOGICAL :: vector_finished
    INTEGER(INTG) :: n
    INTEGER(INTG) :: data_type
    INTEGER(INTG) :: size
    INTEGER(INTG), ALLOCATABLE :: data_intg(:)
    REAL(SP), ALLOCATABLE :: data_sp(:)
    REAL(DP), ALLOCATABLE :: data_dp(:)
    LOGICAL, ALLOCATABLE :: data_l(:)
  END TYPE vector_type

  TYPE matrix_type
    INTEGER(INTG) :: id
    LOGICAL :: matrix_finished
    INTEGER(INTG) :: m
    INTEGER(INTG) :: n
    INTEGER(INTG) :: max_m
    INTEGER(INTG) :: max_n
    INTEGER(INTG) :: data_type
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: number_non_zeros
    INTEGER(INTG) :: size
    INTEGER(INTG) :: maximum_column_indices_per_row
    INTEGER(INTG), ALLOCATABLE :: row_indices(:)
    INTEGER(INTG), ALLOCATABLE :: column_indices(:)
    TYPE(linkedlist),POINTER :: list(:) !\todo Comment
    INTEGER(INTG), ALLOCATABLE :: data_intg(:)
    REAL(SP), ALLOCATABLE :: data_sp(:)
    REAL(DP), ALLOCATABLE :: data_dp(:)
    LOGICAL, ALLOCATABLE :: data_l(:)
  END TYPE matrix_type

  !
  !================================================================================================================================
  !
  ! Domain mapping types
  !
  
  TYPE domain_adjacent_domain_type
    INTEGER(INTG) :: domain_number
    INTEGER(INTG) :: number_of_send_ghosts
    INTEGER(INTG) :: number_of_receive_ghosts
    INTEGER(INTG), ALLOCATABLE :: local_ghost_send_indices(:)
    INTEGER(INTG), ALLOCATABLE :: local_ghost_receive_indices(:)
  END TYPE domain_adjacent_domain_type
  
  TYPE domain_global_mapping_type
    INTEGER(INTG) :: number_of_domains
    INTEGER(INTG), ALLOCATABLE :: local_number(:)
    INTEGER(INTG), ALLOCATABLE :: domain_number(:)
    INTEGER(INTG), ALLOCATABLE :: local_type(:)
  END TYPE domain_global_mapping_type

  TYPE domain_mapping_type
    INTEGER(INTG) :: number_of_local
    INTEGER(INTG) :: total_number_of_local
    INTEGER(INTG), ALLOCATABLE :: number_of_domain_local(:)
    INTEGER(INTG), ALLOCATABLE :: number_of_domain_ghost(:)
    INTEGER(INTG) :: number_of_global
    INTEGER(INTG) :: number_of_domains
    INTEGER(INTG) :: number_of_internal
    INTEGER(INTG) :: number_of_boundary
    INTEGER(INTG) :: number_of_ghost
    INTEGER(INTG) :: internal_start
    INTEGER(INTG) :: internal_finish
    INTEGER(INTG) :: boundary_start
    INTEGER(INTG) :: boundary_finish
    INTEGER(INTG) :: ghost_start
    INTEGER(INTG) :: ghost_finish
    INTEGER(INTG), ALLOCATABLE :: domain_list(:)
    INTEGER(INTG), ALLOCATABLE :: local_to_global_map(:)
    TYPE(domain_global_mapping_type), ALLOCATABLE :: global_to_local_map(:)
    INTEGER(INTG) :: number_of_adjacent_domains
    INTEGER(INTG), ALLOCATABLE :: adjacent_domains_ptr(:)
    INTEGER(INTG), ALLOCATABLE :: adjacent_domains_list(:)
    TYPE(domain_adjacent_domain_type), ALLOCATABLE :: adjacent_domains(:)
  END TYPE domain_mapping_type

  TYPE domain_mappings_type
    TYPE(domain_type), POINTER :: domain
    TYPE(domain_mapping_type), POINTER :: elements
    TYPE(domain_mapping_type), POINTER :: nodes
    TYPE(domain_mapping_type), POINTER :: dofs
  END TYPE domain_mappings_type
  
  TYPE domain_type
    TYPE(decomposition_type), POINTER :: decomposition
    TYPE(mesh_type), POINTER :: mesh
    INTEGER(INTG) :: mesh_component_number
    TYPE(region_type), POINTER :: region
    INTEGER(INTG) :: number_of_dimensions
    INTEGER(INTG), ALLOCATABLE :: node_domain(:)
    TYPE(domain_mappings_type), POINTER :: mappings
    TYPE(domain_topology_type), POINTER :: topology
  END TYPE domain_type

  TYPE domain_ptr_type 
    TYPE(domain_type), POINTER :: ptr
  END TYPE domain_ptr_type

  !
  !================================================================================================================================
  !
  ! Decomposition types
  !
  
  TYPE decomposition_line_type
    INTEGER(INTG) :: number
    INTEGER(INTG) :: xi_direction
    INTEGER(INTG) :: number_of_surrounding_elements
    INTEGER(INTG), ALLOCATABLE :: surrounding_elements(:)
    INTEGER(INTG), ALLOCATABLE :: element_lines(:)
    INTEGER(INTG) :: adjacent_lines(0:1)
    LOGICAL :: boundary_line
  END TYPE decomposition_line_type

  TYPE decomposition_lines_type
    TYPE(decomposition_type), POINTER :: decomposition
    INTEGER(INTG) :: number_of_lines
    TYPE(decomposition_line_type), ALLOCATABLE :: lines(:)
  END TYPE decomposition_lines_type

  TYPE decomposition_face_type
    INTEGER(INTG) :: number
    INTEGER(INTG) :: xi_direction
    INTEGER(INTG) :: number_of_surrounding_elements
    INTEGER(INTG), ALLOCATABLE :: surrounding_elements(:)
    INTEGER(INTG), ALLOCATABLE :: element_faces(:)
!    INTEGER(INTG) :: ADJACENT_FACES(0:1) !<ADJACENT_FACES(0:1). The face number of adjacent faces. ADJACENT_FACES(0) is the face number adjacent in the -xi direction. ADJACENT_FACES(1) is the face number adjacent in the +xi direction. Old CMISS name NPL(2..3,0,nl).
    LOGICAL :: boundary_face
    INTEGER(INTG) :: element_number
  END TYPE decomposition_face_type

  TYPE decomposition_faces_type
    TYPE(decomposition_type), POINTER :: decomposition
    INTEGER(INTG) :: number_of_faces
    TYPE(decomposition_face_type), ALLOCATABLE :: faces(:)
  END TYPE decomposition_faces_type

  TYPE decomposition_adjacent_element_type
    INTEGER(INTG) :: number_of_adjacent_elements
    INTEGER(INTG), ALLOCATABLE :: adjacent_elements(:)
  END TYPE decomposition_adjacent_element_type
  
  TYPE decomposition_element_type
    INTEGER(INTG) :: local_number
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(decomposition_adjacent_element_type), ALLOCATABLE :: adjacent_elements(:)
    !\todo INTEGER(INTG), ALLOCATABLE :: NUMBER_OF_ADJACENT_ELEMENTS(:) !<NUMBER_OF_ADJACENT_ELEMENTS(-ni:ni). The number of elements adjacent to this element in the ni'th xi direction. Note that -ni gives the adjacent element before the element in the ni'th direction and +ni gives the adjacent element after the element in the ni'th direction. The ni=0 index should be 1 for the current element. Old CMISS name NXI(-ni:ni,0:nei,ne).
    !\todo INTEGER(INTG), ALLOCATABLE :: ADJACENT_ELEMENTS(:,:) !<ADJACENT_ELEMENTS(nei,-ni:ni). The local element numbers of the elements adjacent to this element in the ni'th xi direction. Note that -ni gives the adjacent elements before the element in the ni'th direction and +ni gives the adjacent elements after the element in the ni'th direction. The ni=0 index should give the current element number. Old CMISS name NXI(-ni:ni,0:nei,ne).
    INTEGER(INTG), ALLOCATABLE :: element_lines(:)
    INTEGER(INTG), ALLOCATABLE :: element_faces(:)
    LOGICAL :: boundary_element
  END TYPE decomposition_element_type

  TYPE decomposition_elements_type
    TYPE(decomposition_type), POINTER :: decomposition
    INTEGER(INTG) :: number_of_elements
    INTEGER(INTG) :: total_number_of_elements
    INTEGER(INTG) :: number_of_global_elements
    TYPE(decomposition_element_type), POINTER :: elements(:)
    TYPE(tree_type), POINTER :: elements_tree
  END TYPE decomposition_elements_type
  
  TYPE decompositionelementdatapointtype
    INTEGER(INTG) :: usernumber
    INTEGER(INTG) :: globalnumber
    INTEGER(INTG) :: localnumber
  END TYPE decompositionelementdatapointtype
  
  TYPE decompositionelementdatapointstype
    INTEGER(INTG) :: numberofprojecteddata
    INTEGER(INTG) :: globalelementnumber
    TYPE(decompositionelementdatapointtype), ALLOCATABLE :: dataindices(:)
  END TYPE decompositionelementdatapointstype
  
  TYPE decompositiondatapointstype
    TYPE(decomposition_type), POINTER :: decomposition
    INTEGER(INTG) :: numberofdatapoints
    INTEGER(INTG) :: totalnumberofdatapoints
    INTEGER(INTG) :: numberofglobaldatapoints
    INTEGER(INTG), ALLOCATABLE :: numberofdomainlocal(:)
    INTEGER(INTG), ALLOCATABLE :: numberofdomainghost(:)
    INTEGER(INTG), ALLOCATABLE :: numberofelementdatapoints(:)
    TYPE(decompositionelementdatapointstype), ALLOCATABLE :: elementdatapoint(:)
    TYPE(tree_type), POINTER :: datapointstree
  END TYPE decompositiondatapointstype

  TYPE decomposition_topology_type
    TYPE(decomposition_type), POINTER :: decomposition
    TYPE(decomposition_elements_type), POINTER :: elements
    TYPE(decomposition_lines_type), POINTER :: lines
    TYPE(decomposition_faces_type), POINTER :: faces
    TYPE(decompositiondatapointstype), POINTER :: datapoints
  END TYPE decomposition_topology_type

  TYPE decomposition_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    LOGICAL :: decomposition_finished
    TYPE(decompositions_type), POINTER :: decompositions
    TYPE(mesh_type), POINTER :: mesh
    INTEGER(INTG) :: mesh_component_number
    INTEGER(INTG) :: decomposition_type
    INTEGER(INTG) :: number_of_domains
    INTEGER(INTG) :: number_of_edges_cut
    INTEGER(INTG), ALLOCATABLE :: element_domain(:)
    TYPE(decomposition_topology_type), POINTER :: topology
    TYPE(domain_ptr_type), POINTER :: domain(:)
    LOGICAL :: calculate_faces
    LOGICAL :: calculate_lines
  END TYPE decomposition_type

  TYPE decomposition_ptr_type
    TYPE(decomposition_type), POINTER :: ptr
  END TYPE decomposition_ptr_type

  TYPE decompositions_type
    TYPE(mesh_type), POINTER :: mesh
    INTEGER(INTG) :: number_of_decompositions
    TYPE(decomposition_ptr_type), POINTER :: decompositions(:)
  END TYPE decompositions_type

  !
  !================================================================================================================================
  !
  ! Field types
  !

  TYPE field_physical_point_type
    TYPE(field_interpolated_point_type), POINTER :: field_interpolated_point
    TYPE(field_interpolated_point_type), POINTER :: geometric_interpolated_point
    INTEGER(INTG) :: physical_derivative_type
    REAL(DP), ALLOCATABLE :: values(:)
  END TYPE field_physical_point_type

  TYPE field_physical_point_ptr_type
    TYPE(field_physical_point_type), POINTER :: ptr
  END TYPE field_physical_point_ptr_type
  
  TYPE field_interpolated_point_metrics_type
    TYPE(field_interpolated_point_type), POINTER :: interpolated_point
    INTEGER(INTG) :: number_of_x_dimensions
    INTEGER(INTG) :: number_of_xi_dimensions
    REAL(DP), ALLOCATABLE :: gl(:,:)
    REAL(DP), ALLOCATABLE :: gu(:,:)
    REAL(DP), ALLOCATABLE :: dx_dxi(:,:)
    REAL(DP), ALLOCATABLE :: dxi_dx(:,:)
    REAL(DP) :: jacobian
    INTEGER(INTG) :: jacobian_type
  END TYPE field_interpolated_point_metrics_type

  TYPE field_interpolated_point_metrics_ptr_type
    TYPE(field_interpolated_point_metrics_type), POINTER :: ptr
  END TYPE field_interpolated_point_metrics_ptr_type
  
  TYPE field_interpolated_point_type
    TYPE(field_interpolation_parameters_type), POINTER :: interpolation_parameters
    INTEGER(INTG) :: max_partial_derivative_index
    INTEGER(INTG) :: partial_derivative_type
    REAL(DP), ALLOCATABLE :: values(:,:)
  END TYPE field_interpolated_point_type

  TYPE field_interpolated_point_ptr_type
    TYPE(field_interpolated_point_type), POINTER :: ptr
  END TYPE field_interpolated_point_ptr_type
  
  TYPE field_interpolation_parameters_type
    TYPE(field_type), POINTER :: field
    TYPE(field_variable_type), POINTER :: field_variable
    INTEGER(INTG) :: number_of_xi
    TYPE(basis_ptr_type), ALLOCATABLE :: bases(:)
    INTEGER(INTG), ALLOCATABLE :: number_of_parameters(:)
    REAL(DP), ALLOCATABLE :: parameters(:,:)
    REAL(DP), ALLOCATABLE :: scale_factors(:,:)
  END TYPE field_interpolation_parameters_type

  TYPE field_interpolation_parameters_ptr_type
    TYPE(field_interpolation_parameters_type), POINTER :: ptr
  END TYPE field_interpolation_parameters_ptr_type
  
  TYPE field_geometric_parameters_type
    INTEGER(INTG) :: number_of_lines
    INTEGER(INTG) :: number_of_areas
    INTEGER(INTG) :: number_of_volumes
    REAL(DP), ALLOCATABLE :: lengths(:)
    REAL(DP), ALLOCATABLE :: areas(:)
    REAL(DP), ALLOCATABLE :: volumes(:)
    INTEGER(INTG) :: number_of_fields_using
    TYPE(field_ptr_type), POINTER :: fields_using(:)
  END TYPE field_geometric_parameters_type

  TYPE field_scaling_type
    INTEGER(INTG) :: mesh_component_number
    INTEGER(INTG) :: max_number_of_derivatives
    INTEGER(INTG) :: max_number_of_element_parameters
    TYPE(distributed_vector_type), POINTER :: scale_factors
  END TYPE field_scaling_type

  TYPE field_scalings_type
    INTEGER(INTG) :: scaling_type
    INTEGER(INTG) :: number_of_scaling_indices
    TYPE(field_scaling_type), ALLOCATABLE :: scalings(:)
  END TYPE field_scalings_type

  TYPE field_dof_to_param_map_type
    INTEGER(INTG) :: number_of_dofs
    INTEGER(INTG), ALLOCATABLE :: dof_type(:,:)
    INTEGER(INTG) :: number_of_constant_dofs
    INTEGER(INTG) :: number_of_element_dofs
    INTEGER(INTG) :: number_of_node_dofs
    INTEGER(INTG) :: number_of_grid_point_dofs
    INTEGER(INTG) :: number_of_gauss_point_dofs
    INTEGER(INTG) :: number_of_data_point_dofs
    INTEGER(INTG), ALLOCATABLE :: constant_dof2param_map(:)
    INTEGER(INTG), ALLOCATABLE :: element_dof2param_map(:,:)
    INTEGER(INTG), ALLOCATABLE :: node_dof2param_map(:,:)
    INTEGER(INTG), ALLOCATABLE :: grid_point_dof2param_map(:,:)
    INTEGER(INTG), ALLOCATABLE :: gauss_point_dof2param_map(:,:)
    INTEGER(INTG), ALLOCATABLE :: data_point_dof2param_map(:,:)
  END TYPE field_dof_to_param_map_type

  TYPE field_node_param_to_dof_map_derivative_type
    INTEGER(INTG) :: number_of_versions
    INTEGER(INTG), ALLOCATABLE :: versions(:)
  END TYPE field_node_param_to_dof_map_derivative_type

  TYPE field_node_param_to_dof_map_node_type
    INTEGER(INTG) :: number_of_derivatives
    TYPE(field_node_param_to_dof_map_derivative_type), ALLOCATABLE :: derivatives(:) ! The mapping from field node derivative parameter to a dof
  END TYPE field_node_param_to_dof_map_node_type

  TYPE field_node_param_to_dof_map_type
    INTEGER(INTG) :: number_of_node_parameters
    TYPE(field_node_param_to_dof_map_node_type), ALLOCATABLE :: nodes(:)  ! The mapping from field node parameter to a dof
  END TYPE field_node_param_to_dof_map_type

  TYPE field_element_param_to_dof_map_type
    INTEGER(INTG) :: number_of_element_parameters
    INTEGER(INTG), ALLOCATABLE :: elements(:)
  END TYPE field_element_param_to_dof_map_type

  TYPE field_grid_point_param_to_dof_map_type
    INTEGER(INTG) :: number_of_grid_point_parameters
    INTEGER(INTG), ALLOCATABLE :: grid_points(:)
  END TYPE field_grid_point_param_to_dof_map_type

  TYPE field_gauss_point_param_to_dof_map_type
    INTEGER(INTG) :: number_of_gauss_point_parameters
    INTEGER(INTG), ALLOCATABLE :: gauss_points(:,:)
  END TYPE field_gauss_point_param_to_dof_map_type
  
  TYPE field_data_point_param_to_dof_map_type
    INTEGER(INTG) :: number_of_data_point_parameters
    INTEGER(INTG), ALLOCATABLE :: data_points(:)
  END TYPE field_data_point_param_to_dof_map_type

  TYPE field_param_to_dof_map_type
    INTEGER(INTG) :: number_of_constant_parameters
    INTEGER(INTG) :: constant_param2dof_map
    TYPE(field_element_param_to_dof_map_type) :: element_param2dof_map
    TYPE(field_node_param_to_dof_map_type) :: node_param2dof_map
    TYPE(field_grid_point_param_to_dof_map_type) :: grid_point_param2dof_map
    TYPE(field_gauss_point_param_to_dof_map_type) :: gauss_point_param2dof_map
    TYPE(field_data_point_param_to_dof_map_type) :: data_point_param2dof_map
  END TYPE field_param_to_dof_map_type

  TYPE field_variable_component_type
    INTEGER(INTG) :: component_number
    TYPE(field_variable_type), POINTER :: field_variable
    TYPE(varying_string) :: component_label
    INTEGER(INTG) :: interpolation_type
    INTEGER(INTG) :: mesh_component_number
    INTEGER(INTG) :: scaling_index
    TYPE(domain_type), POINTER :: domain
    INTEGER(INTG) :: maxnumberelementinterpolationparameters
    INTEGER(INTG) :: maxnumbernodeinterpolationparameters
    TYPE(field_param_to_dof_map_type) :: param_to_dof_map
  END TYPE field_variable_component_type

  TYPE field_parameter_set_type
    INTEGER(INTG) :: set_index
    INTEGER(INTG) :: set_type
  !###      corresponds to.
    TYPE(distributed_vector_type), POINTER :: parameters
  END TYPE field_parameter_set_type
  
  TYPE field_parameter_set_ptr_type
    TYPE(field_parameter_set_type), POINTER :: ptr
  END TYPE field_parameter_set_ptr_type

  TYPE field_parameter_sets_type    
    TYPE(field_variable_type), POINTER :: field_variable
    INTEGER(INTG) :: number_of_parameter_sets
    TYPE(field_parameter_set_ptr_type), POINTER :: set_type(:)
    TYPE(field_parameter_set_ptr_type), POINTER :: parameter_sets(:)
  END TYPE field_parameter_sets_type
  
  TYPE field_variable_type
    INTEGER(INTG) :: variable_number
    INTEGER(INTG) :: variable_type
    TYPE(varying_string) :: variable_label
    TYPE(field_type), POINTER :: field
    TYPE(region_type), POINTER :: region
    INTEGER(INTG) :: dimension
    INTEGER(INTG) :: data_type
    INTEGER(INTG) :: dof_order_type
    INTEGER(INTG) :: maxnumberelementinterpolationparameters
    INTEGER(INTG) :: maxnumbernodeinterpolationparameters
    INTEGER(INTG) :: number_of_dofs
    INTEGER(INTG) :: total_number_of_dofs
    INTEGER(INTG) :: number_of_global_dofs
    TYPE(domain_mapping_type), POINTER :: domain_mapping
    INTEGER(INTG) :: number_of_components
    TYPE(field_variable_component_type), ALLOCATABLE :: components(:)
    TYPE(field_dof_to_param_map_type) :: dof_to_param_map
    TYPE(field_parameter_sets_type) :: parameter_sets
  END TYPE field_variable_type
  
  TYPE field_variable_ptr_type
    TYPE(field_variable_type), POINTER :: ptr
  END TYPE field_variable_ptr_type

  TYPE field_create_values_cache_type
    LOGICAL :: label_locked
    LOGICAL :: decomposition_locked
    LOGICAL :: dataprojectionlocked
    LOGICAL :: dependent_type_locked
    LOGICAL :: number_of_variables_locked
    LOGICAL :: geometric_field_locked
    LOGICAL :: scaling_type_locked
    LOGICAL :: type_locked
    INTEGER(INTG), ALLOCATABLE :: variable_types(:)
    LOGICAL :: variable_types_locked
    TYPE(varying_string), ALLOCATABLE :: variable_labels(:)
    LOGICAL, ALLOCATABLE :: variable_labels_locked(:)
    INTEGER(INTG), ALLOCATABLE :: DIMENSION(:)
    LOGICAL, ALLOCATABLE :: dimension_locked(:)
    INTEGER(INTG), ALLOCATABLE :: data_types(:)
    LOGICAL, ALLOCATABLE :: data_types_locked(:)
    INTEGER(INTG), ALLOCATABLE :: dof_order_types(:)
    LOGICAL, ALLOCATABLE :: dof_order_types_locked(:)
    INTEGER(INTG), ALLOCATABLE :: number_of_components(:)
    LOGICAL, ALLOCATABLE :: number_of_components_locked(:)
    TYPE(varying_string), ALLOCATABLE :: component_labels(:,:)
    LOGICAL, ALLOCATABLE :: component_labels_locked(:,:)
    INTEGER(INTG), ALLOCATABLE :: interpolation_type(:,:)
    LOGICAL, ALLOCATABLE :: interpolation_type_locked(:,:)
    INTEGER(INTG), ALLOCATABLE :: mesh_component_number(:,:)
    LOGICAL, ALLOCATABLE :: mesh_component_number_locked(:,:)
  END TYPE field_create_values_cache_type

  TYPE field_type
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(varying_string) :: label
    LOGICAL :: field_finished
    TYPE(fields_type), POINTER :: fields
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: TYPE
    INTEGER(INTG) :: dependent_type
    TYPE(decomposition_type), POINTER :: decomposition
    INTEGER(INTG) :: number_of_variables
    TYPE(field_variable_ptr_type), ALLOCATABLE :: variable_type_map(:)
    TYPE(field_variable_type), ALLOCATABLE :: variables(:)
    TYPE(field_scalings_type) :: scalings
    TYPE(field_type), POINTER :: geometric_field
    TYPE(field_geometric_parameters_type), POINTER :: geometric_field_parameters
    TYPE(field_create_values_cache_type), POINTER :: create_values_cache
    TYPE(data_projection_type), POINTER :: dataprojection
  END TYPE field_type

  TYPE field_ptr_type
    TYPE(field_type), POINTER :: ptr
  END TYPE field_ptr_type

  TYPE fields_type
    TYPE(region_type), POINTER :: region
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: number_of_fields
    TYPE(field_ptr_type), POINTER :: fields(:)
  END TYPE fields_type

  !
  !================================================================================================================================
  !
  ! Equations matrices types
  !
  
  TYPE element_matrix_type
    INTEGER(INTG) :: equations_matrix_number
    INTEGER(INTG) :: structure_type
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: number_of_columns
    INTEGER(INTG) :: max_number_of_rows
    INTEGER(INTG) :: max_number_of_columns
    INTEGER(INTG), ALLOCATABLE :: row_dofs(:)
    INTEGER(INTG), ALLOCATABLE :: column_dofs(:)
    REAL(DP), ALLOCATABLE :: matrix(:,:)
  END TYPE element_matrix_type

  TYPE element_vector_type
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: max_number_of_rows
    INTEGER(INTG), ALLOCATABLE :: row_dofs(:)
    REAL(DP), ALLOCATABLE :: vector(:)
  END TYPE element_vector_type

  TYPE nodalmatrixtype
    INTEGER(INTG) :: equationsmatrixnumber
    INTEGER(INTG) :: structuretype
    INTEGER(INTG) :: numberofrows
    INTEGER(INTG) :: numberofcolumns
    INTEGER(INTG) :: maxnumberofrows
    INTEGER(INTG) :: maxnumberofcolumns
    INTEGER(INTG), ALLOCATABLE :: rowdofs(:)
    INTEGER(INTG), ALLOCATABLE :: columndofs(:)
    REAL(DP), ALLOCATABLE :: matrix(:,:)
  END TYPE nodalmatrixtype

  TYPE nodalvectortype
    INTEGER(INTG) :: numberofrows
    INTEGER(INTG) :: maxnumberofrows
    INTEGER(INTG), ALLOCATABLE :: rowdofs(:)
    REAL(DP), ALLOCATABLE :: vector(:)
  END TYPE nodalvectortype

  TYPE equations_matrix_type
    INTEGER(INTG) :: matrix_number
    TYPE(equations_matrices_dynamic_type), POINTER :: dynamic_matrices
    TYPE(equations_matrices_linear_type), POINTER :: linear_matrices
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: structure_type
    LOGICAL :: lumped
    INTEGER(INTG) :: number_of_columns
    LOGICAL :: update_matrix
    LOGICAL :: first_assembly
    TYPE(distributed_matrix_type), POINTER :: matrix
    TYPE(element_matrix_type) :: element_matrix
    TYPE(nodalmatrixtype) :: nodalmatrix
    TYPE(distributed_vector_type), POINTER :: temp_vector
  END TYPE equations_matrix_type

  TYPE equations_matrix_ptr_type
    TYPE(equations_matrix_type), POINTER :: ptr
  END TYPE equations_matrix_ptr_type
 
  TYPE equations_jacobian_type
    INTEGER(INTG) :: jacobian_number
    TYPE(equations_matrices_nonlinear_type), POINTER :: nonlinear_matrices
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: structure_type
    INTEGER(INTG) :: number_of_columns
    LOGICAL :: update_jacobian
    TYPE(distributed_matrix_type), POINTER :: jacobian
    LOGICAL :: first_assembly
    TYPE(element_matrix_type) :: element_jacobian
    TYPE(nodalmatrixtype) :: nodaljacobian
    INTEGER(INTG) :: jacobian_calculation_type
  END TYPE equations_jacobian_type

  TYPE equations_jacobian_ptr_type
    TYPE(equations_jacobian_type), POINTER :: ptr
  END TYPE equations_jacobian_ptr_type

  TYPE equations_matrices_dynamic_type
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    INTEGER(INTG) :: number_of_dynamic_matrices
    TYPE(equations_matrix_ptr_type), ALLOCATABLE :: matrices(:)
    TYPE(distributed_vector_type), POINTER :: temp_vector
  END TYPE equations_matrices_dynamic_type

  TYPE equations_matrices_linear_type
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    INTEGER(INTG) :: number_of_linear_matrices
    TYPE(equations_matrix_ptr_type), ALLOCATABLE :: matrices(:)
  END TYPE equations_matrices_linear_type

  TYPE equations_matrices_nonlinear_type
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    INTEGER(INTG) :: number_of_jacobians
    TYPE(equations_jacobian_ptr_type), ALLOCATABLE :: jacobians(:)
    LOGICAL :: update_residual
    LOGICAL :: first_assembly
    TYPE(distributed_vector_type), POINTER :: residual
    TYPE(element_vector_type) :: element_residual
    TYPE(nodalvectortype) :: nodalresidual
    INTEGER(INTG) :: nodalresidualcalculated
    INTEGER(INTG) :: element_residual_calculated
  END TYPE equations_matrices_nonlinear_type

  TYPE equations_matrices_rhs_type
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    LOGICAL :: update_vector
    LOGICAL :: first_assembly
    TYPE(distributed_vector_type), POINTER :: vector
    TYPE(element_vector_type) :: element_vector
    TYPE(nodalvectortype) :: nodalvector
  END TYPE equations_matrices_rhs_type
  
  TYPE equations_matrices_source_type
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    LOGICAL :: update_vector
    LOGICAL :: first_assembly
    TYPE(distributed_vector_type), POINTER :: vector
    TYPE(element_vector_type) :: element_vector
    TYPE(nodalvectortype) :: nodalvector
  END TYPE equations_matrices_source_type
  
  TYPE equations_matrices_type
    TYPE(equations_type), POINTER :: equations
    LOGICAL :: equations_matrices_finished
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: total_number_of_rows
    INTEGER(INTG) :: number_of_global_rows
    !Equations matrices components
    TYPE(equations_matrices_dynamic_type), POINTER :: dynamic_matrices
    TYPE(equations_matrices_linear_type), POINTER :: linear_matrices
    TYPE(equations_matrices_nonlinear_type), POINTER :: nonlinear_matrices
    TYPE(equations_matrices_rhs_type), POINTER :: rhs_vector
    TYPE(equations_matrices_source_type), POINTER :: source_vector
  END TYPE equations_matrices_type

  !
  !================================================================================================================================
  !
  ! Equations mapping types
  !
  
  TYPE var_to_equations_column_map_type
    INTEGER(INTG), ALLOCATABLE :: column_dof(:)
  END TYPE var_to_equations_column_map_type

  TYPE var_to_equations_matrices_map_type
    INTEGER(INTG) :: variable_index
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG) :: number_of_equations_matrices
    INTEGER(INTG), ALLOCATABLE :: equations_matrix_numbers(:)
    TYPE(var_to_equations_column_map_type), ALLOCATABLE :: dof_to_columns_maps(:)
    INTEGER(INTG), ALLOCATABLE :: dof_to_rows_map(:)
  END TYPE var_to_equations_matrices_map_type

  TYPE equations_matrix_to_var_map_type
    INTEGER(INTG) :: matrix_number
    TYPE(equations_matrix_type), POINTER :: equations_matrix
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG) :: number_of_columns
    REAL(DP) :: matrix_coefficient
    INTEGER(INTG), ALLOCATABLE :: column_to_dof_map(:)
    TYPE(domain_mapping_type), POINTER :: column_dofs_mapping
  END TYPE equations_matrix_to_var_map_type

  TYPE equations_mapping_dynamic_type
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    INTEGER(INTG) :: number_of_dynamic_equations_matrices

    INTEGER(INTG) :: stiffness_matrix_number
    INTEGER(INTG) :: damping_matrix_number
    INTEGER(INTG) :: mass_matrix_number
    INTEGER(INTG) :: dynamic_variable_type
    TYPE(field_variable_type), POINTER :: dynamic_variable
    TYPE(var_to_equations_matrices_map_type), ALLOCATABLE :: var_to_equations_matrices_maps(:)
    TYPE(equations_matrix_to_var_map_type), ALLOCATABLE :: equations_matrix_to_var_maps(:)
    INTEGER(INTG), ALLOCATABLE :: equations_row_to_variable_dof_maps(:)
  END TYPE equations_mapping_dynamic_type

  TYPE equations_mapping_linear_type
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    INTEGER(INTG) :: number_of_linear_equations_matrices
    INTEGER(INTG) :: number_of_linear_matrix_variables
    INTEGER(INTG), ALLOCATABLE :: linear_matrix_variable_types(:)
    TYPE(var_to_equations_matrices_map_type), ALLOCATABLE :: var_to_equations_matrices_maps(:)
    TYPE(equations_matrix_to_var_map_type), ALLOCATABLE :: equations_matrix_to_var_maps(:)
    INTEGER(INTG), ALLOCATABLE :: equations_row_to_variable_dof_maps(:,:)
  END TYPE equations_mapping_linear_type

  TYPE equations_jacobian_to_var_map_type
    INTEGER(INTG) :: jacobian_number
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    TYPE(equations_jacobian_type), POINTER :: jacobian
    INTEGER(INTG) :: number_of_columns
    REAL(DP) :: jacobian_coefficient
    INTEGER(INTG), ALLOCATABLE :: equations_column_to_dof_variable_map(:)
    TYPE(domain_mapping_type), POINTER :: column_dofs_mapping
  END TYPE equations_jacobian_to_var_map_type

  TYPE var_to_equations_jacobian_map_type
    INTEGER(INTG) :: jacobian_number
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG), ALLOCATABLE :: dof_to_columns_map(:)
    INTEGER(INTG), ALLOCATABLE :: dof_to_rows_map(:)
  END TYPE var_to_equations_jacobian_map_type

  TYPE equations_mapping_nonlinear_type
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    INTEGER(INTG) :: number_of_residual_variables
    TYPE(field_variable_ptr_type), ALLOCATABLE :: residual_variables(:)
    TYPE(var_to_equations_jacobian_map_type), ALLOCATABLE :: var_to_jacobian_map(:)
    TYPE(equations_jacobian_to_var_map_type), ALLOCATABLE :: jacobian_to_var_map(:)
    REAL(DP) :: residual_coefficient
    INTEGER(INTG), ALLOCATABLE :: equations_row_to_residual_dof_map(:)
  END TYPE equations_mapping_nonlinear_type

  TYPE equations_mapping_rhs_type
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    INTEGER(INTG) :: rhs_variable_type
    TYPE(field_variable_type), POINTER :: rhs_variable
    TYPE(domain_mapping_type), POINTER :: rhs_variable_mapping
    REAL(DP) :: rhs_coefficient
    INTEGER(INTG), ALLOCATABLE :: rhs_dof_to_equations_row_map(:)
    INTEGER(INTG), ALLOCATABLE :: equations_row_to_rhs_dof_map(:)
  END TYPE equations_mapping_rhs_type

  TYPE equations_mapping_source_type
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    INTEGER(INTG) :: source_variable_type
    TYPE(field_variable_type), POINTER :: source_variable
    TYPE(domain_mapping_type), POINTER :: source_variable_mapping
    REAL(DP) :: source_coefficient
    INTEGER(INTG), ALLOCATABLE :: source_dof_to_equations_row_map(:)
    INTEGER(INTG), ALLOCATABLE :: equations_row_to_source_dof_map(:)
  END TYPE equations_mapping_source_type

  TYPE equations_mapping_create_values_cache_type
    INTEGER(INTG) :: number_of_dynamic_equations_matrices
    INTEGER(INTG) :: dynamic_stiffness_matrix_number
    INTEGER(INTG) :: dynamic_damping_matrix_number
    INTEGER(INTG) :: dynamic_mass_matrix_number
    INTEGER(INTG) :: dynamic_variable_type
    REAL(DP), ALLOCATABLE :: dynamic_matrix_coefficients(:)
    INTEGER(INTG) :: number_of_linear_equations_matrices
    INTEGER(INTG), ALLOCATABLE :: linear_matrix_variable_types(:)
    REAL(DP), ALLOCATABLE :: linear_matrix_coefficients(:)
    INTEGER(INTG) :: number_of_residual_variables
    INTEGER(INTG), ALLOCATABLE :: residual_variable_types(:)
    REAL(DP) :: residual_coefficient
    INTEGER(INTG) :: rhs_variable_type
    REAL(DP) :: rhs_coefficient
    INTEGER(INTG) :: source_variable_type
    REAL(DP) :: source_coefficient
  END TYPE equations_mapping_create_values_cache_type

  TYPE equations_mapping_type
    TYPE(equations_type), POINTER :: equations
    LOGICAL :: equations_mapping_finished
    TYPE(equations_matrices_type), POINTER :: equations_matrices
    !Row mappings
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: total_number_of_rows
    INTEGER(INTG) :: number_of_global_rows
    TYPE(domain_mapping_type), POINTER :: row_dofs_mapping
    !Equations mapping components
    TYPE(equations_mapping_dynamic_type), POINTER :: dynamic_mapping
    TYPE(equations_mapping_linear_type), POINTER :: linear_mapping
    TYPE(equations_mapping_nonlinear_type), POINTER :: nonlinear_mapping
    TYPE(equations_mapping_rhs_type), POINTER :: rhs_mapping
    TYPE(equations_mapping_source_type), POINTER :: source_mapping
    !Create values cache
    TYPE(equations_mapping_create_values_cache_type), POINTER :: create_values_cache
  END TYPE equations_mapping_type

  !
  !================================================================================================================================
  !
  ! Equations types
  !
  
  TYPE equations_interpolation_type
    TYPE(equations_type), POINTER :: equations
    TYPE(field_type), POINTER :: geometric_field
    TYPE(field_type), POINTER :: fibre_field
    TYPE(field_type), POINTER :: dependent_field
    TYPE(field_type), POINTER :: independent_field
    TYPE(field_type), POINTER :: materials_field
    TYPE(field_type), POINTER :: source_field
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: geometric_interp_parameters(:)
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: fibre_interp_parameters(:)
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: dependent_interp_parameters(:)
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: independent_interp_parameters(:)
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: materials_interp_parameters(:)
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: source_interp_parameters(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: geometric_interp_point(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: fibre_interp_point(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: dependent_interp_point(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: independent_interp_point(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: materials_interp_point(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: source_interp_point(:)
    TYPE(field_physical_point_ptr_type), POINTER :: dependent_physical_point(:)
     TYPE(field_interpolated_point_metrics_ptr_type), POINTER :: dependent_interp_point_metrics(:)
    TYPE(field_interpolated_point_metrics_ptr_type), POINTER :: independent_interp_point_metrics(:)
    TYPE(field_interpolated_point_metrics_ptr_type), POINTER :: geometric_interp_point_metrics(:)
    TYPE(field_interpolated_point_metrics_ptr_type), POINTER :: fibre_interp_point_metrics(:)
  END TYPE equations_interpolation_type

  TYPE equations_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: equations_finished
    INTEGER(INTG) :: linearity
    INTEGER(INTG) :: time_dependence
    INTEGER(INTG) :: output_type
    INTEGER(INTG) :: sparsity_type
    INTEGER(INTG) :: lumping_type
    TYPE(equations_interpolation_type), POINTER :: interpolation
    TYPE(equations_mapping_type), POINTER :: equations_mapping
    TYPE(equations_matrices_type), POINTER :: equations_matrices
  END TYPE equations_type

  TYPE equations_ptr_type
    TYPE(equations_type), POINTER :: ptr
  END TYPE equations_ptr_type

  !
  !================================================================================================================================
  !
  ! Boundary conditions types
  !

  TYPE boundary_conditions_variable_type
    TYPE(boundary_conditions_type), POINTER :: boundary_conditions
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG), ALLOCATABLE :: dof_types(:)
    INTEGER(INTG), ALLOCATABLE :: condition_types(:)
    TYPE(boundary_conditions_dirichlet_type), POINTER :: dirichlet_boundary_conditions
    INTEGER(INTG) :: number_of_dirichlet_conditions
    TYPE(boundaryconditionsneumanntype), POINTER :: neumannboundaryconditions
    TYPE(boundary_conditions_pressure_incremented_type), POINTER :: pressure_incremented_boundary_conditions
    INTEGER(INTG), ALLOCATABLE :: dof_counts(:)
    LOGICAL, ALLOCATABLE :: parametersetrequired(:)
    TYPE(boundaryconditionsdofconstraintstype), POINTER :: dofconstraints
  END TYPE boundary_conditions_variable_type

  TYPE boundary_conditions_variable_ptr_type
    TYPE(boundary_conditions_variable_type), POINTER :: ptr
  END TYPE boundary_conditions_variable_ptr_type

  TYPE boundary_conditions_type
    TYPE(solver_equations_type), POINTER :: solver_equations
    LOGICAL :: boundary_conditions_finished
    INTEGER(INTG) :: number_of_boundary_conditions_variables
    TYPE(boundary_conditions_variable_ptr_type), ALLOCATABLE :: boundary_conditions_variables(:)
    INTEGER(INTG) :: neumannmatrixsparsity
  END TYPE boundary_conditions_type

  TYPE boundary_conditions_sparsity_indices_ptr_type
    TYPE(boundary_conditions_sparsity_indices_type), POINTER :: ptr
  END TYPE boundary_conditions_sparsity_indices_ptr_type

  TYPE boundary_conditions_dirichlet_type
    INTEGER(INTG), ALLOCATABLE :: dirichlet_dof_indices(:)
    TYPE(boundary_conditions_sparsity_indices_ptr_type), ALLOCATABLE :: linear_sparsity_indices(:,:)
    TYPE(boundary_conditions_sparsity_indices_ptr_type), ALLOCATABLE :: dynamic_sparsity_indices(:,:)
  END TYPE boundary_conditions_dirichlet_type

  TYPE boundary_conditions_sparsity_indices_type
    INTEGER(INTG), ALLOCATABLE :: sparse_row_indices(:)
    INTEGER(INTG), ALLOCATABLE :: sparse_column_indices(:)
  END TYPE boundary_conditions_sparsity_indices_type

  TYPE boundaryconditionsneumanntype
    INTEGER(INTG), ALLOCATABLE :: setdofs(:)
    TYPE(distributed_matrix_type), POINTER :: integrationmatrix
    TYPE(distributed_vector_type), POINTER :: pointvalues
    TYPE(domain_mapping_type), POINTER :: pointdofmapping
  END TYPE boundaryconditionsneumanntype

  TYPE boundary_conditions_pressure_incremented_type
    INTEGER(INTG), ALLOCATABLE :: pressure_incremented_dof_indices(:)
  END TYPE boundary_conditions_pressure_incremented_type

  TYPE boundaryconditionsdofconstrainttype
    INTEGER(INTG) :: globaldof
    INTEGER(INTG) :: numberofdofs
    INTEGER(INTG), ALLOCATABLE :: dofs(:)
    REAL(DP), ALLOCATABLE :: coefficients(:)
  END TYPE boundaryconditionsdofconstrainttype

  TYPE boundaryconditionsdofconstraintptrtype
    TYPE(boundaryconditionsdofconstrainttype), POINTER :: ptr
  END TYPE boundaryconditionsdofconstraintptrtype

  TYPE boundaryconditionscoupleddofstype
    INTEGER(INTG) :: numberofdofs
    INTEGER(INTG), ALLOCATABLE :: globaldofs(:)
    INTEGER(INTG), ALLOCATABLE :: localdofs(:)
    REAL(DP), ALLOCATABLE :: coefficients(:)
  END TYPE boundaryconditionscoupleddofstype

  TYPE boundaryconditionscoupleddofsptrtype
    TYPE(boundaryconditionscoupleddofstype), POINTER :: ptr
  END TYPE boundaryconditionscoupleddofsptrtype

  TYPE boundaryconditionsdofconstraintstype
    INTEGER(INTG) :: numberofconstraints
    INTEGER(INTG) :: numberofdofs
    TYPE(boundaryconditionsdofconstraintptrtype), ALLOCATABLE :: constraints(:)
    TYPE(boundaryconditionscoupleddofsptrtype), ALLOCATABLE :: dofcouplings(:)
  END TYPE boundaryconditionsdofconstraintstype

  !
  !================================================================================================================================
  !
  ! Equations set types
  !

  TYPE equations_set_setup_type
    INTEGER(INTG) :: setup_type
    INTEGER(INTG) :: action_type
    INTEGER(INTG) :: field_user_number
    TYPE(field_type), POINTER :: field
    INTEGER(INTG) :: analytic_function_type
  END TYPE equations_set_setup_type  
  
  TYPE equations_set_geometry_type
    TYPE(equations_set_type), POINTER :: equations_set
    TYPE(field_type), POINTER :: geometric_field
    TYPE(field_type), POINTER :: fibre_field
  END TYPE equations_set_geometry_type

  TYPE equations_set_materials_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: materials_finished
    LOGICAL :: materials_field_auto_created
    TYPE(field_type), POINTER :: materials_field
  END TYPE equations_set_materials_type

  TYPE equations_set_dependent_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: dependent_finished
    LOGICAL :: dependent_field_auto_created
    TYPE(field_type), POINTER :: dependent_field
  END TYPE equations_set_dependent_type

  TYPE equationssetderivedtype
    TYPE(equations_set_type), POINTER :: equationsset
    LOGICAL :: derivedfinished
    LOGICAL :: derivedfieldautocreated
    TYPE(field_type), POINTER :: derivedfield
    INTEGER(INTG) :: numberofvariables
    INTEGER(INTG), ALLOCATABLE :: variabletypes(:)
  END TYPE equationssetderivedtype

  TYPE equations_set_independent_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: independent_finished
    LOGICAL :: independent_field_auto_created
    TYPE(field_type), POINTER :: independent_field
  END TYPE equations_set_independent_type

  TYPE equations_set_source_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: source_finished
    LOGICAL :: source_field_auto_created
    TYPE(field_type), POINTER :: source_field
  END TYPE equations_set_source_type
  
  TYPE equations_set_analytic_type
    TYPE(equations_set_type), POINTER :: equations_set
    INTEGER(INTG) :: analytic_function_type
    LOGICAL :: analytic_finished
    LOGICAL :: analytic_field_auto_created
    TYPE(field_type), POINTER :: analytic_field
    REAL(DP) :: analytic_time
    REAL(DP) :: analytic_user_params(20)
  END TYPE equations_set_analytic_type

  TYPE equations_set_equations_set_field_type
    TYPE(equations_set_type), POINTER :: equations_set
    LOGICAL :: equations_set_field_finished
    LOGICAL :: equations_set_field_auto_created
    TYPE(field_type), POINTER :: equations_set_field_field
  END TYPE equations_set_equations_set_field_type

  TYPE equations_set_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    LOGICAL :: equations_set_finished
    TYPE(equations_sets_type), POINTER :: equations_sets
    TYPE(region_type), POINTER :: region
    INTEGER(INTG), ALLOCATABLE :: specification(:)
    INTEGER(INTG) :: solution_method
    TYPE(equations_set_geometry_type) :: geometry
    TYPE(equations_set_materials_type), POINTER :: materials
    TYPE(equations_set_source_type), POINTER :: source
    TYPE(equations_set_dependent_type) :: dependent
    TYPE(equations_set_independent_type), POINTER :: independent
    TYPE(equations_set_analytic_type), POINTER :: analytic
    TYPE(equationssetderivedtype), POINTER :: derived
    TYPE(equations_type), POINTER :: equations
    TYPE(boundary_conditions_type), POINTER :: boundary_conditions
    TYPE(equations_set_equations_set_field_type) :: equations_set_field
  END TYPE equations_set_type
  
  TYPE equations_set_ptr_type
    TYPE(equations_set_type), POINTER :: ptr
  END TYPE equations_set_ptr_type
  
  TYPE equations_sets_type
    TYPE(region_type) , POINTER :: region
    INTEGER(INTG) :: number_of_equations_sets
    TYPE(equations_set_ptr_type), POINTER :: equations_sets(:)
  END TYPE equations_sets_type

  !
  !================================================================================================================================
  !
  ! Interface types
  
  TYPE interface_matrix_type
    TYPE(interface_matrices_type), POINTER :: interface_matrices
    INTEGER(INTG) :: matrix_number
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: structure_type
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: total_number_of_rows
    INTEGER(INTG) :: interface_matrix_time_dependence_type
    INTEGER(INTG) :: interface_matrix_transpose_time_dependence_type
    LOGICAL :: update_matrix
    LOGICAL :: first_assembly
    LOGICAL :: has_transpose
    TYPE(distributed_matrix_type), POINTER :: matrix
    TYPE(distributed_matrix_type), POINTER :: matrix_transpose
    TYPE(distributed_vector_type), POINTER :: temp_vector
    TYPE(distributed_vector_type), POINTER :: temp_transpose_vector
    TYPE(element_matrix_type) :: element_matrix
  END TYPE interface_matrix_type

  TYPE interface_matrix_ptr_type
    TYPE(interface_matrix_type), POINTER :: ptr
  END TYPE interface_matrix_ptr_type

  TYPE interface_rhs_type
    TYPE(interface_matrices_type), POINTER :: interface_matrices
    LOGICAL :: update_vector
    LOGICAL :: first_assembly
    TYPE(distributed_vector_type), POINTER :: rhs_vector
    TYPE(element_vector_type) :: element_vector
  END TYPE interface_rhs_type
  
  TYPE interface_matrices_type
    TYPE(interface_equations_type), POINTER :: interface_equations
    LOGICAL :: interface_matrices_finished
    TYPE(interface_mapping_type), POINTER :: interface_mapping
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    INTEGER(INTG) :: number_of_columns
    INTEGER(INTG) :: total_number_of_columns
    INTEGER(INTG) :: number_of_global_columns
    INTEGER(INTG) :: number_of_interface_matrices
    TYPE(interface_matrix_ptr_type), ALLOCATABLE :: matrices(:)
    TYPE(interface_rhs_type), POINTER :: rhs_vector
  END TYPE interface_matrices_type

  TYPE interface_matrix_to_var_map_type
    INTEGER(INTG) :: matrix_number
    TYPE(interface_matrix_type), POINTER :: interface_matrix
    TYPE(equations_set_type), POINTER :: equations_set
    TYPE(interface_equations_type), POINTER :: interface_equations
    INTEGER(INTG) :: variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG) :: mesh_index
    REAL(DP) :: matrix_coefficient
    LOGICAL :: has_transpose
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: total_number_of_rows
    INTEGER(INTG) :: number_of_global_rows
    TYPE(domain_mapping_type), POINTER :: row_dofs_mapping
    INTEGER(INTG), ALLOCATABLE :: variable_dof_to_row_map(:)
  END TYPE interface_matrix_to_var_map_type

  TYPE interface_mapping_rhs_type
    TYPE(interface_mapping_type), POINTER :: interface_mapping
    INTEGER(INTG) :: rhs_variable_type
    TYPE(field_variable_type), POINTER :: rhs_variable
    TYPE(domain_mapping_type), POINTER :: rhs_variable_mapping
    REAL(DP) :: rhs_coefficient
    INTEGER(INTG), ALLOCATABLE :: rhs_dof_to_interface_row_map(:)
    INTEGER(INTG), ALLOCATABLE :: interface_row_to_rhs_dof_map(:)
  END TYPE interface_mapping_rhs_type
  
  TYPE interface_mapping_create_values_cache_type
    INTEGER(INTG) :: number_of_interface_matrices
    INTEGER(INTG) :: lagrange_variable_type
    REAL(DP), ALLOCATABLE :: matrix_coefficients(:)
    LOGICAL, ALLOCATABLE :: has_transpose(:)
    INTEGER(INTG), ALLOCATABLE :: matrix_row_field_variable_indices(:)
    INTEGER(INTG), ALLOCATABLE :: matrix_col_field_variable_indices(:)
    INTEGER(INTG) :: rhs_lagrange_variable_type
    REAL(DP) :: rhs_coefficient
  END TYPE interface_mapping_create_values_cache_type
  
  TYPE interface_mapping_type  
    TYPE(interface_equations_type), POINTER :: interface_equations
    LOGICAL :: interface_mapping_finished
    INTEGER(INTG) :: lagrange_variable_type
    TYPE(field_variable_type), POINTER :: lagrange_variable
    INTEGER(INTG) :: number_of_columns
    INTEGER(INTG) :: total_number_of_columns
    INTEGER(INTG) :: number_of_global_columns
    TYPE(domain_mapping_type), POINTER :: column_dofs_mapping
    INTEGER(INTG), ALLOCATABLE :: lagrange_dof_to_column_map(:)
    INTEGER(INTG) :: number_of_interface_matrices
    TYPE(interface_matrix_to_var_map_type), ALLOCATABLE :: interface_matrix_rows_to_var_maps(:)
    TYPE(interface_mapping_rhs_type), POINTER :: rhs_mapping
    TYPE(interface_mapping_create_values_cache_type), POINTER :: create_values_cache
  END TYPE interface_mapping_type

  TYPE interface_equations_interpolation_set_type
    TYPE(field_interpolation_parameters_ptr_type), POINTER :: interpolation_parameters(:)
    TYPE(field_interpolated_point_ptr_type), POINTER :: interpolated_point(:)
    TYPE(field_interpolated_point_metrics_ptr_type), POINTER :: interpolated_point_metrics(:)
  END TYPE interface_equations_interpolation_set_type

  TYPE interface_equations_domain_interpolation_type
    TYPE(interface_equations_interpolation_type), POINTER :: interpolation
    TYPE(field_type), POINTER :: geometric_field
    INTEGER(INTG) :: number_of_geometric_interpolation_sets
    TYPE(interface_equations_interpolation_set_type), ALLOCATABLE :: geometric_interpolation(:)
    TYPE(field_type), POINTER :: dependent_field
    INTEGER(INTG) :: number_of_dependent_interpolation_sets
    TYPE(interface_equations_interpolation_set_type), ALLOCATABLE :: dependent_interpolation(:)
    TYPE(field_type), POINTER :: penalty_field
    INTEGER(INTG) :: number_of_penalty_interpolation_sets
    TYPE(interface_equations_interpolation_set_type), ALLOCATABLE :: penalty_interpolation(:)
  END TYPE interface_equations_domain_interpolation_type
  
  TYPE interface_equations_interpolation_type
    TYPE(interface_equations_type), POINTER :: interface_equations
    TYPE(interface_equations_domain_interpolation_type) :: interface_interpolation
    TYPE(interface_equations_domain_interpolation_type), ALLOCATABLE :: variable_interpolation(:)
  END TYPE interface_equations_interpolation_type

  TYPE interface_equations_type
    TYPE(interface_condition_type), POINTER :: interface_condition
    LOGICAL :: interface_equations_finished
    INTEGER(INTG) :: output_type
    INTEGER(INTG) :: sparsity_type
    INTEGER(INTG) :: linearity
    INTEGER(INTG) :: time_dependence
    TYPE(interface_equations_interpolation_type), POINTER :: interpolation
    TYPE(interface_mapping_type), POINTER :: interface_mapping
    TYPE(interface_matrices_type), POINTER :: interface_matrices
  END TYPE interface_equations_type
  
  TYPE interface_geometry_type
    TYPE(interface_condition_type), POINTER :: interface_condition
    TYPE(field_type), POINTER :: geometric_field
  END TYPE interface_geometry_type

  TYPE interface_penalty_type
    TYPE(interface_condition_type), POINTER :: interface_condition
    LOGICAL :: penalty_finished
    LOGICAL :: penalty_field_auto_created
    TYPE(field_type), POINTER :: penalty_field
  END TYPE interface_penalty_type

  TYPE interface_lagrange_type
    TYPE(interface_condition_type), POINTER :: interface_condition
    LOGICAL :: lagrange_finished
    LOGICAL :: lagrange_field_auto_created
    TYPE(field_type), POINTER :: lagrange_field
    INTEGER(INTG) :: number_of_components
  END TYPE interface_lagrange_type

  TYPE interface_dependent_type
    TYPE(interface_condition_type), POINTER :: interface_condition
    INTEGER(INTG) :: number_of_dependent_variables
    TYPE(equations_set_ptr_type), POINTER :: equations_sets(:)
    TYPE(field_variable_ptr_type), POINTER :: field_variables(:)
    INTEGER(INTG), POINTER :: variable_mesh_indices(:)
  END TYPE interface_dependent_type

  TYPE interface_condition_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    LOGICAL :: interface_condition_finished
    TYPE(interface_conditions_type), POINTER :: interface_conditions
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: method
    INTEGER(INTG) :: operator
    INTEGER(INTG) :: integrationtype
    TYPE(interface_geometry_type) :: geometry
    TYPE(interface_penalty_type), POINTER :: penalty
    TYPE(interface_lagrange_type), POINTER :: lagrange
    TYPE(interface_dependent_type), POINTER :: dependent
    TYPE(interface_equations_type), POINTER :: interface_equations
    TYPE(boundary_conditions_type), POINTER :: boundary_conditions
  END TYPE interface_condition_type

  TYPE interface_condition_ptr_type
    TYPE(interface_condition_type), POINTER :: ptr
  END TYPE interface_condition_ptr_type
  
  TYPE interface_conditions_type
    TYPE(interface_type), POINTER :: interface
    INTEGER(INTG) :: number_of_interface_conditions
    TYPE(interface_condition_ptr_type), POINTER :: interface_conditions(:)
  END TYPE interface_conditions_type
  
  TYPE interface_element_connectivity_type
    INTEGER(INTG) :: coupled_mesh_element_number
    REAL(DP), ALLOCATABLE :: xi(:,:,:)
    INTEGER(INTG) :: connected_face
    INTEGER(INTG) :: connected_line
  END TYPE interface_element_connectivity_type

  TYPE interface_mesh_connectivity_type
    TYPE(interface_type), POINTER :: interface
    TYPE(mesh_type), POINTER :: interface_mesh
    TYPE(basis_type), POINTER :: basis
    LOGICAL :: mesh_connectivity_finished
    INTEGER(INTG) :: number_of_interface_elements
    INTEGER(INTG) :: number_of_coupled_meshes
    TYPE(interface_element_connectivity_type), ALLOCATABLE :: element_connectivity(:,:)
  END TYPE interface_mesh_connectivity_type
  
  TYPE interfacepointconnectivitytype
    INTEGER(INTG) :: coupledmeshelementnumber
    INTEGER(INTG) :: elementlinefacenumber
    REAL(DP), ALLOCATABLE :: xi(:)
    REAL(DP), ALLOCATABLE :: reducedxi(:)
  END TYPE interfacepointconnectivitytype
  
  !Contains information on coupled mesh elements that are connected to each interface element.
  TYPE interfacecoupledelementstype
    INTEGER(INTG) :: numberofcoupledelements
    INTEGER(INTG), ALLOCATABLE :: elementnumbers(:)
  END TYPE interfacecoupledelementstype
  
  TYPE interfacepointsconnectivitytype
    TYPE(interface_type), POINTER :: interface
    TYPE(mesh_type), POINTER :: interfacemesh
    LOGICAL :: pointsconnectivityfinished
    TYPE(interfacepointconnectivitytype), ALLOCATABLE :: pointsconnectivity(:,:)
    TYPE(interfacecoupledelementstype), ALLOCATABLE :: coupledelements(:,:)
    INTEGER(INTG), ALLOCATABLE :: maxnumberofcoupledelements(:)
  END TYPE interfacepointsconnectivitytype
 
  TYPE interface_type
    INTEGER :: user_number
    INTEGER :: global_number
    LOGICAL :: interface_finished
    TYPE(varying_string) :: label
    TYPE(coordinate_system_type), POINTER :: coordinate_system
    TYPE(interfaces_type), POINTER :: interfaces
    TYPE(region_type), POINTER :: parent_region
    INTEGER(INTG) :: number_of_coupled_meshes
    TYPE(mesh_ptr_type), POINTER :: coupled_meshes(:)
    TYPE(interface_mesh_connectivity_type), POINTER :: mesh_connectivity
    TYPE(interfacepointsconnectivitytype), POINTER :: pointsconnectivity
    TYPE(data_points_type), POINTER :: data_points
    TYPE(nodes_type), POINTER :: nodes
    TYPE(meshes_type), POINTER :: meshes
    TYPE(generated_meshes_type), POINTER :: generated_meshes
    TYPE(fields_type), POINTER :: fields
    TYPE(interface_conditions_type), POINTER :: interface_conditions
  END TYPE interface_type

  TYPE interface_ptr_type
    TYPE(interface_type), POINTER :: ptr
  END TYPE interface_ptr_type
  
  TYPE interfaces_type
    TYPE(region_type), POINTER :: parent_region
    INTEGER(INTG) :: number_of_interfaces
    TYPE(interface_ptr_type), POINTER :: interfaces(:)
  END TYPE interfaces_type

  !
  !================================================================================================================================
  !
  ! CellML types (belongs under field types?)

  TYPE cellml_model_type
    TYPE(cellml_type), POINTER :: cellml
    INTEGER(INTG) :: global_number
    TYPE(varying_string) :: model_id
    TYPE(c_ptr) :: ptr
    INTEGER(INTG) :: number_of_state
    TYPE(varying_string), ALLOCATABLE :: state_variable_id(:)
    INTEGER(INTG) :: number_of_intermediate
    TYPE(varying_string), ALLOCATABLE :: intermediate_variable_id(:)
    INTEGER(INTG) :: number_of_parameters
    TYPE(varying_string), ALLOCATABLE :: parameter_variable_id(:)
 END TYPE cellml_model_type

  TYPE cellml_model_ptr_type
    TYPE(cellml_model_type), POINTER :: ptr
  END TYPE cellml_model_ptr_type

  TYPE cellml_models_field_type
    TYPE(cellml_type), POINTER :: cellml
    LOGICAL :: models_field_finished
    LOGICAL :: models_field_auto_created
    TYPE(field_type), POINTER :: models_field
    INTEGER(INTG) :: only_one_model_index
  END TYPE cellml_models_field_type
  
  TYPE cellml_state_field_type
    TYPE(cellml_type), POINTER :: cellml
    LOGICAL :: state_field_finished
    LOGICAL :: state_field_auto_created
    TYPE(field_type), POINTER :: state_field
  END TYPE cellml_state_field_type
  
  TYPE cellml_intermediate_field_type
    TYPE(cellml_type), POINTER :: cellml
    LOGICAL :: intermediate_field_finished
    LOGICAL :: intermediate_field_auto_created
    TYPE(field_type), POINTER :: intermediate_field
  END TYPE cellml_intermediate_field_type
  
  TYPE cellml_parameters_field_type
    TYPE(cellml_type), POINTER :: cellml
    LOGICAL :: parameters_field_finished
    LOGICAL :: parameters_field_auto_created
    TYPE(field_type), POINTER :: parameters_field
  END TYPE cellml_parameters_field_type
 
  TYPE cellmlpetsccontexttype
    TYPE(solver_type), POINTER :: solver
    TYPE(cellml_type), POINTER :: cellml
    INTEGER(INTG) :: dofidx
    REAL(DP), POINTER :: rates(:)
    INTEGER(INTG), ALLOCATABLE :: ratesindices(:)
  END TYPE cellmlpetsccontexttype
  
  TYPE cellml_model_map_type
    INTEGER(INTG) :: cellml_map_type
    TYPE(field_type), POINTER :: field
    INTEGER(INTG) :: variable_type
    INTEGER(INTG) :: component_number
    INTEGER(INTG) :: field_parameter_set
    TYPE(varying_string) :: variable_id
    INTEGER(INTG) :: cellml_field_type
    INTEGER(INTG) :: cellml_variable_number
    INTEGER(INTG) :: cellml_parameter_set
  END TYPE cellml_model_map_type

  TYPE cellml_model_map_ptr_type
    TYPE(cellml_model_map_type), POINTER :: ptr
  END TYPE cellml_model_map_ptr_type

  TYPE cellml_model_maps_type
    INTEGER(INTG) :: number_of_fields_mapped_to
    TYPE(cellml_model_map_ptr_type), ALLOCATABLE :: fields_mapped_to(:)
    INTEGER(INTG) :: number_of_fields_mapped_from
    TYPE(cellml_model_map_ptr_type), ALLOCATABLE :: fields_mapped_from(:)
  END TYPE cellml_model_maps_type

  TYPE cellml_model_maps_ptr_type
    TYPE(cellml_model_maps_type), POINTER :: ptr
  END TYPE cellml_model_maps_ptr_type
  
  TYPE cellml_field_maps_type
    TYPE(cellml_type), POINTER :: cellml
    LOGICAL :: cellml_field_maps_finished
    TYPE(field_type), POINTER :: source_geometric_field
    TYPE(field_variable_type), POINTER :: source_field_variable
    TYPE(domain_type), POINTER :: source_field_domain
    INTEGER(INTG) :: source_field_interpolation_type
    TYPE(cellml_model_maps_ptr_type), ALLOCATABLE :: model_maps(:)
    !INTEGER(INTG) :: NUMBER_OF_SOURCE_DOFS !<The number of local (excluding ghosts) source dofs.
    !INTEGER(INTG) :: TOTAL_NUMBER_OF_SOURCE_DOFS !<The number of local (including ghosts) source dofs.
    !INTEGER(INTG) :: GLOBAL_NUMBER_OF_SOURCE_DOFS !<The number of global source dofs.
  END TYPE cellml_field_maps_type
  
  TYPE cellml_type
    TYPE(region_type), POINTER :: region
    INTEGER(INTG) :: global_number
    INTEGER(INTG) :: user_number
    TYPE(cellml_environments_type), POINTER :: environments
    LOGICAL :: cellml_finished
    INTEGER(INTG) :: number_of_models
    INTEGER(INTG) :: maximum_number_of_state
    INTEGER(INTG) :: maximum_number_of_parameters
    INTEGER(INTG) :: maximum_number_of_intermediate
    TYPE(cellml_model_ptr_type), ALLOCATABLE :: models(:)
    TYPE(cellml_field_maps_type), POINTER :: field_maps
    TYPE(cellml_models_field_type), POINTER :: models_field
    TYPE(cellml_state_field_type), POINTER :: state_field
    TYPE(cellml_intermediate_field_type), POINTER :: intermediate_field
    TYPE(cellml_parameters_field_type), POINTER :: parameters_field
    LOGICAL :: cellml_generated
  END TYPE cellml_type

  TYPE cellml_ptr_type
    TYPE(cellml_type), POINTER :: ptr
  END TYPE cellml_ptr_type

  TYPE cellml_environments_type
    TYPE(region_type), POINTER :: region
    INTEGER(INTG) :: number_of_environments
    TYPE(cellml_ptr_type), ALLOCATABLE :: environments(:)
  END TYPE cellml_environments_type

  !
  !================================================================================================================================
  !
  ! Solver matrices types
  !
  
  TYPE solver_matrix_type
    INTEGER(INTG) :: matrix_number
    TYPE(solver_matrices_type), POINTER :: solver_matrices
    LOGICAL :: update_matrix
    INTEGER(INTG) :: storage_type
    INTEGER(INTG) :: number_of_columns
    TYPE(distributed_vector_type), POINTER :: solver_vector
    TYPE(distributed_matrix_type), POINTER :: matrix
  END TYPE solver_matrix_type

  TYPE solver_matrix_ptr_type
    TYPE(solver_matrix_type), POINTER :: ptr
  END TYPE solver_matrix_ptr_type
  
  TYPE solver_matrices_type
    TYPE(solver_equations_type), POINTER :: solver_equations
    LOGICAL :: solver_matrices_finished
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: number_of_global_rows
    INTEGER(INTG) :: library_type
    !Linear matrices
    INTEGER(INTG) :: number_of_matrices
    TYPE(solver_matrix_ptr_type), ALLOCATABLE :: matrices(:)
    !Nonlinear matrices and vectors
    LOGICAL :: update_residual
    TYPE(distributed_vector_type), POINTER :: residual
    !Right hand side vector
    LOGICAL :: update_rhs_vector
    TYPE(distributed_vector_type), POINTER :: rhs_vector
  END TYPE solver_matrices_type

  !
  !================================================================================================================================
  !
  ! Solver equations types
  !

  TYPE solver_equations_type
    TYPE(solver_type), POINTER :: solver
    LOGICAL :: solver_equations_finished

    INTEGER(INTG) :: linearity
    INTEGER(INTG) :: time_dependence

    INTEGER(INTG) :: sparsity_type

    TYPE(solver_mapping_type), POINTER :: solver_mapping
    TYPE(solver_matrices_type), POINTER :: solver_matrices

    TYPE(boundary_conditions_type), POINTER :: boundary_conditions

  END TYPE solver_equations_type

  !
  !================================================================================================================================
  !
  ! CellML equations types
  !

  TYPE cellml_equations_type
    TYPE(solver_type), POINTER :: solver
    LOGICAL :: cellml_equations_finished
    INTEGER(INTG) :: number_of_cellml_environments
    TYPE(cellml_ptr_type), ALLOCATABLE :: cellml_environments(:)
  END TYPE cellml_equations_type
  
  !
  !================================================================================================================================
  !
  ! Solver types
  !
  
  TYPE dynamic_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: solver_library
    LOGICAL :: solver_initialised
    INTEGER(INTG) :: linearity
    INTEGER(INTG) :: order
    INTEGER(INTG) :: degree
    INTEGER(INTG) :: scheme
    REAL(DP), ALLOCATABLE :: theta(:)
    LOGICAL :: explicit
    LOGICAL :: restart
    LOGICAL :: ale
    LOGICAL :: fsi
    LOGICAL :: update_bc
    REAL(DP) :: current_time
    REAL(DP) :: time_increment
    TYPE(solver_type), POINTER :: linear_solver
    TYPE(solver_type), POINTER :: nonlinear_solver
  END TYPE dynamic_solver_type
  
  TYPE forward_euler_dae_solver_type
    TYPE(euler_dae_solver_type), POINTER :: euler_dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE forward_euler_dae_solver_type

  TYPE backward_euler_dae_solver_type
    TYPE(euler_dae_solver_type), POINTER :: euler_dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE backward_euler_dae_solver_type

  TYPE improved_euler_dae_solver_type
    TYPE(euler_dae_solver_type), POINTER :: euler_dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE improved_euler_dae_solver_type
  
  TYPE euler_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: euler_type
    TYPE(forward_euler_dae_solver_type), POINTER :: forward_euler_solver
    TYPE(backward_euler_dae_solver_type), POINTER :: backward_euler_solver
    TYPE(improved_euler_dae_solver_type), POINTER :: improved_euler_solver
    INTEGER(INTG) :: solver_library
  END TYPE euler_dae_solver_type

  TYPE crank_nicolson_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE crank_nicolson_dae_solver_type
  
  TYPE runge_kutta_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE runge_kutta_dae_solver_type
  
  TYPE adams_moulton_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE adams_moulton_dae_solver_type
  
  TYPE bdf_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE bdf_dae_solver_type
  
  TYPE rush_larson_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
    INTEGER(INTG) :: solver_library
  END TYPE rush_larson_dae_solver_type
  
  TYPE external_dae_solver_type
    TYPE(dae_solver_type), POINTER :: dae_solver
  END TYPE external_dae_solver_type
  
  TYPE dae_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: dae_type
    INTEGER(INTG) :: dae_solve_type
    REAL(DP) :: start_time
    REAL(DP) :: end_time
    REAL(DP) :: initial_step
    TYPE(euler_dae_solver_type), POINTER :: euler_solver
    TYPE(crank_nicolson_dae_solver_type), POINTER :: crank_nicolson_solver
    TYPE(runge_kutta_dae_solver_type), POINTER :: runge_kutta_solver
    TYPE(adams_moulton_dae_solver_type), POINTER :: adams_moulton_solver
    TYPE(bdf_dae_solver_type), POINTER :: bdf_solver
    TYPE(rush_larson_dae_solver_type), POINTER :: rush_larson_solver
    TYPE(external_dae_solver_type), POINTER :: external_solver
  END TYPE dae_solver_type
  
  TYPE linear_direct_solver_type
    TYPE(linear_solver_type), POINTER :: linear_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    INTEGER(INTG) :: direct_solver_type
    TYPE(petscpctype) :: pc
    TYPE(petscksptype) :: ksp
  END TYPE linear_direct_solver_type

  TYPE linear_iterative_solver_type
    TYPE(linear_solver_type), POINTER :: linear_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    INTEGER(INTG) :: iterative_solver_type
    INTEGER(INTG) :: iterative_preconditioner_type
    INTEGER(INTG) :: solution_initialise_type
    INTEGER(INTG) :: maximum_number_of_iterations
    REAL(DP) :: relative_tolerance
    REAL(DP) :: absolute_tolerance
    REAL(DP) :: divergence_tolerance
    INTEGER(INTG) :: gmres_restart
    TYPE(petscpctype) :: pc
    TYPE(petscksptype) :: ksp
  END TYPE linear_iterative_solver_type
  
  TYPE linear_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: linear_solve_type
    LOGICAL :: linked_newton_petsc_solver
    TYPE(linear_direct_solver_type), POINTER :: direct_solver
    TYPE(linear_iterative_solver_type), POINTER :: iterative_solver
  END TYPE linear_solver_type

  TYPE newton_linesearch_solver_type
    TYPE(newton_solver_type), POINTER :: newton_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    INTEGER(INTG) :: linesearch_type
    REAL(DP) :: linesearch_alpha
    REAL(DP) :: linesearch_maxstep
    REAL(DP) :: linesearch_steptolerance
    TYPE(petscmatcoloringtype) :: jacobianmatcoloring
    TYPE(petsciscoloringtype) :: jacobianiscoloring
    TYPE(petscmatfdcoloringtype) :: jacobianmatfdcoloring
    TYPE(petscsnestype) :: snes
    TYPE(petscsneslinesearchtype) :: sneslinesearch
    LOGICAL :: linesearchmonitoroutput
  END TYPE newton_linesearch_solver_type
  
  TYPE newton_trustregion_solver_type
    TYPE(newton_solver_type), POINTER :: newton_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    REAL(DP) :: trustregion_tolerance
    REAL(DP) :: trustregion_delta0
    TYPE(petscsnestype) :: snes
  END TYPE newton_trustregion_solver_type

  TYPE newtonsolverconvergencetest
    REAL(DP) :: energyfirstiter
    REAL(DP) :: normalisedenergy
  END TYPE newtonsolverconvergencetest

  TYPE newton_solver_type
    TYPE(nonlinear_solver_type), POINTER :: nonlinear_solver
    INTEGER(INTG) :: newton_solve_type
    INTEGER(INTG) :: solution_initialise_type
    INTEGER(INTG) :: total_number_of_function_evaluations
    INTEGER(INTG) :: total_number_of_jacobian_evaluations
    INTEGER(INTG) :: maximum_number_of_iterations
    INTEGER(INTG) :: maximum_number_of_function_evaluations
    INTEGER(INTG) :: jacobian_calculation_type
    INTEGER(INTG) :: convergencetesttype
    REAL(DP) :: absolute_tolerance
    REAL(DP) :: relative_tolerance
    REAL(DP) :: solution_tolerance
    TYPE(newtonsolverconvergencetest), POINTER :: convergencetest
    TYPE(newton_linesearch_solver_type), POINTER :: linesearch_solver
    TYPE(newton_trustregion_solver_type), POINTER :: trustregion_solver
    TYPE(solver_type), POINTER :: linear_solver
    TYPE(solver_type), POINTER :: cellml_evaluator_solver
  END TYPE newton_solver_type

  TYPE quasi_newton_linesearch_solver_type
    TYPE(quasi_newton_solver_type), POINTER :: quasi_newton_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    INTEGER(INTG) :: linesearch_type
    REAL(DP) :: linesearch_maxstep
    REAL(DP) :: linesearch_steptolerance
    TYPE(petscmatcoloringtype) :: jacobianmatcoloring
    TYPE(petsciscoloringtype) :: jacobianiscoloring
    TYPE(petscmatfdcoloringtype) :: jacobianmatfdcoloring
    TYPE(petscsnestype) :: snes
    TYPE(petscsneslinesearchtype) :: sneslinesearch
    LOGICAL :: linesearchmonitoroutput
  END TYPE quasi_newton_linesearch_solver_type
  
  TYPE quasi_newton_trustregion_solver_type
    TYPE(quasi_newton_solver_type), POINTER :: quasi_newton_solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
    REAL(DP) :: trustregion_tolerance
    REAL(DP) :: trustregion_delta0
    TYPE(petscsnestype) :: snes
  END TYPE quasi_newton_trustregion_solver_type
  
  TYPE quasi_newton_solver_type
    TYPE(nonlinear_solver_type), POINTER :: nonlinear_solver
    INTEGER(INTG) :: quasi_newton_solve_type
    INTEGER(INTG) :: quasi_newton_type
    INTEGER(INTG) :: restart_type
    INTEGER(INTG) :: restart
    INTEGER(INTG) :: scale_type
    INTEGER(INTG) :: solution_initialise_type
    INTEGER(INTG) :: total_number_of_function_evaluations
    INTEGER(INTG) :: total_number_of_jacobian_evaluations
    INTEGER(INTG) :: maximum_number_of_iterations
    INTEGER(INTG) :: maximum_number_of_function_evaluations
    INTEGER(INTG) :: jacobian_calculation_type
    INTEGER(INTG) :: convergencetesttype
    REAL(DP) :: absolute_tolerance
    REAL(DP) :: relative_tolerance
    REAL(DP) :: solution_tolerance
    TYPE(newtonsolverconvergencetest), POINTER :: convergencetest
    TYPE(quasi_newton_linesearch_solver_type), POINTER :: linesearch_solver
    TYPE(quasi_newton_trustregion_solver_type), POINTER :: trustregion_solver
    TYPE(solver_type), POINTER :: linear_solver
    TYPE(solver_type), POINTER :: cellml_evaluator_solver
  END TYPE quasi_newton_solver_type

  TYPE nonlinear_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: nonlinear_solve_type
    TYPE(newton_solver_type), POINTER :: newton_solver
    TYPE(quasi_newton_solver_type), POINTER :: quasi_newton_solver
  END TYPE nonlinear_solver_type
  
  TYPE eigenproblem_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
  END TYPE eigenproblem_solver_type
  
  TYPE optimiser_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: solver_library
    INTEGER(INTG) :: solver_matrices_library
  END TYPE optimiser_solver_type

  TYPE cellml_evaluator_solver_type
    TYPE(solver_type), POINTER :: solver
    INTEGER(INTG) :: solver_library
    TYPE(cellml_type), POINTER :: cellml
    REAL(DP) :: current_time
  END TYPE cellml_evaluator_solver_type
  
  TYPE geometrictransformationsolvertype
    TYPE(solver_type), POINTER :: solver
    LOGICAL :: arbitrarypath
    INTEGER(INTG) :: numberofincrements
    REAL(DP), ALLOCATABLE :: scalings(:) !scaling(loadIncrementIdx), the scaling factors for each load increment, apply the full transformation in 1 load increment if unallocated. Only allocated if there are multiple load steps and if the transformation is uni-directional.
    REAL(DP), ALLOCATABLE :: transformationmatrices(:,:,:)
    TYPE(field_type), POINTER :: field
    INTEGER(INTG) :: fieldvariabletype
  END TYPE geometrictransformationsolvertype

  TYPE solver_ptr_type
    TYPE(solver_type), POINTER :: ptr
  END TYPE solver_ptr_type

  TYPE solver_type
    TYPE(solvers_type), POINTER :: solvers
    INTEGER(INTG) :: global_number
    TYPE(solver_type), POINTER :: linking_solver
    INTEGER(INTG) :: number_of_linked_solvers
    TYPE(solver_ptr_type), ALLOCATABLE :: linked_solvers(:)
    TYPE(solver_ptr_type), ALLOCATABLE :: linked_solver_type_map(:)
    LOGICAL :: solver_finished
    TYPE(varying_string) :: label
     
    INTEGER(INTG) :: output_type
    
    INTEGER(INTG) :: solve_type
    
    TYPE(linear_solver_type), POINTER :: linear_solver
    TYPE(nonlinear_solver_type), POINTER :: nonlinear_solver
    TYPE(dynamic_solver_type), POINTER :: dynamic_solver
    TYPE(dae_solver_type), POINTER :: dae_solver
    TYPE(eigenproblem_solver_type), POINTER :: eigenproblem_solver
    TYPE(optimiser_solver_type), POINTER :: optimiser_solver
    TYPE(cellml_evaluator_solver_type), POINTER :: cellml_evaluator_solver
    TYPE(geometrictransformationsolvertype), POINTER :: geometrictransformationsolver
    TYPE(solver_equations_type), POINTER :: solver_equations
    TYPE(cellml_equations_type), POINTER :: cellml_equations
    
  END TYPE solver_type

  TYPE solvers_type
    TYPE(control_loop_type), POINTER :: control_loop
    LOGICAL :: solvers_finished
    INTEGER(INTG) :: number_of_solvers
    TYPE(solver_ptr_type), ALLOCATABLE :: solvers(:)
  END TYPE solvers_type
  
  !
  !================================================================================================================================
  !
  ! Solver mapping types
  !
  
  TYPE equations_col_to_solver_cols_map_type
    INTEGER(INTG) :: number_of_solver_cols
    INTEGER(INTG), ALLOCATABLE :: solver_cols(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE equations_col_to_solver_cols_map_type

  TYPE equations_to_solver_maps_type
    INTEGER(INTG) :: equations_matrix_type
    INTEGER(INTG) :: equations_matrix_number
    INTEGER(INTG) :: solver_matrix_number
    TYPE(equations_matrix_type), POINTER :: equations_matrix
    TYPE(solver_matrix_type), POINTER :: solver_matrix
    TYPE(equations_col_to_solver_cols_map_type), ALLOCATABLE :: equations_col_to_solver_cols_map(:)
  END TYPE equations_to_solver_maps_type

  TYPE equations_to_solver_maps_ptr_type
    TYPE(equations_to_solver_maps_type), POINTER :: ptr
  END TYPE equations_to_solver_maps_ptr_type
  
  TYPE interface_to_solver_maps_type
    INTEGER(INTG) :: interface_matrix_number
    INTEGER(INTG) :: solver_matrix_number
    TYPE(interface_matrix_type), POINTER :: interface_matrix
    TYPE(solver_matrix_type), POINTER :: solver_matrix
    TYPE(equations_col_to_solver_cols_map_type), ALLOCATABLE :: interface_row_to_solver_cols_map(:)
  END TYPE interface_to_solver_maps_type

  TYPE interface_to_solver_maps_ptr_type
    TYPE(interface_to_solver_maps_type), POINTER :: ptr
  END TYPE interface_to_solver_maps_ptr_type
  
  TYPE jacobian_col_to_solver_cols_map_type
    INTEGER(INTG) :: number_of_solver_cols
    INTEGER(INTG), ALLOCATABLE :: solver_cols(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE jacobian_col_to_solver_cols_map_type

  TYPE jacobian_to_solver_map_type
    INTEGER(INTG) :: solver_matrix_number
    TYPE(equations_jacobian_type), POINTER :: jacobian_matrix
    TYPE(solver_matrix_type), POINTER :: solver_matrix
    TYPE(jacobian_col_to_solver_cols_map_type), ALLOCATABLE :: jacobian_col_to_solver_cols_map(:)
  END TYPE jacobian_to_solver_map_type

  TYPE jacobian_to_solver_map_ptr_type
    TYPE(jacobian_to_solver_map_type), POINTER :: ptr
  END TYPE jacobian_to_solver_map_ptr_type

  TYPE variable_to_solver_col_map_type
    INTEGER(INTG), ALLOCATABLE :: column_numbers(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
    REAL(DP), ALLOCATABLE :: additive_constants(:)
  END TYPE variable_to_solver_col_map_type

  TYPE equations_to_solver_matrix_maps_interface_type
    INTEGER(INTG) :: interface_condition_index
    TYPE(interface_condition_type), POINTER :: interface_condition
    INTEGER(INTG) :: interface_matrix_number
  END TYPE equations_to_solver_matrix_maps_interface_type
  
  TYPE equations_to_solver_matrix_maps_sm_type
    INTEGER(INTG) :: solver_matrix_number

    INTEGER(INTG) :: number_of_variables
    INTEGER(INTG), ALLOCATABLE :: variable_types(:)
    TYPE(field_variable_ptr_type), ALLOCATABLE :: variables(:)
    TYPE(variable_to_solver_col_map_type), ALLOCATABLE :: variable_to_solver_col_maps(:)
    INTEGER(INTG) :: number_of_dynamic_equations_matrices
    TYPE(equations_to_solver_maps_ptr_type), ALLOCATABLE :: dynamic_equations_to_solver_matrix_maps(:)
    INTEGER(INTG) :: number_of_linear_equations_matrices
    TYPE(equations_to_solver_maps_ptr_type), ALLOCATABLE :: linear_equations_to_solver_matrix_maps(:)
    INTEGER(INTG) :: number_of_equations_jacobians
    TYPE(jacobian_to_solver_map_ptr_type), ALLOCATABLE :: jacobian_to_solver_matrix_maps(:)
  END TYPE equations_to_solver_matrix_maps_sm_type

  TYPE equations_to_solver_matrix_maps_em_type
    INTEGER(INTG) :: equations_matrix_number
    INTEGER(INTG) :: number_of_solver_matrices
    TYPE(equations_to_solver_maps_ptr_type), ALLOCATABLE :: equations_to_solver_matrix_maps(:)
  END TYPE equations_to_solver_matrix_maps_em_type

  TYPE equations_row_to_solver_rows_map_type
    INTEGER(INTG) :: number_of_solver_rows
    INTEGER(INTG), ALLOCATABLE :: solver_rows(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE equations_row_to_solver_rows_map_type
  
  TYPE equations_set_to_solver_map_type
    INTEGER(INTG) :: equations_set_index
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    TYPE(equations_type), POINTER :: equations
    INTEGER(INTG) :: number_of_interface_conditions
    TYPE(equations_to_solver_matrix_maps_interface_type), ALLOCATABLE :: equations_to_solver_matrix_maps_interface(:)
    TYPE(equations_to_solver_matrix_maps_sm_type), ALLOCATABLE :: equations_to_solver_matrix_maps_sm(:)
    TYPE(equations_to_solver_matrix_maps_em_type), ALLOCATABLE :: equations_to_solver_matrix_maps_em(:)
    TYPE(jacobian_to_solver_map_ptr_type), ALLOCATABLE :: equations_to_solver_matrix_maps_jm(:)
    TYPE(equations_row_to_solver_rows_map_type), ALLOCATABLE :: equations_row_to_solver_rows_maps(:)
  END TYPE equations_set_to_solver_map_type

  TYPE interface_to_solver_matrix_maps_sm_type
    INTEGER(INTG) :: solver_matrix_number
    INTEGER(INTG) :: lagrange_variable_type
    TYPE(field_variable_type), POINTER :: lagrange_variable
    TYPE(variable_to_solver_col_map_type) :: lagrange_variable_to_solver_col_map
    INTEGER(INTG) :: number_of_dependent_variables
    INTEGER(INTG), ALLOCATABLE :: dependent_variable_types(:)
    TYPE(field_variable_ptr_type), ALLOCATABLE :: dependent_variables(:)
    TYPE(variable_to_solver_col_map_type), ALLOCATABLE :: dependent_variable_to_solver_col_maps(:)
    INTEGER(INTG) :: number_of_interface_matrices
    TYPE(interface_to_solver_maps_ptr_type), ALLOCATABLE :: interface_equations_to_solver_matrix_maps(:)
    TYPE(equations_col_to_solver_cols_map_type), ALLOCATABLE :: interface_col_to_solver_cols_map(:)
  END TYPE interface_to_solver_matrix_maps_sm_type

  TYPE interface_row_to_solver_rows_map_type
    INTEGER(INTG) :: number_of_solver_rows
    INTEGER(INTG) :: solver_row
    REAL(DP) :: coupling_coefficient
  END TYPE interface_row_to_solver_rows_map_type
  
  TYPE interface_to_solver_matrix_maps_im_type
    INTEGER(INTG) :: interface_matrix_number
    INTEGER(INTG) :: number_of_solver_matrices
    TYPE(interface_to_solver_maps_ptr_type), ALLOCATABLE :: interface_to_solver_matrix_maps(:)
    TYPE(interface_row_to_solver_rows_map_type), ALLOCATABLE :: interface_row_to_solver_rows_map(:)
  END TYPE interface_to_solver_matrix_maps_im_type

  TYPE interface_to_solver_matrix_maps_equations_type
    INTEGER(INTG) :: equations_set_index
    TYPE(equations_set_type), POINTER :: equations_set
    INTEGER(INTG) :: interface_matrix_index
  END TYPE interface_to_solver_matrix_maps_equations_type
  
  TYPE interface_column_to_solver_rows_map_type
    INTEGER(INTG) :: number_of_solver_rows
    INTEGER(INTG) :: solver_row
    REAL(DP) :: coupling_coefficient
  END TYPE interface_column_to_solver_rows_map_type
  
  TYPE interface_condition_to_solver_map_type
    INTEGER(INTG) :: interface_condition_index
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    TYPE(interface_equations_type), POINTER :: interface_equations
    INTEGER(INTG) :: number_of_equations_sets
    TYPE(interface_to_solver_matrix_maps_equations_type), ALLOCATABLE :: interface_to_solver_matrix_maps_equations(:)
    TYPE(interface_to_solver_matrix_maps_sm_type), ALLOCATABLE :: interface_to_solver_matrix_maps_sm(:)
    TYPE(interface_to_solver_matrix_maps_im_type), ALLOCATABLE :: interface_to_solver_matrix_maps_im(:)
    TYPE(interface_column_to_solver_rows_map_type), ALLOCATABLE :: interface_column_to_solver_rows_maps(:)
  END TYPE interface_condition_to_solver_map_type
  
  TYPE solver_col_to_dynamic_equations_map_type
    INTEGER(INTG) :: number_of_dynamic_equations_matrices
    INTEGER(INTG), ALLOCATABLE :: equations_matrix_numbers(:)
    INTEGER(INTG), ALLOCATABLE :: equations_col_numbers(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE solver_col_to_dynamic_equations_map_type

  TYPE solver_col_to_static_equations_map_type
    INTEGER(INTG) :: number_of_linear_equations_matrices
    INTEGER(INTG), ALLOCATABLE :: equations_matrix_numbers(:)
    INTEGER(INTG), ALLOCATABLE :: equations_col_numbers(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
    INTEGER(INTG) :: jacobian_col_number
    REAL(DP) :: jacobian_coupling_coefficient
  END TYPE solver_col_to_static_equations_map_type

  TYPE solver_dof_to_variable_map_type
    INTEGER(INTG) :: number_of_equation_dofs
    INTEGER(INTG), ALLOCATABLE :: equations_types(:)
    INTEGER(INTG), ALLOCATABLE :: equations_indices(:)
    TYPE(field_variable_ptr_type), ALLOCATABLE :: variable(:)
    INTEGER(INTG), ALLOCATABLE :: variable_dof(:)
    REAL(DP), ALLOCATABLE :: variable_coefficient(:)
    REAL(DP), ALLOCATABLE :: additive_constant(:)
  END TYPE solver_dof_to_variable_map_type
  
  TYPE solver_col_to_equations_set_map_type
    TYPE(equations_type), POINTER :: equations
    LOGICAL :: have_dynamic
    LOGICAL :: have_static
    TYPE(solver_col_to_dynamic_equations_map_type), ALLOCATABLE :: solver_col_to_dynamic_equations_maps(:)
    TYPE(solver_col_to_static_equations_map_type), ALLOCATABLE :: solver_col_to_static_equations_maps(:)
  END TYPE solver_col_to_equations_set_map_type
  
  TYPE solver_col_to_interface_equations_map_type
    INTEGER(INTG) :: number_of_interface_matrices
    INTEGER(INTG), ALLOCATABLE :: interface_matrix_numbers(:)
    INTEGER(INTG), ALLOCATABLE :: interface_col_numbers(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE solver_col_to_interface_equations_map_type

  TYPE solver_col_to_interface_map_type
    TYPE(equations_type), POINTER :: interface_equations
    TYPE(solver_col_to_interface_equations_map_type), ALLOCATABLE :: solver_col_to_interface_equations_maps(:)
  END TYPE solver_col_to_interface_map_type
  
  TYPE solver_col_to_equations_maps_type
    INTEGER(INTG) :: solver_matrix_number
    TYPE(solver_mapping_type), POINTER :: solver_mapping
    TYPE(solver_matrix_type), POINTER :: solver_matrix
    INTEGER(INTG) :: number_of_columns
    TYPE(solver_col_to_equations_set_map_type), ALLOCATABLE :: solver_col_to_equations_set_maps(:)
    TYPE(solver_col_to_interface_map_type), ALLOCATABLE :: solver_col_to_interface_maps(:)
    INTEGER(INTG) :: number_of_dofs
    INTEGER(INTG) :: total_number_of_dofs
    INTEGER(INTG) :: number_of_global_dofs
    TYPE(solver_dof_to_variable_map_type), ALLOCATABLE :: solver_dof_to_variable_maps(:)
    TYPE(domain_mapping_type), POINTER :: column_dofs_mapping
  END TYPE solver_col_to_equations_maps_type

  TYPE solver_row_to_equations_maps_type
    INTEGER(INTG) :: number_of_equations_set_rows
    INTEGER(INTG) :: interface_condition_index
    INTEGER(INTG), ALLOCATABLE :: equations_index(:)
    INTEGER(INTG), ALLOCATABLE :: rowcol_number(:)
    REAL(DP), ALLOCATABLE :: coupling_coefficients(:)
  END TYPE solver_row_to_equations_maps_type

  TYPE solver_mapping_create_values_cache_type
    TYPE(list_ptr_type), POINTER :: equations_variable_list(:)
    INTEGER, ALLOCATABLE :: dynamic_variable_type(:)
    INTEGER(INTG), ALLOCATABLE :: matrix_variable_types(:,:,:)
    INTEGER(INTG), ALLOCATABLE :: residual_variable_types(:,:)
    INTEGER(INTG), ALLOCATABLE :: rhs_variable_type(:)
    INTEGER, ALLOCATABLE :: source_variable_type(:)
    TYPE(list_ptr_type), POINTER :: interface_variable_list(:)
    TYPE(list_ptr_type), POINTER :: interface_indices(:)
  END TYPE solver_mapping_create_values_cache_type

  TYPE solver_mapping_variable_type
    TYPE(field_variable_type), POINTER :: variable
    INTEGER(INTG) :: variable_type
    INTEGER(INTG) :: number_of_equations
    INTEGER(INTG), ALLOCATABLE :: equation_types(:)
    INTEGER(INTG), ALLOCATABLE :: equation_indices(:)
  END TYPE solver_mapping_variable_type

  TYPE solver_mapping_variables_type
    INTEGER(INTG) :: number_of_variables
    TYPE(solver_mapping_variable_type), ALLOCATABLE :: variables(:)
  END TYPE solver_mapping_variables_type

  TYPE solvermappingdofcouplingstype
    INTEGER(INTG) :: numberofcouplings
    INTEGER(INTG) :: capacity
    TYPE(boundaryconditionscoupleddofsptrtype), ALLOCATABLE :: dofcouplings(:)
  END TYPE solvermappingdofcouplingstype

  TYPE solver_mapping_type
    TYPE(solver_equations_type), POINTER :: solver_equations
    LOGICAL :: solver_mapping_finished
    INTEGER(INTG) :: number_of_solver_matrices
    INTEGER(INTG) :: number_of_rows
    INTEGER(INTG) :: number_of_global_rows
    INTEGER(INTG) :: number_of_equations_sets
    TYPE(equations_set_ptr_type), ALLOCATABLE :: equations_sets(:)
    TYPE(equations_set_to_solver_map_type), ALLOCATABLE :: equations_set_to_solver_map(:)
    INTEGER(INTG) :: number_of_interface_conditions
    TYPE(interface_condition_ptr_type), ALLOCATABLE :: interface_conditions(:)
    TYPE(interface_condition_to_solver_map_type), ALLOCATABLE :: interface_condition_to_solver_map(:)
    TYPE(solver_mapping_variables_type), ALLOCATABLE :: variables_list(:)
    TYPE(solver_col_to_equations_maps_type), ALLOCATABLE :: solver_col_to_equations_cols_map(:)
    TYPE(solver_row_to_equations_maps_type), ALLOCATABLE :: solver_row_to_equations_rows_map(:)
    !LOGICAL :: HAVE_JACOBIAN !<Is .TRUE. if the Jacobian exists for nonlinear problems.
    TYPE(domain_mapping_type), POINTER :: row_dofs_mapping
    TYPE(solver_mapping_create_values_cache_type), POINTER :: create_values_cache
  END TYPE solver_mapping_type

  !
  !================================================================================================================================
  !
  ! History types
  !

  TYPE history_type
    TYPE(control_loop_type), POINTER :: control_loop
    LOGICAL :: history_finished
    INTEGER(INTG) :: file_format
    TYPE(varying_string) :: filename
    INTEGER(INTG) :: unit_number
    TYPE(field_type), POINTER :: field
  END TYPE history_type
  
  !
  !================================================================================================================================
  !
  ! Control types
  !
  
  TYPE control_loop_simple_type
    TYPE(control_loop_type), POINTER :: control_loop
  END TYPE control_loop_simple_type

  TYPE control_loop_fixed_type
    TYPE(control_loop_type), POINTER :: control_loop
    INTEGER(INTG) :: iteration_number
    INTEGER(INTG) :: start_iteration
    INTEGER(INTG) :: stop_iteration
    INTEGER(INTG) :: iteration_increment
  END TYPE control_loop_fixed_type

  TYPE control_loop_time_type
    TYPE(control_loop_type), POINTER :: control_loop
    INTEGER(INTG) :: iteration_number
    INTEGER(INTG) :: global_iteration_number
! sebk: is thei usefull?
    INTEGER(INTG) :: output_number
    INTEGER(INTG) :: input_number
    REAL(DP) :: current_time
    REAL(DP) :: start_time
    REAL(DP) :: stop_time
    REAL(DP) :: time_increment
    INTEGER(INTG) :: number_of_iterations
  END TYPE control_loop_time_type

  TYPE control_loop_while_type
    TYPE(control_loop_type), POINTER :: control_loop
    INTEGER(INTG) :: iteration_number
    INTEGER(INTG) :: maximum_number_of_iterations
    REAL(DP) :: absolute_tolerance
    LOGICAL :: continue_loop
  END TYPE control_loop_while_type

  TYPE control_loop_load_increment_type
    TYPE(control_loop_type), POINTER :: control_loop
    INTEGER(INTG) :: iteration_number
    INTEGER(INTG) :: maximum_number_of_iterations
    INTEGER(INTG) :: output_number
  END TYPE control_loop_load_increment_type

  TYPE control_loop_ptr_type
    TYPE(control_loop_type), POINTER :: ptr
  END TYPE control_loop_ptr_type

  TYPE control_loop_type
    TYPE(problem_type), POINTER :: problem
    TYPE(control_loop_type), POINTER :: parent_loop
    LOGICAL :: control_loop_finished
    TYPE(varying_string) :: label
    
    INTEGER(INTG) :: loop_type
    INTEGER(INTG) :: control_loop_level
    INTEGER(INTG) :: sub_loop_index
    INTEGER(INTG) :: output_type
    
    TYPE(control_loop_simple_type), POINTER :: simple_loop
    TYPE(control_loop_fixed_type), POINTER :: fixed_loop
    TYPE(control_loop_time_type), POINTER :: time_loop
    TYPE(control_loop_while_type), POINTER :: while_loop
    TYPE(control_loop_load_increment_type), POINTER :: load_increment_loop
    
    INTEGER(INTG) :: number_of_sub_loops
    TYPE(control_loop_ptr_type), ALLOCATABLE :: sub_loops(:)
    
    TYPE(solvers_type), POINTER :: solvers
    TYPE(history_type), POINTER :: history
  END TYPE control_loop_type  
  
  !
  !================================================================================================================================
  !
  ! Problem types
  !

  TYPE problem_setup_type
    INTEGER(INTG) :: setup_type
    INTEGER(INTG) :: action_type
  END TYPE problem_setup_type
  
  TYPE problem_type
    INTEGER(INTG) :: user_number
    INTEGER(INTG) :: global_number
    LOGICAL :: problem_finished
    TYPE(problems_type), POINTER :: problems
    INTEGER(INTG), ALLOCATABLE :: specification(:)
    TYPE(control_loop_type), POINTER :: control_loop
  END TYPE problem_type
  
  TYPE problem_ptr_type
    TYPE(problem_type), POINTER :: ptr
  END TYPE problem_ptr_type
       
  TYPE problems_type
    INTEGER(INTG) :: number_of_problems
    TYPE(problem_ptr_type), POINTER :: problems(:)
  END TYPE problems_type

  !
  !================================================================================================================================
  !
  ! Region types
  
  TYPE region_ptr_type
    TYPE(region_type), POINTER :: ptr
  END TYPE region_ptr_type
     
  TYPE region_type 
    INTEGER(INTG) :: user_number
    LOGICAL :: region_finished
    TYPE(varying_string) :: label
    TYPE(coordinate_system_type), POINTER :: coordinate_system
    TYPE(data_points_type), POINTER :: data_points
    TYPE(nodes_type), POINTER :: nodes
    TYPE(meshes_type), POINTER :: meshes
    TYPE(generated_meshes_type), POINTER :: generated_meshes
    TYPE(fields_type), POINTER :: fields
    TYPE(equations_sets_type), POINTER :: equations_sets
    TYPE(cellml_environments_type), POINTER :: cellml_environments
    TYPE(region_type), POINTER :: parent_region
    INTEGER(INTG) :: number_of_sub_regions
    TYPE(region_ptr_type), POINTER :: sub_regions(:)
    TYPE(interfaces_type), POINTER :: interfaces
  END TYPE region_type

  TYPE regions_type
    TYPE(region_type), POINTER :: world_region
  END TYPE regions_type

  !
  !================================================================================================================================
  !


END MODULE types