Re: [petsc-dev] git help needed

2014-02-18 Thread Svetlana Tkachenko
You can `git reset --soft previous item in git log`, I think, switch to another branch, and commit twice (once to master and once to that branch). - Original message - From: Barry Smith bsm...@mcs.anl.gov To: petsc-dev petsc-dev@mcs.anl.gov, Karl Rupp r...@mcs.anl.gov Subject:

[petsc-dev] stensils, ex2f.F, openmp

2013-11-07 Thread Svetlana Tkachenko
ex2f.F involves a band matrix with 4 on main diagonal and -1s on the sub-diagonals. I have A_central matrix which stores the -4s, and A_left, A_right, Atop, A_bottom matrices which store the -1s (they actually are more complicated numbers than that). Would it be recommended to simply wrap the

Re: [petsc-dev] stencils, ex2f.F, openmp

2013-11-07 Thread Svetlana Tkachenko
On Fri, 8 Nov 2013, at 14:31, Gryllida wrote: On Fri, 8 Nov 2013, at 14:08, Gryllida wrote: On Fri, 8 Nov 2013, at 10:09, Svetlana Tkachenko wrote: ex2f.F involves a band matrix with 4 on main diagonal and -1s on the sub-diagonals. I have A_central matrix which stores the -4s

Re: [petsc-dev] stencils, ex2f.F, openmp

2013-11-07 Thread Svetlana Tkachenko
On Fri, 8 Nov 2013, at 15:17, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: The example has a line to compute matrix A from Ax = b, and a line to compute source term (vector b). call KSPSetComputeOperators(ksp,ComputeMatrix,ctx,ierr) call KSPSetComputeRHS

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
makefiles take care of it and are portable for different MPIs etc. Barry On Nov 6, 2013, at 12:25 AM, Svetlana Tkachenko svetlana.tkache...@fastmail.fm wrote: I have configured petsc-dev (downloaded it today) with these options, and a small example. It appears to fail to compile

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 9:22, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: Thanks. Could you please send an example makefile for this small example as src/ts/examples/tutorials/makefile is huge and I don't know which part of it manages the -I -L stuff

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 10:22, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: Right. I have spent half of an hour now trying to imagine what to do to link trying everything like a headless chicken and it did not work. You always have to send the error

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 12:34, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: On Thu, 7 Nov 2013, at 10:22, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: Right. I have spent half of an hour now trying to imagine what to do

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 12:56, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: ~/dev/test/petsc $ mv solver.f solver.F ~/dev/test/petsc $ make gfortran -c -fPIC -Wall -Wno-unused-variable -g -fopenmp -I/home/username/petsc/include -I/home/username

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 13:51, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: On Thu, 7 Nov 2013, at 12:56, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: ~/dev/test/petsc $ mv solver.f solver.F ~/dev/test/petsc $ make

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 13:59, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: What do you mean? (I don't think the program name has to be 'main'.). No, it doesn't. The name is meaningless in Fortran, but you need to use the keyword program. ~/dev/test/petsc

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 14:15, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: On Thu, 7 Nov 2013, at 13:59, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: What do you mean? (I don't think the program name has to be 'main

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 14:22, Svetlana Tkachenko wrote: On Thu, 7 Nov 2013, at 14:15, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: On Thu, 7 Nov 2013, at 13:59, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: What do

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 14:24, Svetlana Tkachenko wrote: On Thu, 7 Nov 2013, at 14:22, Svetlana Tkachenko wrote: On Thu, 7 Nov 2013, at 14:15, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: On Thu, 7 Nov 2013, at 13:59, Jed Brown wrote

Re: [petsc-dev] openmp

2013-11-06 Thread Svetlana Tkachenko
On Thu, 7 Nov 2013, at 12:56, Jed Brown wrote: Svetlana Tkachenko svetlana.tkache...@fastmail.fm writes: ~/dev/test/petsc $ mv solver.f solver.F ~/dev/test/petsc $ make gfortran -c -fPIC -Wall -Wno-unused-variable -g -fopenmp -I/home/username/petsc/include -I/home/username/petsc

[petsc-dev] openmp

2013-11-05 Thread Svetlana Tkachenko
I have configured petsc-dev (downloaded it today) with these options, and a small example. It appears to fail to compile without MPI with the error message: ./configure --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack --with-openmp --with-mpi=0 ~/dev/test/petsc $ echo