Merge pull request #67 from gha3mi/devel

Feature Additions, Bug Fixes, and Compiler Compatibility Improvements

Author: gha3mi

README.md

@@ -108,6 +108,8 @@ To install PyVista, run the following command:
 pip install pyvista
 ```
 
+By default PyVista visualization is enabled. To disable it, define the preprocessor flag `NOSHOW_PYVISTA` in the `fpm.toml` file or pass it as a compiler flag.
+
 ### Using fpm
 
 #### Running Examples with fpm

example/nearest_point_1d.f90

@@ -7,7 +7,7 @@ program nearest_point_1d
     type(nurbs_curve) :: shape              !! Declare a NURBS curve object
     real(rk), allocatable :: Xc(:,:), Wc(:) !! Arrays for control points and weights
     real(rk) :: knot(6)                     !! Array for knot vector
-    real(rk), allocatable :: nearest_Xg(:)  !! Array for the nearest point on the curve
+    real(rk) :: nearest_Xg(3)               !! Array for the nearest point on the curve
     real(rk) :: nearest_Xt                  !! Array for the parametric coordinates of the nearest point
     integer :: id                           !! Variable for the id of the nearest point
 
@@ -70,6 +70,6 @@ program nearest_point_1d
 
     !> Finalize the NURBS curve object
     call shape%finalize()
-    deallocate(nearest_Xg, Xc, Wc)
+    deallocate(Xc, Wc)
 
 end program

example/nearest_point_2d.f90

@@ -5,8 +5,8 @@ program nearest_point_2d
     implicit none
 
     type(nurbs_surface) :: shape           !! Declare a NURBS surface object
-    real(rk), allocatable :: nearest_Xg(:) !! Coordinates of the nearest point on the surface
-    real(rk), allocatable :: nearest_Xt(:) !! Corresponding parametric coordinates of the nearest point
+    real(rk) :: nearest_Xg(3)              !! Coordinates of the nearest point on the surface
+    real(rk) :: nearest_Xt(2)              !! Corresponding parametric coordinates of the nearest point
     integer :: id                          !! id of the nearest point
     real(rk) :: Xc(4,3)                    !! Control points
     real(rk) :: Wc(4)                      !! Weights of the control points

example/nearest_point_2d_bench.f90

@@ -6,10 +6,11 @@ program nearest_point_2d_bench
     implicit none
 
     type(nurbs_surface) :: shape           !! Declare a NURBS surface object
-    real(rk), allocatable :: nearest_Xg(:) !! Coordinates of the nearest point on the surface
-    real(rk), allocatable :: nearest_Xt(:) !! Corresponding parametric coordinates of the nearest point
+    real(rk) :: nearest_Xg(3)              !! Coordinates of the nearest point on the surface
+    real(rk) :: nearest_Xt(2)              !! Corresponding parametric coordinates of the nearest point
     integer :: id                          !! id of the nearest point
     real(rk), allocatable :: points(:,:)
+    real(rk) :: pointsi(3)
     integer :: i, j
     type(timer) :: t
 
@@ -42,7 +43,8 @@ program nearest_point_2d_bench
         call random_number(points)
         call t%timer_start()
         do i = 1, size(points,1)
-            call shape%nearest_point(points(i,:), nearest_Xg, nearest_Xt, id)
+            pointsi = points(i,:)
+            call shape%nearest_point(pointsi, nearest_Xg, nearest_Xt, id)
         end do
         call t%timer_stop()
         deallocate(points)
@@ -53,6 +55,5 @@ program nearest_point_2d_bench
 
     !> Finalize the NURBS surface object
     call shape%finalize()
-    deallocate(nearest_Xg, nearest_Xt)
 
 end program

example/nearest_point_3d.f90

@@ -5,8 +5,8 @@ program nearest_point_3d
     implicit none
 
     type(nurbs_volume) :: shape            !! Declare a NURBS volume object
-    real(rk), allocatable :: nearest_Xg(:) !! Coordinates of the nearest point on the volume
-    real(rk), allocatable :: nearest_Xt(:) !! Corresponding parametric coordinates of the nearest point
+    real(rk) :: nearest_Xg(3)              !! Coordinates of the nearest point on the volume
+    real(rk) :: nearest_Xt(3)              !! Corresponding parametric coordinates of the nearest point
     integer :: id                          !! id of the nearest point
     real(rk) :: Xc(8,3)                    !! Control points
     real(rk) :: Wc(8)                      !! Weights of the control points

example/poisson_iga_solver_2d.f90

@@ -164,7 +164,7 @@ program poisson_iga_solver_2d
    end do
    call ti%timer_stop(message="L2 error evaluation                  : ")
 
-   print '(a,1pe10.4)', "L2 error norm                        = ", sqrt(l2_error)
+   print '(a,1pe11.4)', "L2 error norm                        = ", sqrt(l2_error)
    print '(a,a,a,a)', trim(filename)//".vtk", " and ", trim(filename)//"_interp.vtk", " exported"
 
    call surf%finalize()

example/poisson_iga_solver_3d.f90

@@ -166,7 +166,7 @@ program poisson_iga_solver_3d
    end do
    call ti%timer_stop(message="L2 error evaluation                  : ")
 
-   print '(a,1pe10.4)', "L2 error norm                        = ", sqrt(l2_error)
+   print '(a,1pe11.4)', "L2 error norm                        = ", sqrt(l2_error)
    print '(a,a,a,a)', trim(filename)//".vtk", " and ", trim(filename)//"_interp.vtk", " exported"
 
    call surf%finalize()