Re: [petsc-users] MatPtAP

2018-06-04 Thread Samuel Lanthaler
in the profiler with your matrices and see if there is any way to speed up the computation. Barry > On Jun 1, 2018, at 11:07 AM, Samuel Lanthaler mailto:s.lantha...@gmail.com>> wrote: > > On 06/01/2018 03:42 PM, Matthew Knepley wrote: >> On

Re: [petsc-users] MatPtAP

2018-06-01 Thread Samuel Lanthaler
me attach the output of -log_view. Please send the output of a run with -log_view so we can see where the time is being spent. Barry On Jun 1, 2018, at 8:21 AM, Samuel Lanthaler wrote: Hi, I was wondering what the most efficient way to use MatPtAP would be in the following

Re: [petsc-users] MatPtAP

2018-06-01 Thread Samuel Lanthaler
On 06/01/2018 03:42 PM, Matthew Knepley wrote: On Fri, Jun 1, 2018 at 9:21 AM, Samuel Lanthaler mailto:s.lantha...@gmail.com>> wrote: Hi, I was wondering what the most efficient way to use MatPtAP would be in the following situation: I am discretizing a PDE

[petsc-users] MatPtAP

2018-06-01 Thread Samuel Lanthaler
Hi, I was wondering what the most efficient way to use MatPtAP would be in the following situation: I am discretizing a PDE system. The discretization yields a matrix A that has a band structure (with k upper and lower bands, say). In order to implement the boundary conditions, I use a

Re: [petsc-users] configuration option for PETSc on cluster

2018-04-25 Thread Samuel Lanthaler
or alternatively - you can disable sharedlibraries --with-shared-libraries=0 Satish On Tue, 24 Apr 2018, Samuel Lanthaler wrote: Dear Satish, As you suggested, the configure appears to work with ./configure --with-fc=ftn --with-cc=cc --with-cxx=CC --with-debugging=0 --with-scalar-type=complex

Re: [petsc-users] configuration option for PETSc on cluster

2018-04-24 Thread Samuel Lanthaler
s an extra option 'LDFLAGS=-dynamic' Satish On Tue, 24 Apr 2018, Samuel Lanthaler wrote: Dear Satish, I get the following when doing the ls: [lanthale@daint103 mkl]$ ls /opt/intel/compilers_and_libraries_2017.4.196/linux/mkl/lib/intel64/ libmkl_ao_worker.so libmkl_blacs_sgim

Re: [petsc-users] configuration option for PETSc on cluster

