Problems with linsolve: quotient by polynomial of higher degree
Subject: Problems with linsolve: quotient by polynomial of higher degree
From: Stefano Ferri
Date: Sat, 26 May 2012 21:18:42 +0200
2012/5/26 andre maute <andre.maute at gmx.de>
>
> Have you tried to keep the nonzero coefficients abstract,
> solve that system and substitute then the nonzero coefficients
> into the solution of the abstract system?
>
Yes, I've tried but computationally this is really a heavy task. Even
though my coefficient matrices (not reported in my example) are always
symmetric and sparse, and I only need to substitute coefficients that are
not 0 or 1, solving a system in such a way requires a lot of time and
results are very complicated, even for small linear systems. Therefore I
cannot use this method. By the way, even using a working method as Gauss'
elimination results are very complicated, but they are produced...
Anyway, the point is about linsolve... I am curious about its behaviour,
and if this is a bug.
Stefano