X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M-SAXS-homology%2Fminimize_p.F;h=f9faf7c2bf08e37cef6b326257f18868e6801ba0;hb=34d3ad3987785642be58fb2f26557d3314215577;hp=3de233fc11fb583aa35f24ee0f102044d6492ae9;hpb=f690e8b70bab14132839afebf080d4a28363b226;p=unres.git diff --git a/source/unres/src_MD-M-SAXS-homology/minimize_p.F b/source/unres/src_MD-M-SAXS-homology/minimize_p.F index 3de233f..f9faf7c 100644 --- a/source/unres/src_MD-M-SAXS-homology/minimize_p.F +++ b/source/unres/src_MD-M-SAXS-homology/minimize_p.F @@ -24,7 +24,7 @@ external func,gradient,fdum external func_restr,grad_restr logical not_done,change,reduce - common /przechowalnia/ v +c common /przechowalnia/ v icall = 1 @@ -70,7 +70,7 @@ c v(25)=4.0D0 if (rtolf.eq.0.0D0) rtolf=1.0D-4 v(32)=rtolf * controls initial step size - v(35)=1.0D-1 + v(35)=1.0D-1 * large vals of d correspond to small components of step do i=1,nphi d(i)=1.0D-1 @@ -256,7 +256,8 @@ c call MATOUT2(my_ng_count,dimen3,maxres2,maxers2,ginv_block) real*8 urparm(1) dimension x(maxvar) c if (jjj.gt.0) then -c write (iout,'(10f8.3)') (rad2deg*x(i),i=1,n) +c write (iout,*) "in func x" +c write (iout,'(10f8.3)') (rad2deg*x(i),i=1,n) c endif nfl=nf icg=mod(nf,2)+1 @@ -269,8 +270,9 @@ cd write (iout,*) 'ETOTAL called from FUNC' call sum_gradient f=energia(0) c if (jjj.gt.0) then -c write (iout,'(10f8.3)') (rad2deg*x(i),i=1,n) -c write (iout,*) 'f=',etot +c write (iout,*) "upon exit from func" +c write (iout,'(10f8.3)') (rad2deg*x(i),i=1,n) +c write (iout,*) 'f=',f c jjj=0 c endif return