2018-04-24 Thread Samuel Lanthaler
- these should be the ones that get picked up. I do not know why 'cc' is picking up static libraries instead of the shared ones. Perhaps it needs an extra option 'LDFLAGS=-dynamic' Satish On Tue, 24 Apr 2018, Samuel Lanthaler wrote: Dear Satish, I get the following when doing the ls: [lanthale

Re: [petsc-users] configuration option for PETSc on cluster

2018-04-24 Thread Samuel Lanthaler
ish Balay wrote: On Tue, 24 Apr 2018, Samuel Lanthaler wrote: Hi there, I was wondering if you could help me with the correct configuration of PETSc-dev version on a cluster (https://www.cscs.ch/)? I'm not sure which information would be useful to you, but basically the problem seems to b

Re: [petsc-users] VecGetArrayReadF90 with gfortran

2018-03-25 Thread Samuel Lanthaler
) 6.2.0 Copyright (C) 2016 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. Satish On Fri, 23 Mar 2018, Samuel Lanthaler wrote: Thank you very much for your effort

Re: [petsc-users] VecGetArrayReadF90 with gfortran

2018-03-23 Thread Samuel Lanthaler
, 2018, at 9:10 AM, Samuel Lanthaler <s.lantha...@gmail.com> wrote: Hi all, I am having trouble using the function VecGetArrayReadF90 under gfortran. I have created a minimal example and put it in the attachment. Basically, it appears that the input pointer is empty upon retur

[petsc-users] VecGetArrayReadF90 with gfortran

2018-03-23 Thread Samuel Lanthaler
Hi all, I am having trouble using the function VecGetArrayReadF90 under gfortran. I have created a minimal example and put it in the attachment. Basically, it appears that the input pointer is empty upon return from VecGetArrayReadF90. The code runs fine, when compiled with ifort, so I don't know

Re: [petsc-users] make: *** No rule to make target `chk_makej', needed by `all'. Stop.

2018-03-08 Thread Samuel Lanthaler
two days for us to re-synchronize. Jose El 9 mar 2018, a las 3:59, Samuel Lanthaler <s.lantha...@gmail.com> escribió: Hi, I've just been trying to compile SLEPC on a cluster, and am getting the following error message after I type "make": ./configure Checking environment...

Re: [petsc-users] make: *** No rule to make target `chk_makej', needed by `all'. Stop.

2018-03-08 Thread Samuel Lanthaler
and try again. When a change in petsc-master affects SLEPc, you have to allow one or two days for us to re-synchronize. Jose El 9 mar 2018, a las 3:59, Samuel Lanthaler <s.lantha...@gmail.com> escribió: Hi, I've just been trying to compile SLEPC on a cluster, and am getting the following

[petsc-users] make: *** No rule to make target `chk_makej', needed by `all'. Stop.

2018-03-08 Thread Samuel Lanthaler
Hi, I've just been trying to compile SLEPC on a cluster, and am getting the following error message after I type "make": > ./configure Checking environment... done Checking PETSc installation... done Checking LAPACK library... done Writing various configuration files... done Generating

Re: [petsc-users] Using SLEPC: Fortran call to NEPNLEIGSSetSingularitiesFunction (similar to SNESSetFunction?)

2018-01-11 Thread Samuel Lanthaler
Amazing, thank you! Samuel On 01/11/2018 02:48 PM, Jose E. Roman wrote: The Fortran stub is now in master. I have also completed the example. https://bitbucket.org/slepc/slepc/commits/acaefe873 Jose El 10 ene 2018, a las 17:29, Samuel Lanthaler <s.lantha...@gmail.com> escribió: Ah,

Re: [petsc-users] Using SLEPC: Fortran call to NEPNLEIGSSetSingularitiesFunction (similar to SNESSetFunction?)

2018-01-10 Thread Samuel Lanthaler
Ah, I see. Great, thank you, Jose! Cheers, Samuel On 01/10/2018 05:23 PM, Jose E. Roman wrote: This function is not implemented in Fortran. I will add it and let you know. Jose El 10 ene 2018, a las 17:11, Samuel Lanthaler <s.lantha...@gmail.com> escribió: Hi there, I'm stuck g

[petsc-users] Using SLEPC: Fortran call to NEPNLEIGSSetSingularitiesFunction (similar to SNESSetFunction?)

2018-01-10 Thread Samuel Lanthaler
Hi there, I'm stuck getting a call to the SLEPC routine "NEPNLEIGSSetSingularitiesFunction" to work from within Fortran. To be more precise: To get started with the solution of a non-linear eigenvalue-problem, I am first trying to write a Fortran version of

Re: [petsc-users] MatCreateShell, MatShellGetContext, MatShellSetContext in fortran

2017-12-13 Thread Samuel Lanthaler
e question, Barry > On Dec 12, 2017, at 3:52 AM, Samuel Lanthaler <s.lantha...@gmail.com> wrote: > > Let me also add a minimal example (relying only on petsc), which leads to the same error message on my machine. Again, what I'm trying to do is very simple: >•

Re: [petsc-users] MatCreateShell, MatShellGetContext, MatShellSetContext in fortran

2017-12-12 Thread Samuel Lanthaler
r Message -- [0]PETSC ERROR: Null argument, when expecting valid pointer [0]PETSC ERROR: Null Pointer: Parameter # 2 Again, thanks for your help! Cheers, Samuel On 12/11/2017 07:41 PM, Samuel Lanthaler wrote: Dear petsc-/slepc-users, I have been trying to understand matrix-free/shell

[petsc-users] MatCreateShell, MatShellGetContext, MatShellSetContext in fortran

2017-12-11 Thread Samuel Lanthaler
Dear petsc-/slepc-users, I have been trying to understand matrix-free/shell matrices in PETSc for eventual use in solving a non-linear eigenvalue problem using SLEPC. But I seem to be having trouble with calls to MatShellGetContext. As far as I understand, this function should initialize a

Re: [petsc-users] MatPtAP problem after subsequent call to MatDuplicate

2017-12-05 Thread Samuel Lanthaler
Thank you for your swift reply, Hong! Yes, what you said is exactly what I'm doing. No, I don't get an error when doing a MatView of matC. Samuel On 12/05/2017 04:38 PM, Hong wrote: Samuel: You try to do following: 1) Create A; 2) Create P; 3) C = PtAP: CALL

[petsc-users] MatPtAP problem after subsequent call to MatDuplicate

2017-12-05 Thread Samuel Lanthaler
Hi there, I am getting error messages after using MatPtAP to create a new matrix C = Pt*A*P and then trying to assign A=C. The following is a minimal example reproducing the problem: #include "slepc/finclude/slepc.h" USE slepcsys USE slepceps IMPLICIT NONE

Re: [petsc-users] MatZeroRowsColumns

2017-11-24 Thread Samuel Lanthaler
Perfect, PETSC_NULL_VEC works in fortran. Thank you very much for your help! :-) Cheers, Sam On 11/24/2017 03:18 PM, Matthew Knepley wrote: On Fri, Nov 24, 2017 at 6:42 AM, Samuel Lanthaler <s.lantha...@gmail.com <mailto:s.lantha...@gmail.com>> wrote: Ah, great. I didn'

Re: [petsc-users] MatZeroRowsColumns

2017-11-24 Thread Samuel Lanthaler
case from fortran? If not, then I'll just stick to creating vectors and distroying them afterwards, which is fine for me as well. Cheers, Sam On 11/24/2017 01:56 AM, Smith, Barry F. wrote: MatZeroRowsColumns() as two vector arguments you are missing On Nov 23, 2017, at 2:43 PM, Samuel

[petsc-users] MatZeroRowsColumns

2017-11-23 Thread Samuel Lanthaler
Hi there, I'm new to PETSc and have been trying to do some basic manipulation of matrices in fortran, but don't seem to be able to set a row/column to zero using the MatZeroRowsColumns command: The following is a small example program:| | |! initialize PETSc|| ||  CALL