diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-28 11:31:09 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-28 11:31:09 +0000 |
commit | 288605bf7151fda07125fd0bc0e3571ed8bb5516 (patch) | |
tree | 3e7838ec869569e7005cca47a4fd52e2b9a547c6 | |
parent | afcc4e7c56586ac436b81763ab145d1ef3feea65 (diff) |
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
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/interpolate.maple | 5 |
1 files changed, 5 insertions, 0 deletions
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; # |