INNER CODE UNIT · Python
m
XMU-Kuangnan-Fang-Team/GENetLib · GENetLib/BsplineFunc.py:94
m = norder - 1
for j in range(nrowcoef):
polycoef[i - k + j, :, m - j] = CoeffD[j, :]
prodmat = np.zeros((nbasis, nbasis))
convmat = np.zeros((norder, norder, prodorder))
for k in range(ninterval):
Coeff = polycoef[k, :, :]
for i in range(ndegree):
ind = np.arange(i + 1)
if len(ind) == 1:
convmat[:, :, i] = np.outer(
Coeff[ind, :], Coeff[i - ind, :]
)
convmat[:, :, prodorder - i - 1] = np.outer(
Coeff[ndegree - ind, :],
Coeff[ndegree - i + ind, :],
)
else: