INNER CODE UNIT · Python
delta
XMU-Kuangnan-Fang-Team/GENetLib · GENetLib/BsplineFunc.py:123
delta = breaks[k + 1] - breaks[k]
power = delta
prodmati = np.zeros((norder, norder))
for i in range(1, prodorder + 1):
prodmati += power * convmat[:, :, prodorder - i] / i
power *= delta
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]: