From e048daa436c38079c64be376e8c2bc37490ffc99 Mon Sep 17 00:00:00 2001 From: tomsmeding Date: Sat, 20 Apr 2019 22:19:13 +0200 Subject: It's Jacobi-style, not Lagrange-style --- aberth/aberth_kernel.fut | 2 +- aberth/host_aberth.cpp | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/aberth/aberth_kernel.fut b/aberth/aberth_kernel.fut index 79d7cfa..fd4d943 100644 --- a/aberth/aberth_kernel.fut +++ b/aberth/aberth_kernel.fut @@ -64,7 +64,7 @@ module aberth = { let radius = max_root_norm p in {p, deriv, bound, radius} - -- Lagrange-style step where the new elements are computed in parallel from + -- Jacobi-style step where the new elements are computed in parallel from -- the previous values let step (ctx: context) (approx: *approx): *(bool, approx) = let pvals = map (eval_c ctx.p) approx diff --git a/aberth/host_aberth.cpp b/aberth/host_aberth.cpp index c567f30..fcffebf 100644 --- a/aberth/host_aberth.cpp +++ b/aberth/host_aberth.cpp @@ -42,7 +42,7 @@ void AberthState::regenerate() { } } -// Lagrange-style step where the new elements are computed in parallel from the previous values +// Jacobi-style step where the new elements are computed in parallel from the previous values bool AberthState::step() { array pairs; for (int i = 0; i < N - 1; i++) { -- cgit v1.2.3