Feel free to join the next Helmholtz Hacky Hour #30 on Wednesday, June 16, 2021 from 2PM to 3PM!

Commit c0870758 authored by Stephan Thober's avatar Stephan Thober
Browse files

deactivated timer, set solver to dgesv in OK.f90

parent e16afb5e
......@@ -76,7 +76,7 @@ subroutine OK(jd,k)
B(nNmax+1) = 1.0_dp
!NOTE: only the upper triangular matrix is needed!
! call D_LSASF (A, B, X)
call gesv(A, B)
call dgesv(A, B)
X = B
!
! The BLUE of z is then:
......
......@@ -12,15 +12,15 @@ subroutine Timer
use mo_kind, only : i4
implicit none
real(i4) :: DTIME
if (RunTime <= 0.0) then
RunTime = DTIME(TA)
else
open (unit=500, file='timer.out', status='unknown')
RunTime = DTIME(TA)
write (500, 100) RunTime, TA(1), TA(2)
close (500, status='keep')
end if
! if (RunTime <= 0.0) then
! RunTime = DTIME(TA)
! else
! open (unit=500, file='timer.out', status='unknown')
! RunTime = DTIME(TA)
! write (500, 100) RunTime, TA(1), TA(2)
! close (500, status='keep')
! end if
print *, '***WARNING: timing not updated'
! formats
100 format ( '-------------------------------------',&
......
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