Do the linear algebra in scipy, using my own code or packages like pyamg.Get the matrix out of FEniCS and into scipy.There is no built in function to export a sparse matrix to the scipy format (.toarray() converts to a dense matrix!) you have to specify a backend type and then write your own exporting function specifically for it.įor someone who just wants to specify a PDE and get a solution, maybe this is fine, but for someone like you or me who wants to experiment with building different preconditioners, solvers, investigate spectral properties of the matrices, etc, this is frustrating. The available linear algebra backend formats change from version to version, breaking old code. The particular details of matrix storage and matrix elements are considered an implementation detail that only developers should be concerned with. As far as I can tell this is a design decision, as they are trying to create an all-encompassing closed pipeline for discretizing and solving PDES. FEniCS tends to hide the details about the actual matrices it builds, and prevent easy manipulation of them.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |