call fivediaginv_mult(dimen,fric_work, d_af_work)
c write (iout,*) "stochastic acceleratios"
call fivediaginv_mult(dimen,stochforcvec, d_as_work)
-c write (iout,*) "Leaving sddir_precalc"
#else
call ginv_mult(fric_work, d_af_work)
call ginv_mult(stochforcvec, d_as_work)
write (iout,'(3f10.5)') (d_af_work(i),i=1,dimen3)
write (iout,*) "d_as_work"
write (iout,'(3f10.5)') (d_as_work(i),i=1,dimen3)
+ write (iout,*) "Leaving sddir_precalc"
#endif
return
end
double precision xv,sigv,lowb,highb,vec_afm(3),Ek1,Ek2,Ek3,aux
integer i,ii,j,k,l,ind
double precision anorm_distr
- logical lprn /.true./
+ logical lprn /.false./
#ifdef FIVEDIAG
integer ichain,n,innt,inct,ibeg,ierr
double precision work(8*maxres6)
dc(j,0)=dc_work(j)
d_t(j,0)=d_t_work(j)
enddo
- ind=3
- do i=nnt,nct-1
+ ind=3
+ do i=nnt,nct-1
do j=1,3
dc(j,i)=dc_work(ind+j)
d_t(j,i)=d_t_work(ind+j)