From 6978a1dcbf1972ad7713f3511dc2325c1ca111d0 Mon Sep 17 00:00:00 2001 From: jjones Date: Fri, 28 Mar 2003 18:56:11 +0000 Subject: [PATCH] Bug fix for red-black relaxation --- sstruct_ls/node_relax.c | 4 ++++ struct_ls/point_relax.c | 4 ++++ 2 files changed, 8 insertions(+) diff --git a/sstruct_ls/node_relax.c b/sstruct_ls/node_relax.c index 7a1fb7c14..2df6eb32e 100644 --- a/sstruct_ls/node_relax.c +++ b/sstruct_ls/node_relax.c @@ -659,6 +659,10 @@ hypre_NodeRelax( void *relax_vdata, if (zero_guess) { + if (num_nodesets > 1) + { + hypre_SStructPVectorSetConstantValues(x, 0.0); + } nodeset = nodeset_ranks[p]; compute_pkg = compute_pkgs[nodeset]; stride = nodeset_strides[nodeset]; diff --git a/struct_ls/point_relax.c b/struct_ls/point_relax.c index 63ffcc924..94fdc7d8b 100644 --- a/struct_ls/point_relax.c +++ b/struct_ls/point_relax.c @@ -408,6 +408,10 @@ hypre_PointRelax( void *relax_vdata, if (zero_guess) { + if (num_pointsets > 1) + { + hypre_StructVectorSetConstantValues(x, 0.0); + } pointset = pointset_ranks[p]; compute_pkg = compute_pkgs[pointset]; stride = pointset_strides[pointset];