From 66ae32b089700585a43d47dae56a0bff10ee1f0d Mon Sep 17 00:00:00 2001 From: Ralf Hartmann Date: Wed, 10 Jan 2001 17:04:43 +0000 Subject: [PATCH] Fix small bug in the Horner scheme that produced wrong results only very scarce. git-svn-id: https://svn.dealii.org/trunk@3634 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/source/polynomial.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/deal.II/base/source/polynomial.cc b/deal.II/base/source/polynomial.cc index 110b448605..419ebe72be 100644 --- a/deal.II/base/source/polynomial.cc +++ b/deal.II/base/source/polynomial.cc @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 2000 by the deal.II authors +// Copyright (C) 2000, 2001 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -72,7 +72,7 @@ void Polynomial::value (const double x, unsigned int j_faculty=1; for (unsigned int j=0; j=static_cast(j); --k) + for (int k=m-2; k>=static_cast(j); --k) a[k]+=x*a[k+1]; values[j]=j_faculty*a[j]; -- 2.39.5