updating lobatto poly

src/modules/Polynomial/src/LobattoPolynomialUtility.F90

@@ -30,6 +30,7 @@ MODULE LobattoPolynomialUtility
 PUBLIC :: LobattoEval
 PUBLIC :: LobattoEvalAll
 PUBLIC :: LobattoKernelEvalAll
+PUBLIC :: LobattoKernelEvalAll_
 PUBLIC :: LobattoKernelGradientEvalAll
 PUBLIC :: LobattoMonomialExpansionAll
 PUBLIC :: LobattoMonomialExpansion
@@ -244,6 +245,34 @@ MODULE PURE FUNCTION LobattoKernelEvalAll1(n, x) RESULT(ans)
   END FUNCTION LobattoKernelEvalAll1
 END INTERFACE LobattoKernelEvalAll
 
+!----------------------------------------------------------------------------
+!                                                      LobattoKernelEvalAll
+!----------------------------------------------------------------------------
+
+!> author: Vikas Sharma, Ph. D.
+! date: 6 Sept 2022
+! summary: Evaluate Lobatto bubble functions order = 0 to n at several points
+!
+!# Introduction
+!
+! Evaluate Lobatto bubble polynomials from order = 0 to n at several points
+!
+!- N, the highest order polynomial to compute.
+!- x: the point at which the polynomials are to be evaluated.
+
+INTERFACE LobattoKernelEvalAll_
+  MODULE PURE SUBROUTINE LobattoKernelEvalAll1_(n, x, ans, nrow, ncol)
+    INTEGER(I4B), INTENT(IN) :: n
+    !! n
+    REAL(DFP), INTENT(IN) :: x(:)
+    REAL(DFP), INTENT(INOUT) :: ans(1:, 0:)
+    !! ans(1:SIZE(x), 0:n)
+    !! Evaluate Lobatto polynomial of order = 0 to n (total n+1)
+    !! at point x
+    INTEGER(I4B), INTENT(OUT) :: nrow, ncol
+  END SUBROUTINE LobattoKernelEvalAll1_
+END INTERFACE LobattoKernelEvalAll_
+
 !----------------------------------------------------------------------------
 !                                               LobattoKernelGradientEvalAll
 !----------------------------------------------------------------------------

src/submodules/Polynomial/src/LobattoPolynomialUtility@Methods.F90

@@ -170,16 +170,27 @@
 !----------------------------------------------------------------------------
 
 MODULE PROCEDURE LobattoKernelEvalAll1
+INTEGER(I4B) :: nrow, ncol
+CALL LobattoKernelEvalAll1_(n=n, x=x, ans=ans, nrow=nrow, ncol=ncol)
+END PROCEDURE LobattoKernelEvalAll1
+
+!----------------------------------------------------------------------------
+!                                                       LobattoKernelEvalAll
+!----------------------------------------------------------------------------
+
+MODULE PROCEDURE LobattoKernelEvalAll1_
 REAL(DFP) :: m, avar
 INTEGER(I4B) :: ii
 
-ans = UltrasphericalEvalAll(n=n, x=x, lambda=1.5_DFP)
+CALL UltrasphericalEvalAll_(n=n, x=x, lambda=1.5_DFP, ans=ans, nrow=nrow, &
+                            ncol=ncol)
+
 DO ii = 0, n
   m = REAL(ii, KIND=DFP)
   avar = -SQRT(8.0_DFP*(2.0_DFP*m+3.0_DFP))/(m+1.0_DFP)/(m+2.0_DFP)
-  ans(:, ii) = avar * ans(:, ii)
+  ans(1:nrow, ii) = avar * ans(1:nrow, ii)
 END DO
-END PROCEDURE LobattoKernelEvalAll1
+END PROCEDURE LobattoKernelEvalAll1_
 
 !----------------------------------------------------------------------------
 !                                             LobattoKernelGradientEvalAll1
@@ -401,22 +412,21 @@
 
 MODULE PROCEDURE LobattoStiffnessMatrix
 INTEGER(I4B) :: ii
-REAL(DFP) :: m
-!!
+
 ans = 0.0_DFP
-!!
+
 DO ii = 1, n + 1
   ans(ii, ii) = 1.0_DFP
 END DO
-!!
+
 ans(1, 1) = 0.5_DFP
-!!
+
 IF (n .EQ. 0_I4B) RETURN
-!!
+
 ans(2, 2) = 0.5_DFP
 ans(1, 2) = -0.5_DFP
 ans(2, 1) = ans(1, 2)
-!!
+
 END PROCEDURE LobattoStiffnessMatrix
 
 !----------------------------------------------------------------------------