Hi Jose, While I am running from petsc4py import PETSc
print(PETSc.ScalarType) <class 'numpy.float64'> only float64 I am getting not complex128.
Do I have to uninstall and then install it? Otherwise doing F.setScale(-1j*t) I am getting an error File "SLEPc/FN.pyx", line 204, in slepc4py.SLEPc.FN.setScale File "SLEPc/SLEPc.pyx", line 115, in slepc4py.SLEPc.asScalar TypeError: can't convert complex to float Sanku ________________________________ From: Jose E. Roman <[email protected]> Sent: Saturday, September 4, 2021 5:33 AM To: Paul, Sanku <[email protected]> Cc: PETSc <[email protected]> Subject: Re: [petsc-users] Matrix exponential No problem, do F.setScale(-1j*t) What do you mean "still facing problem with complex matrices"? Did you get an error during installation? Did it install for real scalars? Try setting PETSC_CONFIGURE_OPTIONS only, and unset PETSC_DIR PETSC_ARCH. Otherwise pip may not take into account --with-scalar-type=complex Jose
El 4 sept 2021, a las 2:11, Paul, Sanku <[email protected]> escribió:
Hi Jose,
I tried to install ml -* foss/2019b Python SciPy-bundle/2019.10-Python-3.7.4 virtualenv slepc4py cd slepc4py source bin/activate pip install numpy mpi4py cython export PETSC_CONFIGURE_OPTIONS="--with-scalar-type=complex" export PETSC_DIR=/path/to/petsc PETSC_ARCH=your-arch-name pip install petsc petsc4py pip install slepc slepc4py
But still facing problem with complex matrices. Please help me to fix this.
Thanks, Sanku
From: Paul, Sanku <[email protected]> Sent: Friday, September 3, 2021 7:39 PM To: Jose E. Roman <[email protected]> Cc: PETSc <[email protected]> Subject: Re: [petsc-users] Matrix exponential
Hi Jose,
I could now do matrix exponential but facing a problem with a complex matrix. In particular, I want to do \exp(-itH), where H is a Hamiltonian. How to implement this?
Thanks, Sanku From: Jose E. Roman <[email protected]> Sent: Friday, September 3, 2021 2:53 PM To: Paul, Sanku <[email protected]> Cc: PETSc <[email protected]> Subject: Re: [petsc-users] Matrix exponential
Please always reply to the list (Reply-All), not to myself.
You should be able to convert from a scipy sparse matrix to a PETSc matrix via PETSc.Mat().createAIJWithArrays(). Don't know how if there is any example in the petsc4py documentation.
Jose
El 3 sept 2021, a las 20:26, Paul, Sanku <[email protected]> escribió:
Dear Jose,
Thank you very much for your help. I have another question can we just simply pass a sparse.csr.matrix to A. For instance, if B is the sparse.csr.matrix can we do A=B.copy(). Or do I have to do it in a different way?
Best, Sanku From: Jose E. Roman <[email protected]> Sent: Friday, September 3, 2021 2:13 PM To: Paul, Sanku <[email protected]> Cc: [email protected] <[email protected]> Subject: Re: [petsc-users] Matrix exponential
You should either create the FN object and then
E.setFN(F)
or extract the FN object and assign to a variable
F = E.getFN()
You can see an example in $SLEPC_DIR/src/binding/slepc4py/demo/ex6.py
Jose
El 3 sept 2021, a las 19:53, Paul, Sanku <[email protected]> escribió:
Dear Sir/Ma'am,
I am trying to use SLEPc to calculate matrix exponential in my python code but I am not getting the correct result. I have attached the code. Could you let me know what I am doing wrong. This is my first time using SLEPc. So, I would like to ask you if you could send me a tutorial on matrix exponential using SLEPc in python code.
Best, Sanku <ex2.py>