Hacker Newsnew | past | comments | ask | show | jobs | submitlogin

The results are approximate (but that's always true), and, IIUC, the quality of the approximation depends on how well the matrix can be factored into a local component (with values inversely proportional to |(i - j)|) and a sparse global component.

Like petters writes above, "matrices of the form s^2 I + B, where B(i, j) depends on the distance between points i, j." The GMRES + FMM approach is more general, but FMM is less and less Fast as we move away from the structure above towards general matrices.

This structure makes sense because it reflects physical reality… useful in when solving physics problems (: Less so when the matrix corresponds to logical/business constraints over 0/1 variables which may not be so local anymore.



Guidelines | FAQ | Lists | API | Security | Legal | Apply to YC | Contact

Search: