@@ -0,0 +1,2512 @@
+/*%0A * Zorbage: an algebraic data hierarchy for use in numeric processing.%0A *%0A * Copyright (C) 2016-2018 Barry DeZonia%0A * %0A * Redistribution and use in source and binary forms, with or without%0A * modification, are permitted provided that the following conditions are met:%0A * %0A * 1. Redistributions of source code must retain the above copyright notice,%0A *    this list of conditions and the following disclaimer.%0A * 2. Redistributions in binary form must reproduce the above copyright notice,%0A *    this list of conditions and the following disclaimer in the documentation%0A *    and/or other materials provided with the distribution.%0A * %0A * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS %22AS IS%22%0A * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE%0A * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE%0A * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDERS OR CONTRIBUTORS BE%0A * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR%0A * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF%0A * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS%0A * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN%0A * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)%0A * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE%0A * POSSIBILITY OF SUCH DAMAGE.%0A */%0Apackage nom.bdezonia.zorbage.algorithm;%0A%0Aimport nom.bdezonia.zorbage.type.algebra.Group;%0Aimport nom.bdezonia.zorbage.type.algebra.MatrixMember;%0Aimport nom.bdezonia.zorbage.type.algebra.Multiplication;%0A%0A/**%0A * %0A * @author Barry DeZonia%0A *%0A */%0Apublic class MatrixDirectProduct %7B%0A%0A%09private MatrixDirectProduct() %7B %7D%0A%09%0A%09public static %3CT extends Group%3CT,U%3E & Multiplication%3CU%3E,U,W extends MatrixMember%3CU%3E%3E%0A%09%09void compute(T group, W in1, W in2, W out)%0A%09%7B%0A%09%09U tmp1 = group.construct();%0A%09%09U tmp2 = group.construct();%0A%09%09U tmp3 = group.construct();%0A%09%09if (out == in1 %7C%7C out == in2)%0A%09%09%09throw new IllegalArgumentException(%22output matrix must not be one of input matrices%22);%0A%09%09out.alloc(in1.rows() * in2.rows(), in1.cols() * in2.cols());%0A%09%09for (long r1 = 0; r1 %3C in1.rows(); r1++) %7B%0A%09%09%09for (long c1 = 0; c1 %3C in1.cols(); c1++) %7B%0A%09%09%09%09in1.v(r1, c1, tmp1);%0A%09%09%09%09for (long r2 = 0; r2 %3C in2.rows(); r2++) %7B%0A%09%09%09%09%09for (long c2 = 0; c2 %3C in2.cols(); c2++) %7B%0A%09%09%09%09%09%09in2.v(r2, c2, tmp2);%0A%09%09%09%09%09%09group.multiply(tmp1, tmp2, tmp3);%0A%09%09%09%09%09%09out.setV(r1*in2.rows(), c1*in2.cols(), tmp3);%0A%09%09%09%09%09%7D%0A%09%09%09%09%7D%0A%09%09%09%7D%0A%09%09%7D%0A%09%7D%0A%7D%0A

Define a matrix direct product algorithm