Determinant_Tensor2D_t Subroutine

  1. SELF_Tensor_2D_t.f90
  2. SELF_Tensor_2D_t
  3. Determinant_Tensor2D_t

public subroutine Determinant_Tensor2D_t(this, det)

Arguments

TypeIntentOptionalAttributesName
class(Tensor2D_t), intent(in) :: this
real(kind=prec), intent(out) :: det(1:this%N+1,1:this%N+1,1:this%nelem,1:this%nvar)

Contents

Source Code

Determinant_Tensor2D_t

Source Code

  subroutine Determinant_Tensor2D_t(this,det)
    implicit none
    class(Tensor2D_t),intent(in) :: this
    real(prec),intent(out) :: det(1:this%N+1,1:this%N+1,1:this%nelem,1:this%nvar)
    ! Local
    integer :: iEl,iVar,i,j

    do concurrent(i=1:this%N+1,j=1:this%N+1, &
                  iel=1:this%nelem,ivar=1:this%nvar)

      det(i,j,iEl,iVar) = this%interior(i,j,iEl,iVar,1,1)* &
                          this%interior(i,j,iEl,iVar,2,2)- &
                          this%interior(i,j,iEl,iVar,1,2)* &
                          this%interior(i,j,iEl,iVar,2,1)

    enddo

  endsubroutine Determinant_Tensor2D_t