fix non-F90 case

This commit is contained in:
Jeff Hammond 2023-10-11 16:05:19 +03:00 committed by edoapra
parent 5c5c34eff8
commit 8d8ad2e602
No known key found for this signature in database
GPG key ID: E6E392F4E14FB937
2 changed files with 6 additions and 6 deletions

View file

@ -358,9 +358,9 @@ C i0 ( p3 p4 h1 h2 )_vt + = 1/2 * Sum ( p5 p6 ) * t ( p5 p6 h1 h2 )_t * v (
p43d = int_mb(k_range+p3b-1)*int_mb(k_range+p4b-1)
p65d = int_mb(k_range+p5b-1)*int_mb(k_range+p6b-1)
call DGEMM('n','n',h21d,p43d,p65d,
& 0.5d0*alpha,f_a,h21d,
& f_b,p65d,
& 1.0d0, f_c,h21d)
& 0.5d0*alpha,dbl_mb(k_a),h21d,
& dbl_mb(k_b),p65d,
& 1.0d0, dbl_mb(k_c),h21d)
END IF
END IF
END IF

View file

@ -9219,9 +9219,9 @@ C i0 ( p3 p4 h1 h2 )_vt + = 1/2 * Sum ( p5 p6 ) * t ( p5 p6 h1 h2 )_t * v (
p43d = int_mb(k_range+p3b-1)*int_mb(k_range+p4b-1)
p65d = int_mb(k_range+p5b-1)*int_mb(k_range+p6b-1)
call DGEMM('n','n',h21d,p43d,p65d,
& 0.5d0*alpha,f_a,h21d,
& f_b,p65d,
& 1.0d0, f_c,h21d)
& 0.5d0*alpha,dbl_mb(k_a),h21d,
& dbl_mb(k_b),p65d,
& 1.0d0, dbl_mb(k_c),h21d)
END IF
END IF
END IF