//
// Source code recreated from a .class file by IntelliJ IDEA
// (powered by FernFlower decompiler)
//

package com.jisuan.Jama.examples;

//import Jama.EigenvalueDecomposition;
//import Jama.LUDecomposition;
//import Jama.Matrix;
//import Jama.QRDecomposition;
import com.jisuan.Jama.EigenvalueDecomposition;
import com.jisuan.Jama.LUDecomposition;
import com.jisuan.Jama.Matrix;
import com.jisuan.Jama.QRDecomposition;

import java.text.DecimalFormat;
import java.util.Date;

public class MagicSquareExample {
    public MagicSquareExample() {
    }

    public static Matrix magic(int var0) {
        double[][] var1 = new double[var0][var0];
        int var2;
        int var3;
        int var5;
        if (var0 % 2 == 1) {
            var2 = (var0 + 1) / 2;
            var3 = var0 + 1;

            for(int var4 = 0; var4 < var0; ++var4) {
                for(var5 = 0; var5 < var0; ++var5) {
                    var1[var5][var4] = (double)(var0 * ((var5 + var4 + var2) % var0) + (var5 + 2 * var4 + var3) % var0 + 1);
                }
            }
        } else if (var0 % 4 == 0) {
            for(var2 = 0; var2 < var0; ++var2) {
                for(var3 = 0; var3 < var0; ++var3) {
                    if ((var3 + 1) / 2 % 2 == (var2 + 1) / 2 % 2) {
                        var1[var3][var2] = (double)(var0 * var0 - var0 * var3 - var2);
                    } else {
                        var1[var3][var2] = (double)(var0 * var3 + var2 + 1);
                    }
                }
            }
        } else {
            var2 = var0 / 2;
            var3 = (var0 - 2) / 4;
            Matrix var9 = magic(var2);

            int var6;
            double var7;
            for(var5 = 0; var5 < var2; ++var5) {
                for(var6 = 0; var6 < var2; ++var6) {
                    var7 = var9.get(var6, var5);
                    var1[var6][var5] = var7;
                    var1[var6][var5 + var2] = var7 + (double)(2 * var2 * var2);
                    var1[var6 + var2][var5] = var7 + (double)(3 * var2 * var2);
                    var1[var6 + var2][var5 + var2] = var7 + (double)(var2 * var2);
                }
            }

            for(var5 = 0; var5 < var2; ++var5) {
                for(var6 = 0; var6 < var3; ++var6) {
                    var7 = var1[var5][var6];
                    var1[var5][var6] = var1[var5 + var2][var6];
                    var1[var5 + var2][var6] = var7;
                }

                for(var6 = var0 - var3 + 1; var6 < var0; ++var6) {
                    var7 = var1[var5][var6];
                    var1[var5][var6] = var1[var5 + var2][var6];
                    var1[var5 + var2][var6] = var7;
                }
            }

            double var10 = var1[var3][0];
            var1[var3][0] = var1[var3 + var2][0];
            var1[var3 + var2][0] = var10;
            var10 = var1[var3][var3];
            var1[var3][var3] = var1[var3 + var2][var3];
            var1[var3 + var2][var3] = var10;
        }

        return new Matrix(var1);
    }

    private static void print(String var0) {
        System.out.print(var0);
    }

    public static String fixedWidthDoubletoString(double var0, int var2, int var3) {
        DecimalFormat var4 = new DecimalFormat();
        var4.setMaximumFractionDigits(var3);
        var4.setMinimumFractionDigits(var3);
        var4.setGroupingUsed(false);

        String var5;
        for(var5 = var4.format(var0); var5.length() < var2; var5 = " " + var5) {
        }

        return var5;
    }

    public static String fixedWidthIntegertoString(int var0, int var1) {
        String var2;
        for(var2 = Integer.toString(var0); var2.length() < var1; var2 = " " + var2) {
        }

        return var2;
    }

    public static void main(String[] var0) {
        print("\n    Test of Matrix Class, using magic squares.\n");
        print("    See MagicSquareExample.main() for an explanation.\n");
        print("\n      n     trace       max_eig   rank        cond      lu_res      qr_res\n\n");
        Date var1 = new Date();
        double var2 = Math.pow(2.0, -52.0);

        for(int var4 = 3; var4 <= 32; ++var4) {
            print(fixedWidthIntegertoString(var4, 7));
            Matrix var5 = magic(var4);
            int var6 = (int)var5.trace();
            print(fixedWidthIntegertoString(var6, 10));
            EigenvalueDecomposition var7 = new EigenvalueDecomposition(var5.plus(var5.transpose()).times(0.5));
            double[] var8 = var7.getRealEigenvalues();
            print(fixedWidthDoubletoString(var8[var4 - 1], 14, 3));
            int var9 = var5.rank();
            print(fixedWidthIntegertoString(var9, 7));
            double var10 = var5.cond();
            print(var10 < 1.0 / var2 ? fixedWidthDoubletoString(var10, 12, 3) : "         Inf");
            LUDecomposition var12 = new LUDecomposition(var5);
            Matrix var13 = var12.getL();
            Matrix var14 = var12.getU();
            int[] var15 = var12.getPivot();
            Matrix var16 = var13.times(var14).minus(var5.getMatrix(var15, 0, var4 - 1));
            double var17 = var16.norm1() / ((double)var4 * var2);
            print(fixedWidthDoubletoString(var17, 12, 3));
            QRDecomposition var19 = new QRDecomposition(var5);
            Matrix var20 = var19.getQ();
            var16 = var19.getR();
            var16 = var20.times(var16).minus(var5);
            var17 = var16.norm1() / ((double)var4 * var2);
            print(fixedWidthDoubletoString(var17, 12, 3));
            print("\n");
        }

        Date var21 = new Date();
        double var22 = (double)(var21.getTime() - var1.getTime()) / 1000.0;
        print("\nElapsed Time = " + fixedWidthDoubletoString(var22, 12, 3) + " seconds\n");
        print("Adios\n");
    }
}