evaluate_add_spacetime_vectorByCoordinate Subroutine

public recursive subroutine evaluate_add_spacetime_vectorByCoordinate(fun, varSys, point, time, tree, nPnts, res)

This subroutine add multiple spacetime function for given coordinate points. Spacetime function are evaluated only on the coordinate which belong to st_fun shape for vectorial variable.

Arguments

Type IntentOptional Attributes Name
class(tem_varSys_op_type), intent(in) :: fun

Description of the method to obtain the variables, here some preset values might be stored, like the space time function to use or the required variables.
type(tem_varSys_type), intent(in) :: varSys

The variable system to obtain the variable from.
real(kind=rk), intent(in) :: point(:,:)

Three-dimensional coordinates at which the variable should be evaluated. Only useful for variables provided as space-time functions.
type(tem_time_type), intent(in) :: time

Point in time at which to evaluate the variable.
type(treelmesh_type), intent(in) :: tree

global treelm mesh info
integer, intent(in) :: nPnts

Number of values to obtain for this variable (vectorized access).
real(kind=rk), intent(out) :: res(:)

Resulting values for the requested variable.

Dimension: n requested entries x nComponents of this variable Access: (iElem-1)*fun%nComponents + iComp

Calls

proc~~evaluate_add_spacetime_vectorbycoordinate~~CallsGraph proc~evaluate_add_spacetime_vectorbycoordinate evaluate_add_spacetime_vectorByCoordinate interface~tem_spacetime_for tem_spacetime_for proc~evaluate_add_spacetime_vectorbycoordinate->interface~tem_spacetime_for proc~tem_coordofreal tem_CoordOfReal proc~evaluate_add_spacetime_vectorbycoordinate->proc~tem_coordofreal proc~tem_idofcoord tem_IdOfCoord proc~evaluate_add_spacetime_vectorbycoordinate->proc~tem_idofcoord proc~tem_posofid tem_PosOfId proc~evaluate_add_spacetime_vectorbycoordinate->proc~tem_posofid tem_positioninsorted tem_positioninsorted proc~evaluate_add_spacetime_vectorbycoordinate->tem_positioninsorted proc~tem_spacetime_for_coord tem_spacetime_for_coord interface~tem_spacetime_for->proc~tem_spacetime_for_coord proc~tem_spacetime_for_stcoord tem_spacetime_for_stcoord interface~tem_spacetime_for->proc~tem_spacetime_for_stcoord proc~tem_spacetime_for_treeids tem_spacetime_for_treeIDs interface~tem_spacetime_for->proc~tem_spacetime_for_treeids proc~tem_spacetime_scalar_for_index tem_spacetime_scalar_for_index interface~tem_spacetime_for->proc~tem_spacetime_scalar_for_index proc~tem_spacetime_vector_for_coord tem_spacetime_vector_for_coord interface~tem_spacetime_for->proc~tem_spacetime_vector_for_coord proc~tem_spacetime_vector_for_index tem_spacetime_vector_for_index interface~tem_spacetime_for->proc~tem_spacetime_vector_for_index proc~tem_spacetime_vector_for_treeids tem_spacetime_vector_for_treeIDs interface~tem_spacetime_for->proc~tem_spacetime_vector_for_treeids proc~tem_pathcomparison tem_PathComparison proc~tem_posofid->proc~tem_pathcomparison proc~tem_pathof tem_PathOf proc~tem_posofid->proc~tem_pathof interface~tem_spacetime_lua_for tem_spacetime_lua_for proc~tem_spacetime_for_coord->interface~tem_spacetime_lua_for interface~tem_spatial_for tem_spatial_for proc~tem_spacetime_for_coord->interface~tem_spatial_for proc~tem_abort tem_abort proc~tem_spacetime_for_coord->proc~tem_abort proc~tem_eval_acoustic_pulse tem_eval_acoustic_pulse proc~tem_spacetime_for_coord->proc~tem_eval_acoustic_pulse proc~tem_eval_cylindricalwave tem_eval_cylindricalWave proc~tem_spacetime_for_coord->proc~tem_eval_cylindricalwave proc~tem_eval_miescatter_displz tem_eval_miescatter_displz proc~tem_spacetime_for_coord->proc~tem_eval_miescatter_displz proc~tem_eval_miescatter_magnx tem_eval_miescatter_magnx proc~tem_spacetime_for_coord->proc~tem_eval_miescatter_magnx proc~tem_eval_miescatter_magny tem_eval_miescatter_magny proc~tem_spacetime_for_coord->proc~tem_eval_miescatter_magny proc~tem_polygon_material_movement_multi tem_polygon_material_movement_multi proc~tem_spacetime_for_coord->proc~tem_polygon_material_movement_multi proc~tem_polygon_material_movement_single tem_polygon_material_movement_single proc~tem_spacetime_for_coord->proc~tem_polygon_material_movement_single proc~tem_temporal_for tem_temporal_for proc~tem_spacetime_for_coord->proc~tem_temporal_for proc~tem_spacetime_for_stcoord->proc~tem_spacetime_for_coord proc~tem_spacetime_for_treeids->interface~tem_spacetime_lua_for proc~tem_spacetime_for_treeids->interface~tem_spatial_for proc~tem_spacetime_for_treeids->proc~tem_abort proc~tem_spacetime_for_treeids->proc~tem_temporal_for proc~tem_spacetime_scalar_for_index->proc~tem_spacetime_for_coord ids ids proc~tem_spacetime_scalar_for_index->ids proc~tem_spacetime_scalar_for_index->interface~tem_spacetime_lua_for proc~tem_spacetime_scalar_for_index->interface~tem_spatial_for posids posids proc~tem_spacetime_scalar_for_index->posids proc~tem_spacetime_scalar_for_index->proc~tem_temporal_for tem_precice_read tem_precice_read proc~tem_spacetime_scalar_for_index->tem_precice_read proc~tem_spacetime_vector_for_coord->interface~tem_spacetime_lua_for proc~tem_spacetime_vector_for_coord->interface~tem_spatial_for proc~tem_spacetime_vector_for_coord->proc~tem_abort proc~tem_spacetime_vector_for_coord->proc~tem_temporal_for proc~tem_spacetime_vector_for_index->proc~tem_spacetime_vector_for_coord proc~tem_spacetime_vector_for_index->ids proc~tem_spacetime_vector_for_index->interface~tem_spacetime_lua_for proc~tem_spacetime_vector_for_index->interface~tem_spatial_for proc~tem_spacetime_vector_for_index->posids proc~tem_spacetime_vector_for_index->proc~tem_temporal_for proc~tem_spacetime_vector_for_index->tem_precice_read proc~tem_spacetime_vector_for_treeids->interface~tem_spacetime_lua_for proc~tem_spacetime_vector_for_treeids->interface~tem_spatial_for proc~tem_spacetime_vector_for_treeids->proc~tem_abort proc~tem_spacetime_vector_for_treeids->proc~tem_temporal_for

