Re: [petsc-users] setting Sundials options within petsc4py

Veltz Romain writes: > OK will try, hopefully I have that python level. That would be adding changes > mainly to > petsc4py / src / PETSc / TS.pyx right? Yes. pgp7iLSoio0JK.pgp Description: PGP signature

Re: [petsc-users] setting Sundials options within petsc4py

On Apr 24, 2014, at 1:23 AM, Jed Brown wrote: > Veltz Romain writes: > >> Dear pets users, >> >> I have been trying to set some sundials options from petsc4py but I am >> unsuccessful. I hope somebody will have a good pointer. For example I >> can pass -ts_sundials_type bdf to python in the

Re: [petsc-users] KSP breakdown in specific cluster (update)

TAY wee-beng writes: > On 24/4/2014 9:41 AM, Barry Smith wrote: >>The numbers like >> >> sum = 1.9762625833649862e-323 >> >> rho = 1.9762625833649862e-323 >> >> beta = 1.600807474747106e-316 >> omega = 1.6910452843641213e-315 >> d2 = 1.571803252194

Re: [petsc-users] inserting vector into row of dense matrix

Randy, The problem is that MatMPIDense stores each row on a single process; while the solution in DMGetGlobalVector() will be spread over all processes. Thus you could use MatSetValues() as you have done to put the values in but it will be terribly inefficient. What I suggest is to i

[petsc-users] inserting vector into row of dense matrix

I have a 3D rectangular grid over which I am doing some computations. I have created a 3D DA for that grid, and I can put the results of the computations into a global vector obtained via DMGetGlobalVector. Now, I need to do these computations for several realizations, and I want to keep all th

Re: [petsc-users] KSP breakdown in specific cluster (update)

On 24/4/2014 9:41 AM, Barry Smith wrote: The numbers like sum = 1.9762625833649862e-323 rho = 1.9762625833649862e-323 beta = 1.600807474747106e-316 omega = 1.6910452843641213e-315 d2 = 1.5718032521948665e-316 are nonsense. They would genera

Re: [petsc-users] KSP breakdown in specific cluster (update)

The numbers like sum = 1.9762625833649862e-323 rho = 1.9762625833649862e-323 beta = 1.600807474747106e-316 omega = 1.6910452843641213e-315 d2 = 1.5718032521948665e-316 are nonsense. They would generally indicate that something is wrong, but unfortu

Re: [petsc-users] Problem with DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

Hi, Ya that's better. I've been trying to remove the unnecessary codes (also some copyright issues with the university) and to present only the code which shows the error. But strangely, by restructuring and removing all the redundant stuffs, the error is no longer there. I'm now trying to ad

Re: [petsc-users] KSP breakdown in specific cluster (update)

On 23/4/2014 6:00 PM, Matthew Knepley wrote: On Wed, Apr 23, 2014 at 5:55 AM, TAY wee-beng > wrote: Hi, Just to update that I managed to compare the values by reducing the problem size to hundred plus values. The matrix and vector are almost the same com

[petsc-users] Your email account

You have 1 new message Sign In Blackboard | Technology Services

Re: [petsc-users] setting Sundials options within petsc4py

Veltz Romain writes: > Dear pets users, > > I have been trying to set some sundials options from petsc4py but I am > unsuccessful. I hope somebody will have a good pointer. For example I > can pass -ts_sundials_type bdf to python in the command line. > > I would like to call TSSundialsSetToleran

Re: [petsc-users] MATLAB and PetSc

On Apr 23, 2014, at 2:06 PM, Umut Tabak wrote: > On 04/23/2014 08:59 PM, Barry Smith wrote: >> >>No mex is not used; instead the MATLAB function calllib() is used which >> is a million times easier than messing with mex. To use PETSc from Matlab >> you need to ./configure PETSc with the o

Re: [petsc-users] MATLAB and PetSc

On 04/23/2014 08:59 PM, Barry Smith wrote: No mex is not used; instead the MATLAB function calllib() is used which is a million times easier than messing with mex. To use PETSc from Matlab you need to ./configure PETSc with the option —with-matlabIt does not work currently in parallel

Re: [petsc-users] MATLAB and PetSc

On Apr 23, 2014, at 7:36 AM, Sudheer Kumar wrote: > Hi Barry, > > I have looked at code in bin/matlab/classes. It is not clear to me how to use > the code. Examples are in the examples/tutorials subdirectory of that directory. > Great if you could help me how to use these codes. So, with

[petsc-users] setting Sundials options within petsc4py

Dear pets users, I have been trying to set some sundials options from petsc4py but I am unsuccessful. I hope somebody will have a good pointer. For example I can pass -ts_sundials_type bdf to python in the command line. I would like to call TSSundialsSetTolerance to set the tolerance. Hence I ne

Re: [petsc-users] install Slepc wth Arpack

El 23/04/2014, a las 17:02, Lorenzo Zanon escribió: > Hello, > > I would like to install Slepc with Arpack, but the procedure I'm using is not > working. > > I downloaded arpack-ng and configured like this: > > ./configure --prefix=$HOME/Desktop/arpack_dir > > So the libraries can be found i

[petsc-users] install Slepc wth Arpack

Hello, I would like to install Slepc with Arpack, but the procedure I'm using is not working. I downloaded arpack-ng and configured like this: ./configure --prefix=$HOME/Desktop/arpack_dir So the libraries can be found inside arpack_dir. Then I set the PETSC_DIR, PETSC_ARCH and SLEPC_DIR. F

Re: [petsc-users] MATLAB and PetSc

Hi Barry, I have looked at code in bin/matlab/classes. It is not clear to me how to use the code. Great if you could help me how to use these codes. So, with these codes can we run the parallel version of Petsc code from MATLAB. Is mex interface involved in compiling the Petsc code to MATLAB execu

Re: [petsc-users] KSP breakdown in specific cluster (update)

On Wed, Apr 23, 2014 at 5:55 AM, TAY wee-beng wrote: > Hi, > > Just to update that I managed to compare the values by reducing the > problem size to hundred plus values. The matrix and vector are almost the > same compared to my win7 output. > Run in the debugger and get a stack trace, Matt

Re: [petsc-users] KSP breakdown in specific cluster (update)

Hi, Just to update that I managed to compare the values by reducing the problem size to hundred plus values. The matrix and vector are almost the same compared to my win7 output. Also tried valgrind but it aborts almost immediately: valgrind --leak-check=yes ./a.out ==17603== Memcheck, a mem

[petsc-users] KSP breakdown in specific cluster

Hi, My code was found to be giving error answer in one of the cluster, even on single processor. No error msg was given. It used to be working fine. I run the debug version and it gives the error msg: [0]PETSC ERROR: [0