zhangyiss/DSurfTomo
1
0
Fork 0
mirror of https://github.com/HongjianFang/DSurfTomo.git synced 2025-05-11 12:51:19 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

Update main.f90

Browse Source
This commit is contained in:
Hongjian Fang 2018-02-03 12:09:59 -05:00 committed by GitHub
parent 6b00b85ba0
commit e239629df6
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 40 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

60
src/main.f90
View File

@ -382,28 +382,7 @@ program SurfTomo
cbst(i) = obst(i) - dsyn(i) cbst(i) = obst(i) - dsyn(i)
enddo enddo
! write out rw, iw and cbst for testing
!open(44,file='iw.dat')
!write(44,*) nar
!do i = 2,nar+1
!write(44,*) iw(i)
!enddo
!do i = 1,nar
!write(44,*) col(i)
!enddo
!close(44)
!open(44,file='rw.dat')
!do i = 1,nar
!write(44,*) rw(i)
!enddo
!close(44)
!open(44,file='residal.dat')
!do i = 1,dall
!write(44,*) cbst(i)
!enddo
!close(44)
!threshold = threshold0+(maxiter/2-iter)/3*0.5 !threshold = threshold0+(maxiter/2-iter)/3*0.5
do i = 1,dall do i = 1,dall
@ -419,6 +398,45 @@ program SurfTomo
rw(i) = rw(i)*datweight(iw(1+i)) rw(i) = rw(i)*datweight(iw(1+i))
enddo enddo
mean = sum(cbst(1:dall))/dall
std_devs = sqrt(sum(cbst(1:dall)**2)/dall - mean**2)
write(*,'(i2,a)'),iter,'th iteration...'
write(*,'(a,f7.3)'),'weight is:',weight
write(*,'(a,f8.1,a,f8.2,a,f8.3)'),'mean,std_devs and rms of &
residual: ',mean*1000,'ms ',1000*std_devs,'ms ',&
dnrm2(dall,cbst,1)/sqrt(real(dall))
write(66,'(i2,a)'),iter,'th iteration...'
write(66,'(a,f7.3)'),'weight is:',weight
write(66,'(a,f8.1,a,f8.2,a,f8.3)'),'mean,std_devs and rms of &
residual: ',mean*1000,'ms ',1000*std_devs,'ms ',&
dnrm2(dall,cbst,1)/sqrt(real(dall))
! write out rw, iw and cbst for random matrix based inversion. Use nc format latter
if (rdm == 1) then
open(44,file='iw.dat')
write(44,*) nar
do i = 2,nar+1
write(44,*) iw(i)
enddo
do i = 1,nar
write(44,*) col(i)
enddo
close(44)
open(44,file='rw.dat')
do i = 1,nar
write(44,*) rw(i)
enddo
close(44)
open(44,file='residal.dat')
do i = 1,dall
write(44,*) cbst(i)
enddo
close(44)
stop
endif
norm=0 norm=0
do i=1,nar do i=1,nar
norm(col(i))=norm(col(i))+abs(rw(i)) norm(col(i))=norm(col(i))+abs(rw(i))