INNER CODE UNIT · Python

penalty_matrix

XMU-Kuangnan-Fang-Team/GENetLib · GENetLib/BsplineFunc.py:16

    def penalty_matrix(self, btype="spline"):
        if btype == "spline":

            # Initialize parameters
            if type(self.Lfdobj) == int:
                Lfdobj = int2lfd(self.Lfdobj)
            else:
                Lfdobj = self.Lfdobj
            nbasis = self.basisobj["nbasis"]
            params = self.basisobj["params"]
            Rangeval = self.basisobj["rangeval"]
            Rng = self.rng if self.rng is not None else Rangeval
            rng = np.array(Rng, dtype=float)
            rangeval = np.array(Rangeval, dtype=float)
            params = np.array(params, dtype=float)
            breaks = np.concatenate(([rangeval[0]], params, [rangeval[1]]))
            nbreaks = len(breaks)
            ninterval = nbreaks - 1

View source record →