Method deform

pythran/fem/FE_model.py:199–227

dV is a linear list or array of deformations to apply to the nodal coordinates. Even indices are the x component, odd are y. scale = None means autoscale so that deformations are 30% of the model's max dimension.

(self, dV,                                    # {{{
               scale=None)

            print(S)
    # }}}
    def deform(self, dV,                                    # {{{
               scale=None):
        """
        dV is a linear list or array of deformations to
        apply to the nodal coordinates.
        Even indices are the x component, odd are y.

        scale = None means autoscale so that deformations
                are 30% of the model's max dimension.
        """
        if scale is None:
            model_len = np.sqrt((self.xmax - self.xmin)**2 +
                                (self.ymax - self.ymin)**2)
            S = 0.1 * model_len / np.max(np.abs(dV))
        else:
            S = scale

       #print(f'max dV {np.max(np.abs(dV))}, model_len {model_len}  S {S}')
        for dX, dY, N in zip(dV[0::2],
                             dV[1::2],
                             self.node):
            N.x += S*dX
            N.y += S*dY

        # update the bounding box to the deformed shape
        self.xmin = min([_.x for _ in self.node])
        self.xmax = max([_.x for _ in self.node])
        self.ymin = min([_.y for _ in self.node])
        self.ymax = max([_.y for _ in self.node])
    # }}}