linearalgebrautilsmodule Module Reference

Functions/Subroutines
pure real(dp) function, dimension(:, :), allocatable, public	eye (n)
pure real(dp) function, dimension(size(a), size(b)), public	outer_product (A, B)
real(dp) function, dimension(3), public	cross_product (a, b)

Function/Subroutine Documentation

cross_product()

real(dp) function, dimension(3), public linearalgebrautilsmodule::cross_product	(	real(dp), dimension(3), intent(in)	a,
		real(dp), dimension(3), intent(in)	b
	)

Definition at line 39 of file LinearAlgebraUtils.f90.

    ! -- return
    real(DP), dimension(3) :: c
    ! -- dummy
    real(DP), dimension(3), intent(in) :: a
    real(DP), dimension(3), intent(in) :: b
 
    c(1) = a(2) * b(3) - a(3) * b(2)
    c(2) = a(3) * b(1) - a(1) * b(3)
    c(3) = a(1) * b(2) - a(2) * b(1)

eye()

pure real(dp) function, dimension(:, :), allocatable, public linearalgebrautilsmodule::eye

(

integer(i4b), intent(in)

n

)

Definition at line 11 of file LinearAlgebraUtils.f90.

    ! -- dummy
    integer(I4B), intent(in) :: n
    real(DP), dimension(:, :), allocatable :: A
    ! -- locals
    integer(I4B) :: i
 
    allocate (a(n, n))
    a = 0.0_dp
    do i = 1, n
      a(i, i) = 1.0_dp
    end do

Here is the caller graph for this function:

outer_product()

pure real(dp) function, dimension(size(a), size(b)), public linearalgebrautilsmodule::outer_product	(	real(dp), dimension(:), intent(in)	A,
		real(dp), dimension(:), intent(in)	B
	)

Definition at line 25 of file LinearAlgebraUtils.f90.

    ! -- dummy
    real(DP), intent(in) :: A(:), B(:)
    real(DP) :: AB(size(A), size(B))
    ! -- locals
    integer :: nA, nB
 
    na = size(a)
    nb = size(b)
 
    ab = spread(source=a, dim=2, ncopies=nb) * &
         spread(source=b, dim=1, ncopies=na)

Here is the caller graph for this function: