INNER CODE UNIT · Python
index
XMU-Kuangnan-Fang-Team/GENetLib · GENetLib/BsplineFunc.py:129
index = np.arange(k, k + norder)
prodmat[index[:, None], index] += prodmati
penaltymat = prodmat
else:
if uniquebreaks:
prodmat = inprod(
self.basisobj, self.basisobj, Lfdobj, Lfdobj, rng
)
else:
rngvec = [rng[0]]
for i in range(1, nbreaks):
if breaks[i] == breaks[i - 1]:
rngvec.append(breaks[i])
rngvec = list(set(rngvec))
nrng = len(rngvec)
if rngvec[-1] < rng[1]:
rngvec.append(rng[1])
nrng += 1