Source Code

  recursive subroutine evaluate_add_spacetime_vectorByCoordinate( fun, varsys, &
    & point, time, tree, nPnts, res                                            )
    !--------------------------------------------------------------------------!
    !> Description of the method to obtain the variables, here some preset
    !! values might be stored, like the space time function to use or the
    !! required variables.
    class(tem_varSys_op_type), intent(in) :: fun

    !> The variable system to obtain the variable from.
    type(tem_varSys_type), intent(in) :: varSys

    !> Three-dimensional coordinates at which the variable should be
    !! evaluated. Only useful for variables provided as space-time functions.
    real(kind=rk), intent(in) :: point(:,:)

    !> Point in time at which to evaluate the variable.
    type(tem_time_type), intent(in)  :: time

    !> global treelm mesh info
    type(treelmesh_type), intent(in) :: tree

    !> Number of values to obtain for this variable (vectorized access).
    integer, intent(in) :: nPnts

    !> Resulting values for the requested variable.
    !!
    !! Dimension: n requested entries x nComponents of this variable
    !! Access: (iElem-1)*fun%nComponents + iComp
    real(kind=rk), intent(out) :: res(:)
    !--------------------------------------------------------------------------!
    ! --------------------------------------------------------------------------
    integer :: iStFun, iPoint, elemPos, pos
    type(tem_st_fun_listElem_type), pointer :: fPtr
    integer :: coord(4)
    integer(kind=long_k):: treeID
    real(kind=rk) :: st_res_pnt(1, fun%nComponents)
    real(kind=rk) :: st_res(nPnts, fun%nComponents)
    ! --------------------------------------------------------------------------
    !call C_F_POINTER( fun%method_Data, fPtr )
    ! avoid compiler warning about unused varsys
    call C_F_POINTER( varsys%method%val(fun%mypos)%method_Data, fPtr )

    st_res = 0.0_rk

    ! The space time function has a custom shape, thus we have to check
    ! whether the coordinates are covered by the space time function.
    ! As we are dealing with vectorial access here, we have to loop over all
    ! requested points.
    stLoop: do iStFun = 1, fPtr%nVals
      if (fPtr%val(iStFun)%subTree%useGlobalMesh) then
        st_res = st_res + tem_spacetime_for( me      = fPtr%val(iStFun), &
          &                                  coord   = point,            &
          &                                  time    = time,             &
          &                                  n       = nPnts,            &
          &                                  nComp   = fun%nComponents   )
      else
        ! The space time function has a custom shape, thus we have to check
        ! whether the coordinates are covered by the space time function.
        do iPoint = 1, nPnts
          coord =  tem_CoordOfReal(tree, point(iPoint,:), tree%global%maxLevel)
          treeId = tem_IdOfCoord(coord)
          elemPos = tem_PosofId(treeId, tree%treeID)
          pos = tem_PositionInSorted(                      &
            & me    = fPtr%val(iStFun)%subTree%map2global, &
            & val   = elemPos                              )
          ! When the element is covered by the space time function's shape, we
          ! evaluate it at the given coordinates.
          if (pos > 0) then
            st_res_pnt = tem_spacetime_for(                        &
              & me      = fPtr%val(iStFun),                        &
              & coord   = reshape( source = (/ point(iPoint,1),    &
              &                                point(iPoint,2),    &
              &                                point(iPoint,3) /), &
              &                    shape  = (/ 1, 3 /) ),          &
              & time    = time,                                    &
              & n       = 1,                                       &
              & nComp   = fun%nComponents                          )

            st_res(iPoint,:) = st_res(iPoint,:) +  st_res_pnt(1,:)
          end if ! point part of subtree
        end do !iPoint
      end if ! global mesh
    end do stLoop

    do iPoint = 1, nPnts
      res( (iPoint-1)*fun%nComponents + 1 : iPoint*fun%nComponents ) &
        & = st_res(iPoint, :)
    end do

  end subroutine evaluate_add_spacetime_vectorByCoordinate