SeqAggTest.java

package org.openlca.eigen.solvers;

import org.junit.Assert;
import org.junit.Test;
import org.openlca.eigen.HashMatrix;

public class SeqAggTest {

	@Test
	public void testSimpleSystem() {
		HashMatrix a = new HashMatrix(new double[][] {
				{ 2., -5., 0., 0., 0. }, { 0., 1., -2., 0., 0. },
				{ 0., 0., 1., 0., -6. }, { 0., 0., -1., 4., 0. },
				{ 0., 0., 0., 0., 1. }, });
		HashMatrix b = new HashMatrix(new double[][] { { 3, 2, 1, 1, 3 },
				{ 1, 2, 2, 1, 1 } });
		SeqAgg seqAgg = new SeqAgg(a, b, 4, 1);
		double[] g = seqAgg.solve();
		Assert.assertArrayEquals(new double[] { 124.5, 68.5 }, g, 1e-16);
	}

	@Test
	public void testLoopSystem() {
		HashMatrix a = new HashMatrix(
				new double[][] { { 4., -2., 0., 0. }, { 0., 1., -0.1, -2. },
						{ 0., -1., 2., 0. }, { 0., 0., 0., 1. } });
		HashMatrix b = new HashMatrix(new double[][] { { 2., 3., 2., 1. },
				{ 3., 1., 2., 1., } });
		SeqAgg seqAgg = new SeqAgg(a, b, 3, 1);
		double[] g = seqAgg.solve();
		Assert.assertArrayEquals(new double[] { 11.526316, 8.3684211 }, g, 1e-6);
	}

}