Telluride Matrix-Vector Multiplication Routine for Orthogonal Operator

SUBROUTINE Y_EQ_AX_ORTHOGONAL (nunk, X, Y, status)

  [snip]

  type(JT_vector_type), intent(in) :: X

  [snip]

  ! Gather element face neighbors
  call EE_GATHER (X_Neighbors, X%Values, Mesh, &
                  BOUNDARY=X%Aux1)

  ! Perform the matrix-vector multiply
  Y = Ortho_Coeff(0,:) * X%Values
  do f = 1,nfc
     Y = Y + Ortho_Coeff(f,:) * X_Neighbors(f,:)
  end do

  return
END SUBROUTINE Y_EQ_AX_ORTHOGONAL