예제 #1
0
        public static double normF(DMatrix5 M)
        {
            double scale = CommonOps_DDF5.elementMaxAbs(M);

            if (scale == 0.0)
            {
                return(0.0);
            }

            double a1 = M.a1 / scale, a2 = M.a2 / scale, a3 = M.a3 / scale, a4 = M.a4 / scale, a5 = M.a5 / scale;
            double sum = a1 * a1 + a2 * a2 + a3 * a3 + a4 * a4 + a5 * a5;

            return(scale * Math.Sqrt(sum));
        }
예제 #2
0
        public static double normF(DMatrix5x5 M)
        {
            double scale = CommonOps_DDF5.elementMaxAbs(M);

            if (scale == 0.0)
            {
                return(0.0);
            }

            double a11 = M.a11 / scale,
                   a12 = M.a12 / scale,
                   a13 = M.a13 / scale,
                   a14 = M.a14 / scale,
                   a15 = M.a15 / scale;
            double a21 = M.a21 / scale,
                   a22 = M.a22 / scale,
                   a23 = M.a23 / scale,
                   a24 = M.a24 / scale,
                   a25 = M.a25 / scale;
            double a31 = M.a31 / scale,
                   a32 = M.a32 / scale,
                   a33 = M.a33 / scale,
                   a34 = M.a34 / scale,
                   a35 = M.a35 / scale;
            double a41 = M.a41 / scale,
                   a42 = M.a42 / scale,
                   a43 = M.a43 / scale,
                   a44 = M.a44 / scale,
                   a45 = M.a45 / scale;
            double a51 = M.a51 / scale,
                   a52 = M.a52 / scale,
                   a53 = M.a53 / scale,
                   a54 = M.a54 / scale,
                   a55 = M.a55 / scale;

            double sum = 0;

            sum += a11 * a11 + a12 * a12 + a13 * a13 + a14 * a14 + a15 * a15;
            sum += a21 * a21 + a22 * a22 + a23 * a23 + a24 * a24 + a25 * a25;
            sum += a31 * a31 + a32 * a32 + a33 * a33 + a34 * a34 + a35 * a35;
            sum += a41 * a41 + a42 * a42 + a43 * a43 + a44 * a44 + a45 * a45;
            sum += a51 * a51 + a52 * a52 + a53 * a53 + a54 * a54 + a55 * a55;

            return(scale * Math.Sqrt(sum));
        }
예제 #3
0
        public static void normalizeF(DMatrix5 M)
        {
            double val = normF(M);

            CommonOps_DDF5.divide(M, val);
        }