Commit f1082e5e authored by Luca's avatar Luca

minor

parent cebe2dba
...@@ -30,12 +30,6 @@ CONTAINS ...@@ -30,12 +30,6 @@ CONTAINS
complex(kind=prec) :: res(product((/(i,i=1,size(v1)+size(v2))/))/ & complex(kind=prec) :: res(product((/(i,i=1,size(v1)+size(v2))/))/ &
(product((/(i,i=1,size(v1))/))*product((/(i,i=1,size(v2))/))), & (product((/(i,i=1,size(v1))/))*product((/(i,i=1,size(v2))/))), &
size(v1) + size(v2)) size(v1) + size(v2))
complex(kind=prec) :: p1(product((/(i,i=1,size(v1)+size(v2)-1)/))/ &
(product((/(i,i=1,size(v1)-1)/))*product((/(i,i=1,size(v2))/))), &
size(v1) + size(v2) - 1)
complex(kind=prec) :: p2(product((/(i,i=1,size(v1)+size(v2)-1)/))/ &
(product((/(i,i=1,size(v1))/))*product((/(i,i=1,size(v2)-1)/))), &
size(v1) + size(v2) - 1)
complex(kind=prec) :: alpha, beta, w1(size(v1)-1), w2(size(v2)-1) complex(kind=prec) :: alpha, beta, w1(size(v1)-1), w2(size(v2)-1)
res = 0 res = 0
......
...@@ -102,7 +102,7 @@ CONTAINS ...@@ -102,7 +102,7 @@ CONTAINS
end subroutine do_GPL_tests end subroutine do_GPL_tests
subroutine do_shuffle_tests() subroutine do_shuffle_tests()
complex(kind=prec) :: v(3) = cmplx((/1,2/)) complex(kind=prec) :: v(2) = cmplx((/1,2/))
complex(kind=prec) :: w(2) = cmplx((/3,4/)) complex(kind=prec) :: w(2) = cmplx((/3,4/))
call print_matrix(shuffle_product(v,w)) call print_matrix(shuffle_product(v,w))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment