Christian Clason
Christian Clason
Julia 1.0, FEniCS.jl master, miniconda 4.5.10 and FEniCS 2018.1 works for me (linux x86_64). Are you sure you have set `conda activate fenicsproject`? If I forget that, I get exactly...
Aha, found the culprit! Whenever PyCall rebuilds (e.g., after an update), it resets to the default `python`, which is the system python unless you have activated a conda environment. So...
@aterenin Have you installed `mshr` by any chance? I can reproduce the problem if I install `mshr`. Curiously, it will say ``` The following packages will be UPDATED: fenics: 2016.1.0-py35_1...
@aterenin Since `FEniCS.jl` just wraps the `PyCall` (and, apparently installs FEniCS via `Conda.jl`), it will run into the same issues. Here's what (sort of) works (the order of the steps...
Exactly, the Julia equivalent to the last line is `sparse(row,col,val)`. The trick is getting at the Eigen backend data using PyCall (note the updated snippet, I forgot that this doesn't...
For PETSc matrices (assuming you do not want to (provide the option to) switch the backend to Eigen nor make use of the C++ PETSc interface), you could use a...
@MiroK I was not aware of the `mat()` function! Yes, if `A` is a `dolfin.cpp.la.Matrix` with PETSc backend, ```python row,col,val = as_backend_type(A).mat().getValuesCSR() A_csr = sp.csr_matrix((val,col,row)) ``` works nicely in Python....
Ah, I see the problem -- unlike the Eigen `data()`, PETSc's `getValuesCSR` does not return the `(row,col,val)` structure but the actual `ind,indptr,val` structure, which Julia's `sparse` doesn't understand (yet --...
I stand corrected; it's not in the docs, but `SparseMatrixCSC` *does* have a constructor for the standard structure. So you don't need SciPy at all and can instead just do...
How do you insert these characters?