Outer product

src/Makefile.manual

@@ -6,7 +6,7 @@ SRCFYPP =\
  stdlib_io.fypp \
  stdlib_linalg.fypp \
  stdlib_linalg_diag.fypp \
-stdlib_linalg_outer_product.fypp \
+ stdlib_linalg_outer_product.fypp \
  stdlib_optval.fypp \
  stdlib_quadrature.fypp \
  stdlib_quadrature_trapz.fypp \
@@ -129,3 +129,4 @@ stdlib_stats_distribution_PRNG.o: \
 stdlib_string_type.o: stdlib_ascii.o stdlib_kinds.o
 stdlib_strings.o: stdlib_ascii.o stdlib_string_type.o
 stdlib_math.o: stdlib_kinds.o
+stdlib_linalg_outer_product.o: stdlib_linalg.o

src/stdlib_linalg.fypp

@@ -73,11 +73,10 @@ module stdlib_linalg
  !! Computes outer product of two vectors, returning a matrix
  !! ([Specification](../page/specs/stdlib_linalg.html#description_3))
  #:for k1, t1 in RCI_KINDS_TYPES
- module function outer_product_${t1[0]}$${k1}$(u,v) result(res)
- ${t1}$, intent(in) :: u(:), v(:)
- ${t1}$ :: res(size(u),size(v))
- integer :: col
- end function outer_product_${t1[0]}$${k1}$
+ module function outer_product_${t1[0]}$${k1}$(u, v) result(res)
+ ${t1}$, intent(in) :: u(:), v(:)
+ ${t1}$ :: res(size(u),size(v))
+ end function outer_product_${t1[0]}$${k1}$
  #:endfor
  end interface outer_product
 

src/stdlib_linalg_outer_product.fypp

@@ -6,15 +6,15 @@ submodule (stdlib_linalg) stdlib_linalg_outer_product
 
 contains
 
-  #:for k1, t1 in RCI_KINDS_TYPES
-  module function outer_product_${t1[0]}$${k1}$(u,v) result(res)
-  ${t1}$, intent(in) :: u(:), v(:)
-  ${t1}$ :: res(size(u),size(v))
-  integer :: col
-  do col=1,size(v)
-  res(:,col) = v(col) * u
-  end do
-  end function outer_product_${t1[0]}$${k1}$
-  #:endfor
+ #:for k1, t1 in RCI_KINDS_TYPES
+ module function outer_product_${t1[0]}$${k1}$(u, v) result(res)
+ ${t1}$, intent(in) :: u(:), v(:)
+ ${t1}$ :: res(size(u),size(v))
+ integer :: col
+ do col = 1, size(v)
+ res(:,col) = v(col) * u
+ end do
+ end function outer_product_${t1[0]}$${k1}$
+ #:endfor
 
 end submodule