woggioni/jwo
1
0
Fork 0
Code Actions Packages Activity

added unit test

Browse Source
This commit is contained in:
Walter Oggioni
2023-05-28 19:59:42 +08:00
parent 0663ca200a
commit 5b3e279c20
1 changed files with 20 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
jmath/src/test/java/net/woggioni/jmath/MatrixTest.java
View File

@@ -4,6 +4,9 @@ import org.junit.jupiter.api.Assertions;
import org.junit.jupiter.api.BeforeAll; import org.junit.jupiter.api.BeforeAll;
import org.junit.jupiter.api.Test; import org.junit.jupiter.api.Test;
import java.util.Optional;
import java.util.Random;
import java.util.function.IntFunction;
import java.util.stream.IntStream; import java.util.stream.IntStream;
public class MatrixTest { public class MatrixTest {
@@ -203,4 +206,21 @@ public class MatrixTest {
Assertions.assertEquals(expectedSolution, inverse.mmul(b)); Assertions.assertEquals(expectedSolution, inverse.mmul(b));
} }
@Test
public void randomLinearSystemTest() {
int size = 10;
Random rnd = new Random(101325);
NumericTypeFactory<Rational> numericTypeFactory = RationalFactory.getInstance();
Matrix.ValueGenerator<Rational> init = (i, j) ->
Rational.of(rnd.nextInt(-1000 * size, 1000 * size), (long) 1000 * size);
Matrix<Rational> mtx = Matrix.of(numericTypeFactory, size, size, init);
Matrix<Rational> lu = mtx.clone();
Matrix.Pivot pivot = lu.lup();
IntFunction<Rational> initVector = (i) -> Rational.of(rnd.nextInt(0, size), size);
Vector<Rational> b = Vector.of(numericTypeFactory, size, initVector);
Vector<Rational> x = lu.luSolve(b, pivot);
Vector<Rational> error = mtx.mmul(x).sub(b);
Assertions.assertEquals(Rational.ZERO, error.norm());
}
} }