def trace(self): """ Trace of the smoothing matrix S(pen) """ if self.pen > 0: _invband = _bspline.invband(self.chol.copy()) tr = _trace_symbanded(_invband, self.btb, lower=1) return tr else: return self.rank
def trace(self): """ Trace of the smoothing matrix S(pen) TODO: addin a reference to Wahba, and whoever else I used. """ if self.pen > 0: _invband = _bspline.invband(self.chol.copy()) tr = _trace_symbanded(_invband, self.btb, lower=1) return tr else: return self.rank