org.la4j.matrix.dense.AbstractBasicMatrix Maven / Gradle / Ivy

Go to download

Show more of this group Show more artifacts with this name
Show all versions of la4j Show documentation

A 100% Java sparse and dense matrix library.

There is a newer version: 0.6.0

/*
 * Copyright 2011-2013, by Vladimir Kostyukov and Contributors.
 * 
 * This file is part of la4j project (http://la4j.org)
 * 
 * Licensed under the Apache License, Version 2.0 (the "License");
 * You may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 * 
 *      http://www.apache.org/licenses/LICENSE-2.0
 * 
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 * 
 * Contributor(s): -
 * 
 */

package org.la4j.matrix.dense;

import org.la4j.factory.Factory;
import org.la4j.matrix.AbstractMatrix;
import org.la4j.matrix.Matrix;

public abstract class AbstractBasicMatrix extends AbstractMatrix 
    implements DenseMatrix {

    private static final int BLOCKSIZE = 64;

    public AbstractBasicMatrix(Factory factory, int rows, int columns) {
        super(factory, rows, columns);
    }

    @Override
    public Matrix multiply(Matrix matrix, Factory factory) {

        if (matrix instanceof DenseMatrix && (rows % BLOCKSIZE == 0) 
                && (columns % BLOCKSIZE == 0) 
                && (matrix.columns() % BLOCKSIZE == 0)) {

            return multiplyBlockedWith64(matrix, factory);
        }

        return super.multiply(matrix, factory);
    }

    private Matrix multiplyBlockedWith64(Matrix matrix, Factory factory) {

        Matrix result = factory.createMatrix(rows, matrix.columns());

        for (int i = 0; i < rows; i += BLOCKSIZE) {
            for (int k = 0; k < columns; k += BLOCKSIZE) {
                for (int j = 0; j < matrix.columns(); j += BLOCKSIZE) {
                    for(int u = 0; u < BLOCKSIZE; u++) {
                        for(int w = 0; w < BLOCKSIZE; w++) {
                            for(int v = 0; v < BLOCKSIZE; v++) {
                                result.set(i + u, j + v, 
                                   result.get(i + u, j + v) 
                                   + (get(i + u, k + w)) 
                                   * matrix.get(k + w, j + v));
                            }
                        }
                    }
                }
            }
        }

        return result;
    }
}