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.doublematrix.calculation.entrywise.creators;
025    
026    import java.util.Arrays;
027    
028    import org.ujmp.core.Coordinates;
029    import org.ujmp.core.Matrix;
030    import org.ujmp.core.MatrixFactory;
031    import org.ujmp.core.doublematrix.calculation.AbstractDoubleCalculation;
032    import org.ujmp.core.enums.ValueType;
033    import org.ujmp.core.exceptions.MatrixException;
034    
035    public class Eye extends AbstractDoubleCalculation {
036            private static final long serialVersionUID = 2547827499345834225L;
037    
038            public Eye(Matrix matrix) {
039                    super(matrix);
040            }
041    
042            public double getDouble(long... coordinates) {
043                    return coordinates[ROW] == coordinates[COLUMN] ? 1.0 : 0.0;
044            }
045    
046            public static Matrix calc(Matrix source) throws MatrixException {
047                    Matrix ret = Matrix.factory.zeros(source.getSize());
048                    long[] c = Coordinates.copyOf(source.getSize());
049                    for (int i = 0; Coordinates.isSmallerThan(c, source.getSize()); i++) {
050                            Arrays.fill(c, i);
051                            ret.setAsDouble(1.0, c);
052                    }
053                    return ret;
054            }
055    
056            public static Matrix calcOrig(Matrix source) throws MatrixException {
057                    long[] c = new long[source.getSize().length];
058                    for (int i = 0; Coordinates.isSmallerThan(c, source.getSize()); i++) {
059                            source.setAsDouble(1.0, c);
060                            Arrays.fill(c, i + 1);
061                    }
062                    return source;
063            }
064    
065            public static Matrix calc(long... size) throws MatrixException {
066                    return calc(ValueType.DOUBLE, size);
067            }
068    
069            public static Matrix calc(ValueType valueType, long... size) throws MatrixException {
070                    Matrix ret = MatrixFactory.zeros(valueType, size);
071                    long[] c = new long[size.length];
072                    for (int i = 0; Coordinates.isSmallerThan(c, size); i++) {
073                            ret.setAsDouble(1.0, c);
074                            Arrays.fill(c, i + 1);
075                    }
076                    return ret;
077            }
078    }