예제 #1
0
        public Basis Execute(Matrix matrix)
        {
            var basis = new Basis();
            var row = matrix.GetRow(0);
            for (int i = 0; i < row.rows; i++)
            {
                if (row[i] == 0)
                {
                    var v = new Vector(row.rows);
                    v[i] = 1;
                    basis.Add(v);
                }
                else
                {
                    bool sign = row[i] > 0;
                    for (int j = i + 1; j < row.rows; j++)
                    {
                        if (row[j] != 0 && row[j] > 0 ^ sign)
                        {
                            var lcm = MathExtensions.LCM(Math.Abs((int)row[i]), Math.Abs((int)row[j]));
                            var v = new Vector(row.rows);
                            v[i] = (lcm / Math.Abs((int)row[i]));
                            v[j] = (lcm / Math.Abs((int)row[j]));
                            basis.Add(v);
                        }
                    }
                }
            }

            return ExcecuteRound(matrix, basis, 1);
        }
예제 #2
0
 private Basis Compose(Basis basis, Vector items)
 {
     var newBasis = new Basis();
     for (int i = 0; i < items.rows; i++)
     {
         if (items[i] == 0)
             newBasis.Add(basis[i]);
         else
         {
             bool sign = items[i] > 0;
             for (int j = i + 1; j < items.rows; j++)
             {
                 if (items[j] != 0 && items[j] > 0 ^ sign)
                 {
                     var lcm = MathExtensions.LCM(Math.Abs((int)items[i]), Math.Abs((int)items[j]));
                     newBasis.Add((lcm / Math.Abs((int)items[i])) * basis[i] + (lcm / Math.Abs((int)items[j])) * basis[j]);
                 }
             }
         }
     }
     return newBasis;
 }