Comment #4 on issue 3493 by [email protected]: bad solution for system of linear equations
http://code.google.com/p/sympy/issues/detail?id=3493
syms = 'Edge082', 'Edge080', 'Edge081', 'Edge068', 'Edge069', 'Edge084', 'Edge085', 'Edge066', 'Edge089', 'Edge087', 'Edge078', 'Edge008', 'Edge033', 'Edge032', 'Edge031', 'Edge030', 'Edge019', 'Edge018', 'Edge015', 'Edge014', 'Edge017', 'Edge016', 'Edge011', 'Edge010', 'Edge013', 'Edge012', 'Edge095', 'Edge094', 'Edge096', 'Edge091', 'Edge090', 'Edge079', 'Edge092', 'Edge077', 'Edge076', 'Edge075', 'Edge067', 'Edge072', 'Edge071', 'Point006_x', 'Edge002', 'Edge', 'Point006_y', 'Edge029', 'Edge073', 'Edge070', 'Point024_z', 'Point024_y', 'Point024_x', 'Edge021', 'Edge022', 'Edge009', 'Edge024', 'Edge025', 'Edge026', 'Edge028', 'Edge003', 'Point006_z', 'Edge001', 'Edge006', 'Edge007', 'Edge004', 'Edge005' Edge082, Edge080, Edge081, Edge068, Edge069, Edge084, Edge085, Edge066, Edge089, Edge087, Edge078, Edge008, Edge033, Edge032, Edge031, Edge030, Edge019, Edge018, Edge015, Edge014, Edge017, Edge016, Edge011, Edge010, Edge013, Edge012, Edge095, Edge094, Edge096, Edge091, Edge090, Edge079, Edge092, Edge077, Edge076, Edge075, Edge067, Edge072, Edge071, Point006_x, Edge002, Edge, Point006_y, Edge029, Edge073, Edge070, Point024_z, Point024_y, Point024_x, Edge021, Edge022, Edge009, Edge024, Edge025, Edge026, Edge028, Edge003, Point006_z, Edge001, Edge006, Edge007, Edge004, Edge005 = symbols(syms) eqs = [-0.0967822504713366*Edge010 - 0.384615371592977*Edge011 + 1.0*Edge024 - 0.0967822504713366*Edge072, 0.967822504713366*Edge010 + 1.0*Edge026 - 0.967822504713366*Edge072 - 1.0*Edge089, 0.232277410361082*Edge010 + 0.923076928502926*Edge011 + 0.232277410361082*Edge072, 0.384615371592977*Edge069 - 0.94537299195069*Edge071 + 0.107628721224152*Edge085 + 0.0967822504713366*Edge090 - 1.0*Edge095, -1.0*Edge084 - 0.960048201530856*Edge085 + 0.967822504713366*Edge090 + 1.0*Edge091, 0.923076928502926*Edge069 + 0.325990653378589*Edge071 + 0.25830894120224*Edge085 + 0.232277410361082*Edge090, -0.384615371592977*Edge - 1.0*Edge002 - 0.107628721224152*Edge021, -1.0*Edge019 - 0.960048201530856*Edge021, -0.923076928502926*Edge - 0.25830894120224*Edge021, -1.0*Edge003 - 0.384615371592977*Edge004 - 0.0967822504713366*Edge028 - 0.0967822504713366*Edge090, 0.967822504713366*Edge028 + 1.0*Edge030 - 0.967822504713366*Edge090 - 1.0*Edge092, -0.923076928502926*Edge004 - 0.232277410361082*Edge028 - 0.232277410361082*Edge090, 0.94537299195069*Edge066 + 1.0*Edge068 + 0.384615371592977*Edge078 + 0.665816322167291*Edge082, 1.0*Edge079 + 0.746115691522181*Edge082, -0.325990653378589*Edge066 - 0.923076928502926*Edge078 - 100.0, 1.0*Edge075 - 0.384615371592977*Edge078 - 0.107628721224152*Edge080, 1.0*Edge076 + 0.960048201530856*Edge080, 0.923076928502926*Edge078 + 0.25830894120224*Edge080, 1.0*Edge006 + 0.94537299195069*Edge008 + 0.0967822504713366*Edge010 + 0.107628721224152*Edge017 + 0.384615371592977*Edge032, -1.0*Edge009 - 0.967822504713366*Edge010 + 1.0*Edge015 + 0.960048201530856*Edge017, -0.325990653378589*Edge008 - 0.232277410361082*Edge010 - 0.25830894120224*Edge017 - 0.923076928502926*Edge032, -1.0*Edge006 - 0.384615371592977*Edge007 - 0.665816322167291*Edge018 - 0.62278495383676*Edge031, 0.746115691522181*Edge018 + 1.0*Edge019 - 1.0*Edge030 - 0.782393060599686*Edge031, -0.923076928502926*Edge007, 0.384615371592977*Edge - 0.94537299195069*Edge001 - 1.0*Edge013 - 0.555560404773211*Edge014 + Point006_x, -0.831476179242813*Edge014 - 1.0*Edge022 + Point006_y, 0.923076928502926*Edge + 0.325990653378589*Edge001 + Point006_z, 1.0*Edge013 - 0.384615371592977*Edge016 - 0.107628721224152*Edge017, -1.0*Edge012 - 0.960048201530856*Edge017, 0.923076928502926*Edge016 + 0.25830894120224*Edge017, 0.555560404773211*Edge014 + 0.511967099532406*Edge025 - 0.384615371592977*Edge032 + 1.0*Edge033, 1.0*Edge012 + 0.831476179242813*Edge014 - 0.859005057608147*Edge025 - 1.0*Edge026, 0.923076928502926*Edge032, 0.384615371592977*Edge007 - 0.94537299195069*Edge008 + 0.107628721224152*Edge021 + 0.0967822504713366*Edge028 - 1.0*Edge033, 0.960048201530856*Edge021 + 1.0*Edge022 - 0.967822504713366*Edge028 - 1.0*Edge029, 0.923076928502926*Edge007 + 0.325990653378589*Edge008 + 0.25830894120224*Edge021 + 0.232277410361082*Edge028, 0.94537299195069*Edge001 + 1.0*Edge002 + 0.384615371592977*Edge016 + 0.665816322167291*Edge018 + Point024_x, -1.0*Edge015 - 0.746115691522181*Edge018 + Point024_y, -0.325990653378589*Edge001 - 0.923076928502926*Edge016 + Point024_z, 1.0*Edge003 + 0.94537299195069*Edge005 + 0.384615371592977*Edge011 + 0.62278495383676*Edge031 + 0.62278495383676*Edge094, 1.0*Edge009 + 0.782393060599686*Edge031 - 1.0*Edge073 - 0.782393060599686*Edge094, -0.325990653378589*Edge005 - 0.923076928502926*Edge011, -0.384615371592977*Edge069 - 1.0*Edge070 - 0.665816322167291*Edge082 - 0.62278495383676*Edge094, -1.0*Edge081 - 0.746115691522181*Edge082 + 1.0*Edge092 + 0.782393060599686*Edge094, -0.923076928502926*Edge069, -0.94537299195069*Edge066 + 0.384615371592977*Edge067 - 1.0*Edge075 - 0.555560404773211*Edge077, 0.831476179242813*Edge077 + 1.0*Edge084, 0.325990653378589*Edge066 + 0.923076928502926*Edge067, 0.384615371592977*Edge004 - 0.94537299195069*Edge005 - 1.0*Edge024 - 0.511967099532406*Edge025 - 0.511967099532406*Edge087, 0.859005057608147*Edge025 + 1.0*Edge029 - 0.859005057608147*Edge087 - 1.0*Edge091, 0.923076928502926*Edge004 + 0.325990653378589*Edge005, 1.0*Edge070 + 0.94537299195069*Edge071 + 0.0967822504713366*Edge072 + 0.107628721224152*Edge080 + 0.384615371592977*Edge096, 0.967822504713366*Edge072 + 1.0*Edge073 - 1.0*Edge079 - 0.960048201530856*Edge080, -0.325990653378589*Edge071 - 0.232277410361082*Edge072 - 0.25830894120224*Edge080 - 0.923076928502926*Edge096, -0.384615371592977*Edge067 - 1.0*Edge068 - 0.107628721224152*Edge085, 1.0*Edge081 + 0.960048201530856*Edge085, -0.923076928502926*Edge067 - 0.25830894120224*Edge085, 0.555560404773211*Edge077 + 0.511967099532406*Edge087 + 1.0*Edge095 - 0.384615371592977*Edge096, -1.0*Edge076 - 0.831476179242813*Edge077 + 0.859005057608147*Edge087 + 1.0*Edge089, 0.923076928502926*Edge096, Edge094, Edge095, Point024_z]
s=solve(eqs) for i in eqs: print i.subs(s) -- 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.
