From 288605bf7151fda07125fd0bc0e3571ed8bb5516 Mon Sep 17 00:00:00 2001 From: jthorn Date: Wed, 28 Aug 2002 11:31:09 +0000 Subject: In Hermite_polynomial_interpolant(), add explicit check for the case when the given conditions don't uniquely determine the interpolant git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@100 df1f8a13-aa1d-4dd4-9681-27ded5b42416 --- src/GeneralizedPolynomial-Uniform/interpolate.maple | 5 +++++ 1 file changed, 5 insertions(+) (limited to 'src/GeneralizedPolynomial-Uniform') diff --git a/src/GeneralizedPolynomial-Uniform/interpolate.maple b/src/GeneralizedPolynomial-Uniform/interpolate.maple index 58a2c97..71280fe 100644 --- a/src/GeneralizedPolynomial-Uniform/interpolate.maple +++ b/src/GeneralizedPolynomial-Uniform/interpolate.maple @@ -208,6 +208,11 @@ deriv_eqnset := `union`(op(%)); # in terms of DATA() and DERIV() values # coeff_eqns := solve[linear](fn_eqnset union deriv_eqnset, coeff_set); +if (indets(map(rhs,%)) <> {}) + then error "no unique solution for coefficients -- %1 eqns for %2 coeffs", + nops(fn_eqnset union deriv_eqnset), + nops(coeff_set); +fi; # -- cgit v1.2.3