Scheduled maintenance on Monday 2019-06-24 between 10:00-11:00 CEST

Commit 08265889 authored by Pavel Kus's avatar Pavel Kus

changing the mat-vec multiplication to ommit the diagonal entry

parent 4604501f
......@@ -108,11 +108,6 @@
-1.0_rk8, h_mat, 4, &
0.0_rk8, s_mat, 4)
s_mat(1,1) = 1
s_mat(2,2) = 1
s_mat(3,3) = 1
s_mat(4,4) = 1
! w_comb = - matmul(q(1:nq, 1:nb+3), transpose(h_mat))
call DGEMM('N', 'T', nq, 4, nb+3, &
-1.0_rk8, q, ldq, &
......@@ -121,7 +116,7 @@
! Rank-1 update
do i = 1, 4
w_comb(1:nq,i) = matmul(w_comb(1:nq,1:i), hh(1,i) * s_mat(i,1:i))
w_comb(1:nq,i) = matmul(w_comb(1:nq,1:i-1), hh(1,i) * s_mat(i,1:i-1)) + hh(1,i) * w_comb(1:nq, i)
enddo
!q(1:nq, 1:nb+3) = matmul(w_comb, h_mat) + q(1:nq, 1:nb+3)
......
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