This function invokes the Lua function, which takes barycentric coordinates of an element.
Lua function defined in the script is connected to the conf handle and return the result of the function. The Lua function takes barycentric coordinate as input argument i.e fun_name(x,y,z)
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | fun_ref |
Lua reference to the function to evaluate. |
||
type(flu_State) | :: | conf |
lua state |
|||
real(kind=rk), | intent(in) | :: | coord(n,3) |
barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates |
||
integer, | intent(in) | :: | n |
number of return values |
return value
function tem_spatial_lua_for_coord( fun_ref, conf, coord, n ) result(res)
! -------------------------------------------------------------------- !
!> Lua reference to the function to evaluate.
integer, intent(in) :: fun_ref
!> number of return values
integer, intent(in) :: n
!> barycentric Ids of an elements.
!! 1st index goes over number of elements and
!! 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: coord(n,3)
!> lua state
type(flu_State) :: conf
!> return value
real(kind=rk) :: res(n)
! -------------------------------------------------------------------- !
type(aot_fun_type) :: fun
integer :: iError
integer :: iDir, jDir
! -------------------------------------------------------------------- !
call aot_fun_open(L=conf, fun=fun, ref=fun_ref)
do iDir = 1, n
do jDir = 1, 3
call aot_fun_put( L=conf, fun=fun, arg=coord(iDir,jDir) )
end do
call aot_fun_do(L=conf, fun=fun, nresults=1)
call aot_top_get_val(L=conf, val=res(iDir), ErrCode=iError)
if ( btest(iError,aoterr_Fatal) ) then
write(logunit(0),*) "ERROR Obtaining a spatial function"
write(logunit(0),*) "Probably wrong number of components returned"
write(logunit(0),*) "or a scalar was return as a lua table"
write(logunit(0),*) 'Expected nComp: 1'
write(logunit(0),*) 'ErrorCodes: ', iError
write(logunit(0),*) "Check return values of your Lua functions!"
call tem_abort()
end if
end do
call aot_fun_close(L=conf, fun=fun)
end function tem_spatial_lua_for_coord