test: stdlib + external blas

Author: perazz

example_packages/metapackage_stdlib_extblas/app/main.f90

@@ -0,0 +1,22 @@
+! fortran-lang stdlib test case
+! This test program will only run if stdlib is properly built and linked to this project. 
+program test_stdlib_metapackage
+
+    ! These USEs would not be possible if stdlib is not found 
+    use metapackage_stdlib, only: external_blas_test,external_lapack_test
+    implicit none
+
+    logical :: ext_blas,ext_lapack
+
+    call external_blas_test(ext_blas)
+    call external_lapack_test(ext_lapack)
+    
+    if (.not.ext_blas) then 
+        stop 1
+    elseif (.not.ext_lapack) then 
+        stop 2
+    else
+        stop 0
+    end if
+
+end program test_stdlib_metapackage

example_packages/metapackage_stdlib_extblas/fpm.toml

@@ -0,0 +1,3 @@
+name = "test_stdlib_ext_blas"
+dependencies.blas = "*"
+dependencies.stdlib = "*" 

example_packages/metapackage_stdlib_extblas/src/metapackage_stdlib.f90

@@ -0,0 +1,81 @@
+module metapackage_stdlib
+    use stdlib_linalg_constants, only: sp,dp,ilp
+    implicit none
+    private
+    
+    public :: external_blas_test
+    public :: external_lapack_test
+
+    contains
+  
+    !> Test availability of the external BLAS interface
+    subroutine external_blas_test(external_blas)
+        !> Error handling
+        logical, intent(out) :: external_blas
+
+#ifdef STDLIB_EXTERNAL_BLAS           
+        interface 
+            subroutine saxpy(n,sa,sx,incx,sy,incy)
+                 import sp,ilp 
+                 implicit none(type,external) 
+                 real(sp), intent(in) :: sa,sx(*)
+                 integer(ilp), intent(in) :: incx,incy,n
+                 real(sp), intent(inout) :: sy(*)
+            end subroutine saxpy        
+        end interface
+        
+        integer(ilp), parameter :: n = 5, inc=1
+        real(sp) :: a,x(n),y(n)
+             
+        x = 1.0_sp
+        y = 2.0_sp
+        a = 3.0_sp
+        
+        call saxpy(n,a,x,inc,y,inc)
+        
+        ! Result must also be correct
+        external_blas = all(abs(y-5.0_sp)<sqrt(epsilon(0.0_sp)))
+        
+#else   
+        external_blas = .false.
+#endif        
+        
+    end subroutine external_blas_test        
+                
+    !> Test availability of the external BLAS interface
+    subroutine external_lapack_test(external_lapack)
+        !> Error handling
+        logical, intent(out) :: external_lapack
+
+#ifdef STDLIB_EXTERNAL_LAPACK       
+        interface 
+           subroutine dgetrf( m, n, a, lda, ipiv, info )
+                import dp,ilp
+                implicit none(type,external) 
+                integer(ilp), intent(out) :: info,ipiv(*)
+                integer(ilp), intent(in) :: lda,m,n
+                real(dp), intent(inout) :: a(lda,*)
+           end subroutine dgetrf       
+        end interface
+        
+        integer(ilp), parameter :: n = 3
+        real(dp) :: A(n,n)
+        integer(ilp) :: ipiv(n),info
+
+
+        A = reshape([1,0,0, 0,1,0, 0,0,1],[3,3])
+        info = 123
+
+        ! Factorize matrix 
+        call dgetrf(n,n,A,n,ipiv,info)
+
+        ! Result must be correct
+        external_lapack = info==0
+        
+#else
+        external_lapack = .false.
+#endif        
+        
+    end subroutine external_lapack_test  
+  
+end module metapackage_stdlib

example_packages/preprocess_cpp/app/main.f90

@@ -0,0 +1,4 @@
+program cpp
+    use preprocess_cpp
+    call say_hello() 
+end program