Created function submodule and interface.

Author: ghbrown

src/stdlib_linalg.fypp

@@ -11,6 +11,7 @@ module stdlib_linalg
   public :: diag
   public :: eye
   public :: trace
+  public :: outer_product
 
   interface diag
     !! version: experimental
@@ -52,6 +53,7 @@ module stdlib_linalg
     #:endfor
   end interface
 
+
   ! Matrix trace
   interface trace
     !! version: experimental
@@ -63,6 +65,22 @@ module stdlib_linalg
     #:endfor
   end interface
 
+
+  ! Outer product (of two vectors)
+  interface outer_product
+    !! version: experimental
+    !!
+    !! 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}$
+    #:endfor
+  end interface outer_product
+
 contains
 
     function eye(n) result(res)

src/stdlib_linalg_outer_product.fypp

@@ -0,0 +1,20 @@
+#:include "common.fypp"
+#:set RCI_KINDS_TYPES = REAL_KINDS_TYPES + CMPLX_KINDS_TYPES + INT_KINDS_TYPES
+submodule (stdlib_linalg) stdlib_linalg_outer_product
+
+  implicit none
+
+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
+
+end submodule