Updates in LineInterpolationUtility

- fixing minor issues
- setting proper number of gauss points for Radau type quadratures
- sorting with explicit range for "INC" case in getting interpolation
points

Author: shishiousan

src/submodules/Line/src/[email protected]

@@ -391,7 +391,7 @@ SUBROUTINE handle_increasing
   !! layout INCREASING
   IF (layout(1:1) .EQ. "I") THEN
     DO ii = 1, nrow
-      CALL HeapSort(ans(ii, :))
+      CALL HeapSort(ans(ii, 1:ncol))
     END DO
   END IF
 END SUBROUTINE

src/submodules/Line/src/[email protected]

@@ -42,6 +42,11 @@
 CASE (qpopt%GaussLegendre, qpopt%GaussChebyshev, &
       qpopt%GaussJacobi, qpopt%GaussUltraspherical)
   ans = 1_I4B + INT(order / 2, kind=I4B)
+CASE (qpopt%GaussLegendreRadauRight, qpopt%GaussLegendreRadauLeft, &
+      qpopt%GaussChebyshevRadauLeft, qpopt%GaussChebyshevRadauRight, &
+      qpopt%GaussJacobiRadauLeft, qpopt%GaussJacobiRadauRight, &
+      qpopt%GaussUltraSphericalRadauLeft, qpopt%GaussUltraSphericalRadauRight)
+  ans = 2_I4B + INT((order - 1) / 2, kind=I4B)
 CASE DEFAULT
   ans = 2_I4B + INT(order / 2, kind=I4B)
 END SELECT