gserdyuk/gnucap-a — Issue #14: linear solver fails if diagonal is zero
Author: gserdyuk
Created: 2011-11-13 23:31:19+00:00
State: open
Labels: low_priority

Lets take system:

[ gmin -1 ] * X = [ -1 ]
10 -10 0

if gmin = 0 this system has exact solution X= 1; 1

but current solver behaves wrong (gives X= {1.0} ) in wide range of gmin.

shortly :
if gmin <=1.e-17 solver returns X={0, 1}
if gmin = 1e-6 X={1,1} (correct)

pls see attached
a) circuit file with test ckt
b) test functiion for "boost test"

now can use current solver, keep in mid that gmin ~ 1.e-6

in future - substitute linear solver with another one.