FEAP User Forum

FEAP => Programming => Topic started by: weiyongtao on February 28, 2021, 05:17:05 PM

Title: is there a bug in the subroutine settfb in the file of framf3d.f
Post by: weiyongtao on February 28, 2021, 05:17:05 PM

In the subroutine of settfb, for the array tmp(3,3), it is orginal set as following:

      tmp(1,1) =  0.0d0
      tmp(2,2) =  0.0d0
      tmp(3,3) =  0.0d0

      tmp(1,2) = -x(3)*c4
      tmp(2,1) =  tmp(1,2)

      tmp(1,3) =  x(2)*c4
      tmp(3,1) =  tmp(1,3)

      tmp(2,3) = -x(1)*c4
      tmp(3,2) =  tmp(2,3)
I guess those statements above set tmp as a skew-symmetric array, if it is true, then statements setting tmp(1,3) and tmp(3,1) may be a bug.

Title: Re: is there a bug in the subroutine settfb in the file of framf3d.f
Post by: Prof. R.L. Taylor on February 28, 2021, 09:41:15 PM

It is a bug