Merge pull request #501 from pdziekan/init_laplacian_edges_fix

pdziekan · web-flow · commit 1c701f2306fa · 2026-03-03T09:40:51.000+01:00
prs solvers laplacian calc (with err_init): move set_edges before G m…
diff --git a/libmpdata++/solvers/detail/mpdata_rhs_vip_prs_common.hpp b/libmpdata++/solvers/detail/mpdata_rhs_vip_prs_common.hpp
@@ -63,6 +63,7 @@ namespace libmpdataxx
               lap_tmp[d](this->ijk) -= tmp_uvw[d](this->ijk);
             }
           }
+          this->set_edges(lap_tmp, this->ijk, err_init ? -1 : 0);
           if (this->mem->G)
           {
             for (int d = 0; d < parent_t::n_dims; ++d)
@@ -71,7 +72,6 @@ namespace libmpdataxx
             }
           }
           if (!simple) this->normalize_vip(lap_tmp);
-          this->set_edges(lap_tmp, this->ijk, err_init ? -1 : 0);
           for (int d = 0; d < parent_t::n_dims; ++d)
           {
             this->xchng_pres(lap_tmp[d], ijk);

Original file line number	Diff line number	Diff line change
`@@ -63,6 +63,7 @@ namespace libmpdataxx`
`63`	`63`	`lap_tmp[d](this->ijk) -= tmp_uvw[d](this->ijk);`
`64`	`64`	`}`
`65`	`65`	`}`
	`66`	`+ this->set_edges(lap_tmp, this->ijk, err_init ? -1 : 0);`
`66`	`67`	`if (this->mem->G)`
`67`	`68`	`{`
`68`	`69`	`for (int d = 0; d < parent_t::n_dims; ++d)`
`@@ -71,7 +72,6 @@ namespace libmpdataxx`
`71`	`72`	`}`
`72`	`73`	`}`
`73`	`74`	`if (!simple) this->normalize_vip(lap_tmp);`
`74`		`- this->set_edges(lap_tmp, this->ijk, err_init ? -1 : 0);`
`75`	`75`	`for (int d = 0; d < parent_t::n_dims; ++d)`
`76`	`76`	`{`
`77`	`77`	`this->xchng_pres(lap_tmp[d], ijk);`