Contents
functionpetsc-3.4.2 2013-07-02[xl,xu] = MSA_Plate(xl,xu,user)
This file only implements a rectangular plate.
xl.Set(-10000000); xu.Set(100000000); xl_2d = user.dm.VecGetArray(xl); bmx = user.bmx; bmy = user.bmy; mx = user.mx; my = user.my; bheight = user.bheight; for (j = 1:my) for (i = 1:mx) if ((i >(mx-bmx)/2) && (i <= mx-(mx-bmx)/2) && (j > (my-bmy)/2) && (j <= my-(my-bmy)/2)) xl_2d(i,j) = bheight; end end end