FEAP User Forum

FEAP => Programming => Topic started by: Ziyu on May 03, 2013, 11:29:02 AM

Title: Bugs in pmacr2.f?
Post by: Ziyu on May 03, 2013, 11:29:02 AM

Hi,

If you look at the line 871-872 under commands [newf] and [newf, zero]:

            do n = i,nneq-1,ndf
              hr(np(28)+i+2*nneq) = 0.0d0  ! f0(n,1) = fixed forces
              hr(np(28)+i+3*nneq) = 0.0d0  ! f0(n,2) = fixed displacements
            end do ! n

It seems like the loop index i  in hr(np(28)+i+2*nneq) and hr(np(28)+i+3*nneq) = 0.0d0 should be replaced with n.

Best,

Ziyu

Title: Re: Bugs in pmacr2.f?
Post by: Prof. R.L. Taylor on May 04, 2013, 05:21:55 PM

You are correct, thank you