Comment #7 on issue 3493 by [email protected]: bad solution for system of
linear equations
http://code.google.com/p/sympy/issues/detail?id=3493
I have tried the same using solve_linear_system. I'm not sure in the
underlining math, but it seems we have a problem here as well.
self.matrix
[ 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0.555560404773211, 0, 0,
0, 0, 0, 0, 0, 0,
0.511967099532406, 0, 0, 0, 0,
0, -0.384615371592977, 1.0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 1.0, 0,
0.831476179242813, 0, 0,
0, 0, 0, 0, 0, 0,
-0.859005057608147, -1.0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0.923076928502926, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0]
[ 0.384615371592977, -0.94537299195069, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, -1.0,
-0.555560404773211, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1.0, 0,
0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
-0.831476179242813, 0, 0,
0, 0, 0, 0, -1.0,
0, 0, 0, 0, 0,
0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1.0,
0, 0]
[ 0.923076928502926, 0.325990653378589, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 1.0, 0]
[ 0, 0, 0, 0,
0, 0, 1.0, 0, 0.94537299195069, 0,
0.0967822504713366, 0, 0, 0, 0,
0, 0, 0.107628721224152, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0.384615371592977, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, 0, -1.0,
-0.967822504713366, 0, 0, 0, 0,
1.0, 0, 0.960048201530856, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, -0.325990653378589, 0,
-0.232277410361082, 0, 0, 0, 0,
0, 0, -0.25830894120224, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
-0.923076928502926, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, -1.0, -0.384615371592977, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, -0.665816322167291,
0, 0, 0, 0, 0,
0, 0, 0, 0, -0.62278495383676,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0.746115691522181,
1.0, 0, 0, 0, 0,
0, 0, 0, -1.0, -0.782393060599686,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0,
0, 0, 0, -0.923076928502926, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0]
[ 0, 0, 0, -1.0,
-0.384615371592977, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, 0, 0,
0, 0, 0, -0.0967822504713
ss
[Edge, Edge001, Edge002, Edge003, Edge004, Edge005, Edge006, Edge007,
Edge008, Edge009, Edge010, Edge011, Edge012, Edge013, Edge014, Edge015,
Edge016, Edge017, Edge018, Edge019, Edge021, Edge022, Edge024, Edge025,
Edge026, Edge028, Edge029, Edge030, Edge031, Edge032, Edge033, Edge066,
Edge067, Edge068, Edge069, Edge070, Edge071, Edge072, Edge073, Edge075,
Edge076, Edge077, Edge078, Edge079, Edge080, Edge081, Edge082, Edge084,
Edge085, Edge087, Edge089, Edge090, Edge091, Edge092, Edge094, Edge095,
Edge096, Point006_x, Point006_y, Point006_z]
s = solve_linear_system(self.matrix,*ss)
s
{Edge009: 5.95801239645908e+17, Edge085: 389.065695247966, Edge090:
-1844.94618007660, Edge014: -1.27815984767627e+18, Edge005:
5.85833635394108e+17, Edge075: -1.42108547152020e-14, Edge010:
8.22190540511026e+17, Edge069: 0, Edge076: -745.188476077054, Edge017:
8.28848496356923e+17, Edge002: -1478.17773586637, Edge013:
-127.126887077369, Edge: -2.31940557606891e+17, Edge011:
-2.06890979154912e+17, Edge031: -7.61511406030670e+17, Point006_y:
-6303.58481499278, Edge026: -7.95734508269023e+17, Edge019:
-7.95734508269019e+17, Edge078: -217.207408837306, Edge066:
308.288433291077, Edge071: 1006.28705076761, Edge021: 8.28848496356915e+17,
Point006_x: -790.327727054299, Edge095: -1088.00000000000, Edge018:
-7.98537340007404e+17, Edge084: 498.864693094590, Edge001:
6.56764159622644e+17, Edge068: 2.84217094304040e-14, Edge067:
-108.874076140770, Edge080: 776.199023016558, Edge072: -2293.75210348430,
Edge030: -7.95734508269021e+17, Edge094: 1568.00000000000, Edge092:
-1833.29448392453, Edge087: 2776.19832038824, Edge003: 2.55804330630912,
Edge089: -3628.82156730837, Edge029: 1.06275946660744e+18, Edge012:
-7.95734508269018e+17, Edge007: 0, Edge096: 0, Edge022:
1.06275946660744e+18, Edge070: -768.620557314571, Edge015:
5.95801239645907e+17, Point006_z: -1051.30994714079, Edge004:
-2.06890979154912e+17, Edge006: 1.00593704068796e+18, Edge008:
-1.24259779501675e+18, Edge033: 1.34349974984468e+18, Edge024:
9.57062997990679, Edge028: 8.22190540511026e+17, Edge081:
-373.521821000162, Edge077: -599.974726334159, Edge073: 1933.42705667270,
Edge032: 0, Edge082: -312.257665334360, Edge025: -1.23719814824680e+18,
Edge079: 232.980343904048, Edge016: -2.31940557606892e+17, Edge091:
2657.96694715784}
(ss*self.matrix).subs(s)
[5.40769512773295e+17, 1.9097628957328e+17, -1568.0, 7.95734508269019e+17,
2.16084861050219e+18, 1.09871556578434e+18, 4.1013580104205e+17,
-3.81780256086895e+16, 1.3560609769535e+18, -2293.7521034843,
2.55957319279494e+17, -5.16730655521584e+17, 1.41827556565331e+18,
-7.95734508269021e+17, 5.89251684906775e+17, 0, 3.06051723587245e+17,
1.24597338145617e+18, 2.16755073572629e+17, 8.22190540511027e+17,
-5532.5993783771, 2.06890979154906e+17, 2.1392342581137e+18,
4.3074605149934e+17, -6.56764159622644e+17, 3.73900635933744e+17,
8.22190540511032e+17, -8.22190540511027e+17, -1.01433222091529e+18,
1.62312010844337e+18, -2.3194055760689e+17, -1.0251124971864e+18,
7.33275439469994e+17, -1.86129680661484e+18, 1848.0436743177,
2234.01187532457, 1355.1879610548, -1.30026929297642e+17, 5069.95042387254,
5.95801239645907e+17, -2.31940557606892e+17, -366.510318606426,
-5.36284050939776e+17, -7.95734508269022e+17, -7.27004632083905e+16,
-2648.39631717794, -1.12538919774218e+18, 432.930810742694,
2.05195901371184e+17, -2.9887594442021e+17, -75.3080893870295,
3.73900635933738e+17, -8.22190540511027e+17, 2785.09383423521,
5649.89577931083, -217.207408837306, 4299.35483172386, 2.55804330630912,
-2.06890979154912e+17, 5.85833635394108e+17, 8.28848496356915e+19]
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.
