write (iout,*) i,(duscdiffx(j,i),j=1,3)
enddo
#endif
-c if(usampl.and.totT.gt.eq_time) then
- if(usampl.and.totT.gt.eq_time .or. constr_homology.gt.0) then
+ if(usampl.and.totT.gt.eq_time) then
+c if(usampl.and.totT.gt.eq_time .or. constr_homology.gt.0) then
c
c Setting suited bounds for HM restrs
c
do i=1,nct
do j=1,3
- gradc(j,i,icg)=gradc(j,i,icg)+dudconst(j,i)+duscdiff(j,i)
- gradx(j,i,icg)=gradx(j,i,icg)+dudxconst(j,i)+duscdiffx(j,i)
+ gradc(j,i,icg)=gradc(j,i,icg)+dudconst(j,i)
+ gradx(j,i,icg)=gradx(j,i,icg)+dudxconst(j,i)
enddo
#ifdef DEBUG
write (iout,*) "i",i," gradc",(gradc(j,i,icg),j=1,3)