Changeset 5f4e0e in git

Timestamp:

Mar 10, 2014, 3:28:59 PM (10 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

0aa5d5341e6004e6426fabd0fb0a1f19f6bd7eb4

Parents:

d21c2cb52c17b8ca70c95d559d91d0d074c00a7e

git-author:

Martin Lee <martinlee84@web.de>2014-03-10 15:28:59+01:00

git-committer:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2014-03-17 18:31:54+01:00

Message:

chg: added 'subst' to substitute alg. variables in f by their
 image in the absolute extension

File:

• factory/libfac/charset/alg_factor.cc (modified) (1 diff)

factory/libfac/charset/alg_factor.cc

@@ -696 +696 @@
   //assert(f.inCoeffDomain());
   return f;
+}
+
+CanonicalForm
+subst (const CanonicalForm& f, const CFList& a, const CFList& b,
+       const CanonicalForm& Rstar, bool isFunctionField)
+{
+  if (isFunctionField)
+    ASSERT (2*a.length() == b.length(), "wrong length of lists");
+  else
+    ASSERT (a.length() == b.length(), "lists of equal length expected");
+  CFListIterator j= b;
+  CanonicalForm result= f, tmp, powj;
+  for (CFListIterator i= a; i.hasItem() && j.hasItem(); i++, j++)
+  {
+    if (!isFunctionField)
+      result= result (j.getItem(), i.getItem().mvar());
+    else
+    {
+      tmp= j.getItem();
+      j++;
+      powj= power (j.getItem(), degree (result, i.getItem().mvar()));
+      result= evaluate (result, tmp, j.getItem(), powj, i.getItem().mvar());
+
+      if (fdivides (powj, result, tmp))
+        result= tmp;
+
+      result /= vcontent (result, Variable (i.getItem().level() + 1));
+    }
+  }
+  result= reduce (result, Rstar);
+  result /= vcontent (result, Variable (Rstar.level() + 1));
+  return result;
 }