[cig-commits] r22688 - seismo/3D/ADJOINT_TOMO/ADJOINT_TOMOGRAPHY_TOOLKIT/ITERATION_UPDATE/X04_SRC_DIRECTION_LBFGS

[hejunzhu at geodynamics.org]

Author: hejunzhu
Date: 2013-07-30 13:07:39 -0700 (Tue, 30 Jul 2013)
New Revision: 22688

Modified:
   seismo/3D/ADJOINT_TOMO/ADJOINT_TOMOGRAPHY_TOOLKIT/ITERATION_UPDATE/X04_SRC_DIRECTION_LBFGS/compute_direction_lbfgs.f90
Log:
change the initial approximate Hessian

Modified: seismo/3D/ADJOINT_TOMO/ADJOINT_TOMOGRAPHY_TOOLKIT/ITERATION_UPDATE/X04_SRC_DIRECTION_LBFGS/compute_direction_lbfgs.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/ADJOINT_TOMOGRAPHY_TOOLKIT/ITERATION_UPDATE/X04_SRC_DIRECTION_LBFGS/compute_direction_lbfgs.f90	2013-07-30 19:19:21 UTC (rev 22687)
+++ seismo/3D/ADJOINT_TOMO/ADJOINT_TOMOGRAPHY_TOOLKIT/ITERATION_UPDATE/X04_SRC_DIRECTION_LBFGS/compute_direction_lbfgs.f90	2013-07-30 20:07:39 UTC (rev 22688)
@@ -105,6 +105,16 @@
      q_vector=q_vector-a(istore)*gradient_diff 
   end do 
 
+  !
+  istore=iter_current-1
+  call get_gradient(istore+1,gradient1)
+  call get_gradient(istore,gradient0)
+  call get_model(istore+1,model1)
+  call get_model(istore,model0)
+  gradient_diff=gradient1-gradient0 
+  model_diff=model1-model0
+  !
+
   p_k_up=sum(gradient_diff*model_diff)
   p_k_down=sum(gradient_diff*gradient_diff)
   call mpi_barrier(MPI_COMM_WORLD,ier)