public static string FormatNarrowPowers(this Rational r, Rational[] narrows = null) { if (narrows == null) { narrows = NarrowUtils.GetDefault(r.GetInvolvedPowerCount()); } Pow[] pows = GetNarrowPowers(r, narrows); if (pows == null) { return(null); //!!! where invalid narrowPrimes from? } return(Powers.ToString(pows, "|}")); }
/* * public static int[] FindCoordinates(Rational[] basis, Rational vector, int vectorLength) * { * int basisSize = basis.Length; * * // get prime powers * int[][] b = new int[basisSize][]; * for (int i = 0; i < basisSize; ++i) { * b[i] = basis[i].GetPrimePowers(); * } * int[] v = vector.GetPrimePowers(); * * // * int[] coordinates = FindCoordinates(b, v, vectorLength); * #if DEBUG * // check result * if (coordinates != null) { * Rational r = Rational.One; * for (int i = 0; i < coordinates.Length; ++i) { * r *= basis[i].Pow(coordinates[i]); * } * if (!r.Equals(vector)) throw new Exception( * String.Format("FindCoordinates failed: {0} != {1}", r, vector) * ); * } #endif * * return coordinates; * } */ #region Tests #if DEBUG private static void CheckVector(Rational[] basis, Rational vector, int vectorLength, bool addStandard = false) { // add standard basis int basisSize = basis.Length; if (addStandard) { Array.Resize(ref basis, basisSize + vectorLength); for (int i = 0; i < vectorLength; ++i) { basis[basisSize + i] = Rational.Prime(i); } } // for (int i = 0; i < basis.Length; ++i) { Debug.WriteLine(String.Format("Basis {0}. {1,-15} {2}", i, basis[i].FormatFraction(), basis[i].FormatMonzo())); } Debug.WriteLine(String.Format("Vector {0,-15} {1}", vector.FormatFraction(), vector.FormatMonzo())); // int[] coords = FindCoordinates(basis, vector, vectorLength); if (coords == null) { Debug.WriteLine("Invalid Basis or Out of basis subspace"); } else { Debug.WriteLine("Coordinates: " + Powers.ToString(coords, "()")); // Debug.WriteLine(vector.FormatFraction() + " = "); for (int i = 0; i < basisSize; ++i) { int e = Powers.SafeAt(coords, i); Debug.Print(" * ({0})^{1}", basis[i].FormatFraction(), e); } if (addStandard) { Rational r = new Rational(coords.Skip(basisSize).ToArray()); Debug.Print(" * {0} {1} {2}", r.FormatFraction(), r.FormatMonzo(), r.FormatNarrowPowers()); } } }
public string FormatMonzo() { return(Powers.ToString(pows)); }