Skip to content

Commit f1d70ca

Browse files
ktkonradktkonrad
committed
WIP on (no branch): 3c8043e rewrote algorithm to interpolate along all nodal domain boundaries moved some util functions around simplified masking to just use counted matrix implemented a bit array to store signs of upsampled grid create interpolation workspace struct to avoid repetitive mallocs
2 parents 3c8043e + 167efc8 commit f1d70ca

20 files changed

+97945
-96511
lines changed

‎c/nodal_domain_driver.c

+1-1
Original file line numberDiff line numberDiff line change
@@ -277,7 +277,7 @@ int main(int argc, char **argv) {
277277
exit(IO_ERR);
278278
}
279279

280-
counted = createScaledMaskFromBilliard(bil, dx/upsample_ratio, k/k_0, &counted_y, &counted_x);
280+
counted = createScaledMaskFromBilliard(bil, 1.0/(((ny-1)*upsample_ratio)+0), k/k_0, &counted_y, &counted_x);
281281

282282
if (counted_x != ((nx-1)*upsample_ratio)+1 || counted_y != ((ny-1)*upsample_ratio)+1) {
283283
ERROR("mask dimensions do not match grid dimensions\n \

‎c/nodal_domain_driver_no_main.c

+1-1
Original file line numberDiff line numberDiff line change
@@ -278,7 +278,7 @@ int count_main(int argc, char **argv) {
278278
exit(IO_ERR);
279279
}
280280

281-
counted = createScaledMaskFromBilliard(bil, dx/upsample_ratio, k/k_0, &counted_y, &counted_x);
281+
counted = createScaledMaskFromBilliard(bil, 1.0/(((ny-1)*upsample_ratio)+0), k/k_0, &counted_y, &counted_x);
282282

283283
if (counted_x != ((nx-1)*upsample_ratio)+1 || counted_y != ((ny-1)*upsample_ratio)+1) {
284284
ERROR("mask dimensions do not match grid dimensions\n \

‎documents/thesis/figs/interpolation/error_norms_2.eps

+43-43
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.

‎documents/thesis/figs/interpolation/error_norms_3.eps

+43-43
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.

0 commit comments

Comments
 (0)