{
let n = 40;
let x = new Float64Array(n);
sparseLU.Solve(b, x);
return test("Solution is correct", () => {
for (let i = 0; i < n; i++) {
if (assertAlmostEqual(x[i] - xExpected[i], 1e-10)) {
throw `invalid result at position ${i}:\n actual = ${x[i]}\n expected = ${xExpected[i]}\n error = ${x[i] - xExpected[i]}`;
}
}
});
}