Commit 329abd33 authored by Pavel Kus's avatar Pavel Kus

included beginning of the q array to mult

parent 847b8b81
......@@ -91,7 +91,7 @@
real(kind=C_DATATYPE_KIND) :: w_comb(nq, 4)
real(kind=C_DATATYPE_KIND) :: h_comb(4)
real(kind=C_DATATYPE_KIND) :: h_mat(4, nb-4)
real(kind=C_DATATYPE_KIND) :: h_mat(4, nb-1)
real(kind=C_DATATYPE_KIND) :: tau1, tau2, tau3, tau4
#endif /* REALCASE==1 */
......@@ -271,26 +271,24 @@
w(1:nq) = w(1:nq) * h1 - ( z(1:nq) * h4 + y(1:nq) * h3 + x(1:nq) * h2)
q(1:nq,1) = q(1:nq,1) - w(1:nq)
h4 = hh(2,4)
q(1:nq,2) = q(1:nq,2) - (w(1:nq) * h4 + z(1:nq))
h3 = hh(2,3)
h4 = hh(3,4)
q(1:nq,2) = q(1:nq,2) - z(1:nq)
q(1:nq,3) = q(1:nq,3) - y(1:nq)
q(1:nq,3) = -( z(1:nq) * h3) + q(1:nq,3)
q(1:nq,3) = -( w(1:nq) * h4) + q(1:nq,3)
h2 = hh(2,2)
h3 = hh(3,3)
h4 = hh(4,4)
q(1:nq,4) = q(1:nq,4) - x(1:nq)
q(1:nq,4) = -(y(1:nq) * h2) + q(1:nq,4)
q(1:nq,4) = -(z(1:nq) * h3) + q(1:nq,4)
q(1:nq,4) = -(w(1:nq) * h4) + q(1:nq,4)
! h4 = hh(2,4)
! q(1:nq,2) = q(1:nq,2) - (w(1:nq) * h4)
!
! h3 = hh(2,3)
! h4 = hh(3,4)
! q(1:nq,3) = -( z(1:nq) * h3) + q(1:nq,3)
! q(1:nq,3) = -( w(1:nq) * h4) + q(1:nq,3)
!
! h2 = hh(2,2)
! h3 = hh(3,3)
! h4 = hh(4,4)
! q(1:nq,4) = -(y(1:nq) * h2) + q(1:nq,4)
! q(1:nq,4) = -(z(1:nq) * h3) + q(1:nq,4)
! q(1:nq,4) = -(w(1:nq) * h4) + q(1:nq,4)
w_comb(:,1) = x
w_comb(:,2) = y
......@@ -308,12 +306,12 @@
!
! enddo
!h_mat(:,:) = 0.0
h_mat(1,1:nb-4) = -hh(2:nb-3, 1)
h_mat(2,1:nb-4) = -hh(3:nb-2, 2)
h_mat(3,1:nb-4) = -hh(4:nb-1, 3)
h_mat(4,1:nb-4) = -hh(5:nb, 4)
q(1:nq, 5:nb) = matmul(w_comb, h_mat) + q(1:nq, 5:nb)
h_mat(:,:) = 0.0
h_mat(1,4:nb-1) = -hh(2:nb-3, 1)
h_mat(2,3:nb-1) = -hh(2:nb-2, 2)
h_mat(3,2:nb-1) = -hh(2:nb-1, 3)
h_mat(4,1:nb-1) = -hh(2:nb, 4)
q(1:nq, 2:nb) = matmul(w_comb, h_mat) + q(1:nq, 2:nb)
......
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