Skip to content
Snippets Groups Projects
Commit 53324c34 authored by ronesy's avatar ronesy
Browse files

bug fix in calculating cost for lognormal in m1qn3_interface

parent 98d1bd5c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
source/m1qn3_interface.f90
+ 2
− 2
View file @ 53324c34
...@@ -190,11 +190,11 @@ subroutine m1qn3_interface(iter, grad, files, config, fluxes, obs, states, covar ...@@ -190,11 +190,11 @@ subroutine m1qn3_interface(iter, grad, files, config, fluxes, obs, states, covar
if ( config%lognormal ) then if ( config%lognormal ) then
! lognormal ! lognormal
! 2Jp = 2*sum(z) + x^Tx (where x = state vector in chi space and z is state vector in physical space) ! 2Jp = 2*sum(z) + x^Tx (where x = state vector in chi space and z is state vector in physical space)
cost_p = sum(states%px) + dot_product(x,x) cost_p = 2.*sum(states%px) + real(dot_product(x,x),kind=4)
else else
! normal ! normal
! 2Jp = x^Tx (where x = state vector in chi space) ! 2Jp = x^Tx (where x = state vector in chi space)
cost_p = dot_product(x,x) cost_p = real(dot_product(x,x),kind=4)
endif endif
! J = 0.5*(Jp + Jo) ! J = 0.5*(Jp + Jo)
cost = 0.5*(cost_p + cost_o) cost = 0.5*(cost_p + cost_o)
......
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