Skip to content

use values from FEniCS rather than half of values from Bempp-cl #27

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
mscroggs merged 2 commits into from
Jan 23, 2026
Merged

use values from FEniCS rather than half of values from Bempp-cl #27

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
323ac2f
use values from FEniCS rather than half of values from Bempp-cl
mscroggs Jan 23, 2026
a298fca
fmt
mscroggs Jan 23, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 5 additions & 8 deletions ndfunctionspace/examples/test_mass_matrix.rs
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
//! Test values in mass matrices against values computing using Bempp
//! Test values in mass matrices against values computing using Bempp--cl and FEniCS

use approx::assert_relative_eq;
use ndelement::{
Expand Down Expand Up @@ -69,6 +69,7 @@ fn test_lagrange_mass_matrix() {
}
}

// Compare matrix entries to values from Bempp
for i in 0..6 {
assert_relative_eq!(mass_matrix[[i, i]], 0.5773502691896255, epsilon = 1e-10);
}
Expand Down Expand Up @@ -146,20 +147,16 @@ fn test_rt_mass_matrix() {
}
}

// The values from Bempp-cl here are divided by 2 as the basis functions in -cl are not scaled by edge length
// Compare matrix entries to FEniCS
for i in 0..12 {
assert_relative_eq!(
mass_matrix[[i, i]],
0.9622504486493761 / 2.0,
epsilon = 1e-10
);
assert_relative_eq!(mass_matrix[[i, i]], 0.4811252243246884, epsilon = 1e-10);
}
for i in 0..12 {
for j in 0..12 {
if i != j && mass_matrix[[i, j]].abs() > 0.001 {
assert_relative_eq!(
mass_matrix[[i, j]].abs(),
0.09622504486493733 / 2.0,
0.0481125224324689,
epsilon = 1e-10
);
}
Expand Down