cholesky: allocatable -> assumed-shape return value

Author: perazz

src/stdlib_linalg.fypp

@@ -74,7 +74,7 @@ module stdlib_linalg
          !> [optional] should the unused half of the return matrix be zeroed out? Default: yes
          logical(lk), optional, intent(in) :: other_zeroed
          !> Output matrix with Cholesky factors c[n,n]
-         ${rt}$, allocatable :: c(:,:)                  
+         ${rt}$ :: c(size(a,1),size(a,2))                  
      end function stdlib_linalg_${ri}$_cholesky_fun
      #:endfor    
   end interface chol  

src/stdlib_linalg_cholesky.fypp

@@ -158,7 +158,7 @@ submodule (stdlib_linalg) stdlib_linalg_cholesky
          !> [optional] should the unused half of the return matrix be zeroed out? Default: yes
          logical(lk), optional, intent(in) :: other_zeroed
 
-         allocate(c,source=a)
+         c = a
 
          call stdlib_linalg_${ri}$_cholesky_inplace(c,lower,other_zeroed)