zhangyiss/DSurfTomo
1
0
Fork 0
mirror of https://github.com/HongjianFang/DSurfTomo.git synced 2025-05-07 15:41:14 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge e239629df6 into 6b943f362a

Browse Source
This commit is contained in:
Hongjian Fang 2018-07-04 17:45:28 +00:00 committed by GitHub
commit 8dad104667
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

62
src/main.f90
View File

@ -382,28 +382,7 @@ program SurfTomo
cbst(i) = obst(i) - dsyn(i)
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
do i = 1,dall
@ -418,6 +397,45 @@ program SurfTomo
do i = 1,nar
rw(i) = rw(i)*datweight(iw(1+i))
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
do i=1,nar
@ -534,7 +552,7 @@ program SurfTomo
do i =1,dall
cbst(i)=cbst(i)/datweight(i)
enddo
mean = sum(cbst(1:dall))/dall
std_devs = sqrt(sum(cbst(1:dall)**2)/dall - mean**2)
write(*,'(i2,a)'),iter,'th iteration...'