001    /*
002     * Copyright (C) 2008-2010 by Holger Arndt
003     * 
004     * This file is part of the Universal Java Matrix Package (UJMP).
005     * See the NOTICE file distributed with this work for additional
006     * information regarding copyright ownership and licensing.
007     *
008     * UJMP is free software; you can redistribute it and/or modify
009     * it under the terms of the GNU Lesser General Public License as
010     * published by the Free Software Foundation; either version 2
011     * of the License, or (at your option) any later version.
012     *
013     * UJMP is distributed in the hope that it will be useful,
014     * but WITHOUT ANY WARRANTY; without even the implied warranty of
015     * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
016     * GNU Lesser General Public License for more details.
017     *
018     * You should have received a copy of the GNU Lesser General Public
019     * License along with UJMP; if not, write to the
020     * Free Software Foundation, Inc., 51 Franklin St, Fifth Floor,
021     * Boston, MA  02110-1301  USA
022     */
023    
024    package org.ujmp.core.benchmark;
025    
026    import org.ujmp.core.Coordinates;
027    import org.ujmp.core.Matrix;
028    import org.ujmp.core.doublematrix.DoubleMatrix2D;
029    import org.ujmp.core.matrix.DenseMatrix;
030    import org.ujmp.core.util.SerializationUtil;
031    
032    public class InvBenchmarkTask extends AbstractBenchmarkTask {
033    
034            public InvBenchmarkTask(long benchmarkSeed, Class<? extends DoubleMatrix2D> matrixClass,
035                            BenchmarkConfig config) {
036                    super(benchmarkSeed, matrixClass, config.getInvSizes(), config);
037            }
038    
039            @Override
040            public BenchmarkResult task(Class<? extends Matrix> matrixClass, long benchmarkSeed, int run,
041                            long[] size) {
042                    final long t0, t1, m0, m1;
043                    final DoubleMatrix2D m;
044                    Matrix r;
045                    try {
046                            m = BenchmarkUtil.createMatrix(matrixClass, size);
047                            if (!m.getClass().getName().startsWith("org.ujmp.core")
048                                            && m.getClass().getDeclaredMethod("inv") == null) {
049                                    System.out.print("-");
050                                    System.out.flush();
051                                    return BenchmarkResult.NOTAVAILABLE;
052                            }
053                            BenchmarkUtil.rand(benchmarkSeed, run, 0, m);
054                            BenchmarkUtil.purgeMemory(getConfig());
055                            m0 = Runtime.getRuntime().totalMemory() - Runtime.getRuntime().freeMemory();
056                            t0 = System.nanoTime();
057                            r = m.inv();
058                            t1 = System.nanoTime();
059                            m1 = Runtime.getRuntime().totalMemory() - Runtime.getRuntime().freeMemory();
060                            if (r == null) {
061                                    System.out.print("e");
062                                    System.out.flush();
063                                    return BenchmarkResult.ERROR;
064                            }
065                            Matrix result = m.mtimes(r);
066                            double diff = BenchmarkUtil.difference(result, DenseMatrix.factory.eye(m.getSize()));
067                            result = null;
068                            long mem = m1 - m0 - SerializationUtil.sizeOf(r);
069                            mem = mem > 0 ? mem : 0;
070                            return new BenchmarkResult((t1 - t0) / 1000000.0, diff, mem);
071                    } catch (Throwable e) {
072                            System.out.print("e");
073                            System.out.flush();
074                            return BenchmarkResult.ERROR;
075                    }
076            }
077    
078            @Override
079            public String getTaskName() {
080                    return "inv";
081            }
082    
083    }