{"action":"closed","issue":{"url":"https://api.github.com/repos/AtheMathmo/rulinalg/issues/74","repository_url":"https://api.github.com/repos/AtheMathmo/rulinalg","labels_url":"https://api.github.com/repos/AtheMathmo/rulinalg/issues/74/labels{/name}","comments_url":"https://api.github.com/repos/AtheMathmo/rulinalg/issues/74/comments","events_url":"https://api.github.com/repos/AtheMathmo/rulinalg/issues/74/events","html_url":"https://github.com/AtheMathmo/rulinalg/issues/74","id":184999281,"number":74,"title":"Erroneous singular matrix in .solve?","user":{"login":"bright-star","id":5308179,"avatar_url":"https://avatars.githubusercontent.com/u/5308179?v=3","gravatar_id":"","url":"https://api.github.com/users/bright-star","html_url":"https://github.com/bright-star","followers_url":"https://api.github.com/users/bright-star/followers","following_url":"https://api.github.com/users/bright-star/following{/other_user}","gists_url":"https://api.github.com/users/bright-star/gists{/gist_id}","starred_url":"https://api.github.com/users/bright-star/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/bright-star/subscriptions","organizations_url":"https://api.github.com/users/bright-star/orgs","repos_url":"https://api.github.com/users/bright-star/repos","events_url":"https://api.github.com/users/bright-star/events{/privacy}","received_events_url":"https://api.github.com/users/bright-star/received_events","type":"User","site_admin":false},"labels":[],"state":"closed","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":7,"created_at":"2016-10-25T01:54:41Z","updated_at":"2016-11-12T05:47:30Z","closed_at":"2016-11-12T05:47:30Z","body":"I'm verifying the results of the FD port I have written using rulinalg, which consists of the following two files:\n- https://gitlab.com/bright-star/electro/blob/master/examples/empty_rect.rs\n- https://gitlab.com/bright-star/electro/blob/master/src/lib.rs\n\nRunning that example `cargo run --example empty_rect` with nightly gives an LUP-related error:\n\n```\nthread 'main' panicked at 'called `Result::unwrap()` on an `Err` value: \nError { kind: DivByZero, error: StringError(\"Singular matrix found in LUP decomposition. A value in the diagonal of U == 0.0.\") }', ../src/libcore/result.rs:799\n```\n\nIt also prints out the matrix and vector it was trying to solve with, which I verified against the MATLAB source I am porting with the same params: https://github.com/bright-star/computational-em/blob/master/finite-difference/e1_uniform.m So I know that this a solvable system.\n\nFor my sanity's sake: these are identical matrices, right?\n\n```\n-4   1   0   1   0   0   0   0   0\n 1  -4   1   0   1   0   0   0   0\n 0   1  -4   0   0   1   0   0   0\n 1   0   0  -4   1   0   1   0   0\n 0   1   0   1  -4   1   0   1   0\n 0   0   1   0   1  -4   0   0   1\n 0   0   0   1   0   0  -4   1   0\n 0   0   0   0   1   0   1  -4   1\n 0   0   0   0   0   1   0   1  -4\n\n-4   1   0   1   0   0   0   0   0\n 1  -4   1   0   1   0   0   0   0\n 0   1  -4   0   0   1   0   0   0\n 1   0   0  -4   1   0   1   0   0\n 0   1   0   1  -4   1   0   1   0\n 0   0   1   0   1  -4   0   0   1\n 0   0   0   1   0   0  -4   1   0\n 0   0   0   0   1   0   1  -4   1\n 0   0   0   0   0   1   0   1  -4\n```\n"}}
{"actor":{"display_login":"bright-star"}}