make test errors

2009-08-04 Thread Satish Balay
Try running manually - and see if the example works. These messages could be verbose - and misleading. i.e cd src/ksp/ksp/examples/tutorials make ex2 mpiexec -n 2 ./ex2 Or you can avoid this by installing mpich - instead of openmpi. Satish On Tue, 4 Aug 2009, Yixun Liu wrote: I'm not

make test errors

2009-08-05 Thread Satish Balay
file. Contact your local System Administrator to confirm the availability of the interfaces in the dat.conf file. Satish Balay wrote: Try running manually - and see if the example works. These messages could be verbose - and misleading. i.e cd src/ksp/ksp/examples/tutorials make

PETSC_LIB

2009-08-09 Thread Satish Balay
If the examples don't work [with the provided makefile] - you'll have to check with the person who built this petsc for you.. Perhaps he needs to reinstall petsc - or just install blas/lapack for this to work.. Satish On Sun, 9 Aug 2009, irfan.khan at gatech.edu wrote: Thank for your reply

How to get information about the off-diagonal submatrix

2009-08-10 Thread Satish Balay
alternatively - just access the private datastructe and print.. #include ../src/mat/impls/aij/mpi/mpiaij.h Mat_MPIAIJ *mpiaij = (Mat_MPIAIJ*)mat-data; MatView(mpiaij-B,...); Satish On Mon, 10 Aug 2009, Matthew Knepley wrote: If this is just for debugging, the easiest way I think is to

cannot find -lustdio with matlab on

2009-08-10 Thread Satish Balay
You don't need to install petsc with matlab for this. You can just use '-mat_view_matlab' - or PETSC_VIEWER_MATLAB in your code. Satish On Mon, 10 Aug 2009, Weiwei Wang wrote: Dear All, I want to check the determinant of some sparse matrix in petsc but I didn't know how to do it using

How to get information about the off-diagonal submatrix

2009-08-10 Thread Satish Balay
for debugging!) I've expected to have an empty off-diagonal matrix. But the matrix is not empty. But how is it possible, when all entries, which are inserted to the matrix via MatSetValues, have indices for row and col within the ownership range of the processors matrix? Thomas Satish

porting petsc to a nec sx9 computer running under unix super ux

2009-08-11 Thread Satish Balay
we had an old version - perhaps 2.2.x working on sx6. It has standard compilers with mpi and blas - so I'm guessing current petsc should work on sx9. If you encounter problems building on sx9 - then send us the error logs at petsc-maint at mcs.anl.gov Satish On Tue, 11 Aug 2009, Joel Guerrero

porting petsc to a nec sx9 computer running under unix super ux

2009-08-11 Thread Satish Balay
how to cross compile petsc on these platforms?. I am using sx9 compilers (sxcc, sxc++, sxmpicc and so on). Thanks for your help, Joel - Original Message From: Satish Balay balay at mcs.anl.gov To: PETSc users list petsc-users at mcs.anl.gov Sent: Wednesday, August 12

porting petsc to a nec sx9 computer running under unix super ux

2009-08-12 Thread Satish Balay
lapack/blas libraries. Hope you can help me, Joel - Original Message From: Satish Balay balay at mcs.anl.gov To: PETSc users list petsc-users at mcs.anl.gov Sent: Wednesday, August 12, 2009 4:54:40 AM Subject: Re: porting petsc to a nec sx9 computer running

cast PetscScalar to double

2009-08-19 Thread Satish Balay
you can replace (double)(var) with one of the following [approprate to your code] thingy. PetscAbsScalar(var) PetscRealPart(var) Satish On Wed, 19 Aug 2009, Fredrik Bengzon wrote: Hi I'm using a Petsc 3 installation compiled with support for complex numbers. However, in the code I have made

Petsc with Intel MKL for blas/lapack

2009-08-25 Thread Satish Balay
Actually we now use a notation with MKL library names that appears to work consistantly across MKL versions. So there could be other issues. configure.log will clear things up [so send it over to us at petsc-maint at mcs.anl.gov] Satish On Tue, 25 Aug 2009, Matthew Knepley wrote: The

Create With Array

2009-08-31 Thread Satish Balay
No new memory is allocated. It uses the array specified by the user.. Satish On Mon, 31 Aug 2009, tsjb00 wrote: Hi! I have a quick question about the 'CreateWith Array' subroutines. Do they increase memory requirement? I am wondering if new memory space is allocated for the vector/matrix

HDF5 and DA ndof 1

2009-09-02 Thread Satish Balay
Easiest thing to do is isolate the changeset - and grab the diff - and apply to petsc-3.0.0. If you know the changeset id [or the app lines that belong to this changeset] - I can take a look. If its multiple changesets for this fix - then we'll have to find/extract them all.. [finding all might

memory reporting question

2009-09-18 Thread Satish Balay
There is no destroy call between PetscMemoryGetCurrentUsage(), PetscMemoryGetMaximumUsage - so MaxUsage isn't updated. Perhaps it should be updated in PetscMemoryGetCurrentUsage() call aswell? Satish On Fri, 18 Sep 2009, Matt Funk wrote: Ok, i am more confused now. The following is part

where to realize petscsetcommonblock_ and petscgetcommoncomm_

2009-09-22 Thread Satish Balay
petscgetcommoncomm_ should be in somefort.F Can you verify if you've configured/compiled PETSc with a fortran compiler. And also verify that there are no errors during petsc build [i.e in make.log] asterix:/home/balay/tmp/petsc-dist-test/asterix64/libnm -Ao *.a |grep -i petscgetcommoncomm

Where is sles in PETSc 3.0.0?

2009-09-24 Thread Satish Balay
SLES has been renamed/merged-into KSP a few releases back.. Satish On Thu, 24 Sep 2009, I.Gutheil wrote: Hello all, we installed the new PETSc 3.0.0-p6 version of PETSc on a new computer and now a user is missing petscsles.h. When I looked into the src directory I saw that there is no

Building Issues on BGL

2009-10-13 Thread Satish Balay
You have to remove the defines or rename them to something else in petscconf.h. Just commenting or undefining doesn't work. [this is because our build system does a grep on this file - so its doesn't exactly follow preprocessing directives] Regarding the original issue - I'm not sure why it

PaStiX crash

2009-10-20 Thread Satish Balay
Perhaps you can try running in valgrind to see here the problem is. You can also try --with-debugging=0 COPTFLAGS='-g -O' - and see if it crashes. If so - run in a debugger to determine the problem. Satish On Tue, 20 Oct 2009, Andreas Grassl wrote: Hello, I wanted to use PaStix and have

How to get more detailed output in make.log

2009-10-20 Thread Satish Balay
Usually the header at the begining lists these compiler options separately. But if you want the extra info - you can do: make ACTION=lib tree Note: this is not exactly the same build as the default build ['libfast' attempts to compile multiple files with a sigle command - whereas 'lib' does one

[petsc-maint #36400] Re: C++ error! MPI_Finalize() could not be located!

2009-10-21 Thread Satish Balay
Since you are using hypre - also use the additional option -with-cxx=icpc [by default its picking up g++ - which is causing this incompatibility with icc etc..] Satish On Wed, 21 Oct 2009, Dominik Szczerba wrote: configure.log is attached. Looking forward your feedback. Dominik

One simple question about complex number

2009-10-22 Thread Satish Balay
we use 'std:complex' with --with-clanguage=cxx [and with c - we use c99 complex support]. I'm not sure what stl::complex is. Satish On Thu, 22 Oct 2009, Yujie wrote: Hi, PETSc Developers, What is the difference between PETSc complex number and stl::complex. If I define a variable var using

set a column of a matrix

2009-11-01 Thread Satish Balay
On Sun, 1 Nov 2009, Umut Tabak wrote: On Sun, Nov 01, 2009 at 07:44:58AM -0600, Barry Smith wrote: MatSetValues() so long as you have good matrix preallocation this will work fine. Doing a transpose is very expensive. Dear Barry and Matt, Thanks for the replies, I am

A problem on compiling a petsc code.

2009-11-11 Thread Satish Balay
Please send the complete compile command that gave this error. Also do PETSc examples compile/run correctly? What do you get for: cd src/ksp/ksp/examples/tutorials make ex2 make ex2f Satish On Tue, 10 Nov 2009, Yin Feng wrote: When I compiled petsc code, I encourntered following error code.

A problem on compiling a petsc code.

2009-11-11 Thread Satish Balay
} and I compiled ex2 successfully. Thank you! On Wed, Nov 11, 2009 at 12:49 AM, Satish Balay balay at mcs.anl.gov wrote: Please send the complete compile command that gave this error. Also do PETSc examples compile/run correctly? What do you get for: cd src/ksp/ksp/examples

scaling in 4-core machine

2009-11-18 Thread Satish Balay
Just want to add one more point to this. Most multicore machines do not provide scalable hardware. [yeah - the FPUs cores are scalable - but the memory subsystem is not]. So one should not expect scalable performance out of them. You should take the 'max' performance you can get out out them -

suprelu

2009-11-20 Thread Satish Balay
works fine for me. Perhaps you can retry after 'rm -rf externalpackages' If the problem persists - send configure.log to petsc-maint at mcs.anl.gov Satish On Fri, 20 Nov 2009, Lukasz Kaczmarczyk wrote: Hello I try to configure petsc (p9) with SuperLU and I get error message

how to run petsc on different machine

2009-12-01 Thread Satish Balay
--download-mpich=1 defaults to gforker i.e a single machine. If you wish to run on multiple machines - then you'll have to figureout what mechanism you'd like to use - and then direct the MPI install for that method. For eg: - the common method is to use MPD. So you can install petsc with

optimized complex number kernel

2009-12-02 Thread Satish Balay
The performance difference between c/c++ kernels and fortran kernels - or c-complex/c++-complex can be compiler specific. You can build all relavent versions [each with a different PETSC_ARCH] and see if it makes a difference for you code. The relavent configure options are:

How to visualize a matrix object from a remote machine?

2009-12-15 Thread Satish Balay
On Tue, 15 Dec 2009, hxie at umn.edu wrote: Hi, I use X11 to connect the server and I want to draw the nonzero structure of matrix. The option '-mat_view_draw' gives some error messages. How to fix it? Thanks. What error message? Normally you can do the following: ssh -Y remotemachine

How to visualize a matrix object from a remote machine?

2009-12-15 Thread Satish Balay
) From: Satish Balay balay at mcs.anl.gov Subject: Re: How to visualize a matrix object from a remote machine? To: PETSc users list petsc-users at mcs.anl.gov Message-ID: alpine.LFD.2.00.0912151056190.20137 at asterix Content-Type: TEXT/PLAIN; charset=US-ASCII On Tue, 15 Dec 2009, hxie

Problems with VecGetArray in Fortran

2009-12-15 Thread Satish Balay
In fortran - try and avoid passing constants with implicit types - and explicitly define them. i.e PetscScalar one one = 1.0d0 call VecSet(tmp, one, ierr) etc.. Also - VecGetArrayF90() can be prefered over the old VecGetArray() [assuming f77 is nolonger needed] Satish

installation problems in mac machine

2009-12-21 Thread Satish Balay
On Mon, 21 Dec 2009, Matthew Knepley wrote: Installation is almost always done with root permission. If you want to be pendatic, you can create other users for specific tasks, like the webserver, but you STILL sudo to that user. In most cases - the part thats done as 'root' is minimised [to

[petsc-users] tag in subjects

2009-12-21 Thread Satish Balay
added now. Satish On Mon, 21 Dec 2009, Dominik Szczerba wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Would it be desirable to have [Petsc] automatically prepended to all list messages? I am receiving loads of emails daily from multiple mailing lists and it would greatly help

[petsc-users] matrix visualization problem

2009-12-23 Thread Satish Balay
well if you need -mat_view_draw - then suggest editing src/mat/impls/aij/mpi/mpiaij.c and increate this cut-off size. [there is now a command-line option to over-ride this limit in petsc-dev] For binary viewer - check mat/examples/tests/ex31.c Satish On Wed, 23 Dec 2009, hxie at umn.edu wrote:

[petsc-users] matrix visualization problem

2009-12-23 Thread Satish Balay
On Wed, 23 Dec 2009, hxie at umn.edu wrote: Hi, Thanks. When I try to run with the option -ksp_xmonitor, no window pops up. Do I need to reinstall the petsc with the option --with-x=1? use -ksp_monitor_draw BTW, if I want to use some external packages, do I have to reinstall the petsc?

[petsc-users] matrix visualization problem

2009-12-23 Thread Satish Balay
Binary viewer is to dump the matrix into a file - not a mat_view_draw wrt -mat_view_draw - I've already indicated what to do.. satish On Wed, 23 Dec 2009, hxie at umn.edu wrote: Hi, I put the following in my code: -- PetscViewer viewer ... call

[petsc-users] petsc 3.0 vs petsc 2.3

2010-01-13 Thread Satish Balay
Is this application making PETSc calls from fortran? You'll have to go through the changes file - and make sure the prototype changes are applied to your code [a fortran compiler doesn't help with prototype checks here]. http://www.mcs.anl.gov/petsc/petsc-as/documentation/changes/300.html [If the

[petsc-users] On what condition is useful MPI-based solution?

2010-01-13 Thread Satish Balay
On Wed, 13 Jan 2010, Jose E. Roman wrote: On 13/01/2010, Takuya Sekikawa wrote: I have a few more questions. Just in case I want to ask. [Q1] In general, like ex1, Don't we need any change to source code to run program with MPI-based system? in other words, to run

[petsc-users] x64 build problem

2010-01-17 Thread Satish Balay
Actually its --with-cc='win32fe cl --nodetect' --with-fc='win32fe ifort --nodetect' etc.. check 'Using Visual Studio 2005:' section of installation instructions http://www.mcs.anl.gov/petsc/petsc-as/documentation/installation.html Satish On Sun, 17 Jan 2010, Barry Smith wrote: Try

[petsc-users] using PETSC with visual studio 2005?

2010-01-21 Thread Satish Balay
Check the installation instructions - specifically - the section: 'Using Visual Studio 2005' http://www.mcs.anl.gov/petsc/petsc-as/documentation/installation.html Satish On Thu, 21 Jan 2010, Zhifeng Sheng - EWI wrote: Dear all Does anyone use Petsc in visual studio 2005? I have used petsc

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
Do you have to use proxy for your http connections? If so - perhaps you can try adding the settings to ~/.hgrc - as this url indicates and see if it fixes the problem. Satish On Mon, 25 Jan 2010, (Rebecca) Xuefei YUAN wrote: Hi,All, I try to install a pets-dev locally according to the

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
Sorry - forgot to mention the url http://mail.openjdk.java.net/pipermail/build-dev/2007-December/000628.html Satish On Mon, 25 Jan 2010, Satish Balay wrote: Do you have to use proxy for your http connections? If so - perhaps you can try adding the settings to ~/.hgrc - as this url indicates

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
! Rebecca Quoting Satish Balay balay at mcs.anl.gov: Sorry - forgot to mention the url http://mail.openjdk.java.net/pipermail/build-dev/2007-December/000628.html Satish On Mon, 25 Jan 2010, Satish Balay wrote: Do you have to use proxy for your http connections

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
On Mon, 25 Jan 2010, (Rebecca) Xuefei YUAN wrote: Dear Satish, I got these from those two commands: rebecca at YuanWork:~/soft/petsc-dev$ env |grep -i proxy http_proxy=http://:8080/ no_proxy=localhost,127.0.0.0/8,*.local rebecca at YuanWork:~/soft/petsc-dev$ env |grep -i http

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
for this build. Satish On Mon, 25 Jan 2010, Satish Balay wrote: On Mon, 25 Jan 2010, (Rebecca) Xuefei YUAN wrote: Dear Satish, I switched to install it manually What do you mean? Did you get the nightly tarball? and encounter an error when make test. Did the build go through

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
Gald it works now. BTW: What OS/version do you have? And what version of mercurial is this? [just for completeness] Satish On Mon, 25 Jan 2010, (Rebecca) Xuefei YUAN wrote: Dear Satish, There is no such an error showed up this time. Thanks very much! Rebecca Quoting Satish

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
/c2html.tar.gz The Fortran compiler I have is g77 and gfortran. If I do not use Fortran, will this be fine? Cheers, Rebecca Quoting Satish Balay balay at mcs.anl.gov: Ok - is your manual install '2. access tarballs directly from our mercurial server.? http://www.mcs.anl.gov

[petsc-users] Errors while install pets-dev

2010-01-25 Thread Satish Balay
, (Rebecca) Xuefei YUAN wrote: I do not know why the connection time is out, but for completeness, I am using Ubuntu 8.04, and the mercurial version is Mercurial Distributed SCM (version 0.9.5). Cheers, Rebecca Quoting Satish Balay balay at mcs.anl.gov: Gald it works now

[petsc-users] Compile Petsc in visual studio 2003 with /MD option

2010-02-02 Thread Satish Balay
Please do not cross post to multiple lists. If you have installation issues - send the relavent configure.log to petsc-maint at mcs.anl.gov You should be able to specify 'CFLAGS=-MT' option to configure and both PETSc and f2cblaslapack should use it. [assuimg you are using --with-mpi=0] Satish

[petsc-users] Mesh class (fwd)

2010-02-02 Thread Satish Balay
-- Forwarded message -- Date: Tue, 2 Feb 2010 15:15:20 + From: Denner, Fabian f.denne...@imperial.ac.uk To: petsc-users-bounces at mcs.anl.gov petsc-users-bounces at mcs.anl.gov Subject: Mesh class Hi everybody, are there any examples available, demonstrating the use of

[petsc-users] Compile Petsc in visual studio 2003 with /MD option

2010-02-02 Thread Satish Balay
- send configure.log, make.log test.log from PETSC_DIR/PETSC_ARCH/conf to petsc-maint] Satish On Tue, 2 Feb 2010, Satish Balay wrote: Please do not cross post to multiple lists. If you have installation issues - send the relavent configure.log to petsc-maint at mcs.anl.gov You should

[petsc-users] A Question about Vec

2010-02-03 Thread Satish Balay
On Wed, 3 Feb 2010, Ryan Yan wrote: Hi all, Can someone help me to understand the Vec? My question is: in the following piece of code, will a be changed, if b is changed by the subroutine SNESComputeFunction? I did see that both a and b are both pointers to an object. Yes - all petsc

[petsc-users] query on meaning of libmpiuni.*

2010-03-01 Thread Satish Balay
On Mon, 1 Mar 2010, Toby D. Young wrote: One final question please: I understand that this library can be built as static or dynamic; right? Then I look at my build, and only libmpiuni.a appears. Why is that? Can a libmpiuni.so exist? I build --with-shared=1. Yes - is is the

[petsc-users] make error from the ex1.c

2010-03-20 Thread Satish Balay
Can't flood petsc-users mailing list with huge configure.log. petsc-maint is more apporpriate for such issues. asterix:/home/balay/download-pinecat makefile.test ALL: ex1 CFLAGS = ${PETSC_CC_INCLUDES} FFLAGS = ${PETSC_FC_INCLUDES} include ${PETSC_DIR}/conf/variables ex1: ex1.o

[petsc-users] Linking error with C++ code: undefined symbols

2010-03-23 Thread Satish Balay
For one we recommend using the latest version i.e petsc-3.0.0 [or the upcoming petsc-3.1 - via petsc-dev] - but not this old version 2.3.3. Secondly - if you want to use PETSc from c++ - the recommended way is to build it with the configure option '--with-clanguage=cxx' And then [after making

[petsc-users] Configuring PETSc with Intel compiler.

2010-03-30 Thread Satish Balay
You should use CXX=icpc Satish On Tue, 30 Mar 2010, Pedro Torres wrote: Hello, I'm trying to compile PETSc with intel compiler and get the error below. I can compile a c/c++ source code without problem. What am I doing wrong?. Thanks a lot ./config/configure.py

[petsc-users] configure failed on MPI and Intel Compiler environment

2010-04-07 Thread Satish Balay
petsc-3.0 [with latest patches] configure should be able to handle MLK-10. Don't remember MLK-11. without MPI (--with-mpi=0), this configure succeeded with no problem. Claims MLK-10 was detected properly without mpi. This is weird - but again configure.log for both builds will clarify the

[petsc-users] Improve of win32fe

2010-04-13 Thread Satish Balay
We have now added a copyright notice in the source tree. Its similar to petsc license - and is cygwin.dll compatible [but not GPL] http://petsc.cs.iit.edu/petsc/win32fe/raw-file/71989470a6d3/readme.html Satish On Sat, 27 Mar 2010, Gong Ding wrote: Dear Petsc developers For long time I am

[petsc-users] Error during compiling my own code

2010-04-16 Thread Satish Balay
/app1/mvapich2/current/include/mpif.h(9): error #6401: The attributes of this name conflict with those made accessible by a USE statement. [MPI_SOURCE] INTEGER MPI_SOURCE, MPI_TAG, MPI_ERROR You are using some module - based on the 'use' statement refered by the error above. And

[petsc-users] array and vec

2010-04-16 Thread Satish Balay
On Fri, 16 Apr 2010, Matthew Knepley wrote: 2010/4/15 Leping Chen? chenleping at yahoo.cn hello, petsc teams, when I use FormFuncion(), I need use a value of x; but x is defined by Vec, I don't know how to get a value of x, for example, the 3rd value. Exactly as Jed

[petsc-users] about Vecgetarray()

2010-04-17 Thread Satish Balay
On Fri, 16 Apr 2010, Leping Chen? wrote: PetscOffset xx_i Vec x double precision u(6) call VecGetArray(x,u,xx_i,ierr) do 30 i=1,n u(i) = 1000.0*i u(i+xx_i) = 1000.0*i 30 continue call

[petsc-users] configure script for petsc

2010-04-23 Thread Satish Balay
On Fri, 23 Apr 2010, Stephan Kramer wrote: Hi all, I'm trying to write a portable configure script for our software that uses petsc, and needs to deal with petsc installations on a number of different platforms. In petsc 3.0 I could do make -f $(PETSC_DIR)/conf/base getlinklibs and

[petsc-users] configure script for petsc

2010-04-23 Thread Satish Balay
On Fri, 23 Apr 2010, Stephan Kramer wrote: Satish Balay wrote: On Fri, 23 Apr 2010, Stephan Kramer wrote: Hi all, I'm trying to write a portable configure script for our software that uses petsc, and needs to deal with petsc installations on a number of different platforms

[petsc-users] [Supercomputing Lab #847] VecView() error in BlueGene

2010-04-28 Thread Satish Balay
! Rebecca Quoting Satish Balay balay at mcs.anl.gov: It should be used at runtime. For ex on our bgl - we use qsub as: qsub --env BG_MAXALIGNEXP=-1 -t 00:15:00 -n 4 ./ex2 Satish On Tue, 27 Apr 2010, (Rebecca) Xuefei YUAN wrote: Dear Satish, I think

[petsc-users] any ways to run PETSc4py on several CPU?

2010-05-10 Thread Satish Balay
Note: this might not work with some older MPI impls like mpich1. However newer ones like mpich2 should work. mpiexec -n 4 python foo.py Satish On Mon, 10 May 2010, Matthew Knepley wrote: You use mpirun. Matt On Sun, May 9, 2010 at 10:42 AM, Dmitrey dmitrey15 at ukr.net wrote: hi

[petsc-users] Error with file linking

2010-05-11 Thread Satish Balay
, Satish Balay wrote: Why not use petsc makefiles - and avoid these hasseles?

[petsc-users] Error with file linking

2010-05-12 Thread Satish Balay
On Wed, 12 May 2010, Satish Balay wrote: Attaching the modified makefile. It might need a couple of iterations of fixing. Also should remove '-c' from FFLAGS Satish

[petsc-users] Error with file linking

2010-05-12 Thread Satish Balay
-Beng Tay On 13/5/2010 1:27 AM, Satish Balay wrote: Attaching the modified makefile. It might need a couple of iterations of fixing. Hopefully its clear why it would work. Note - if not for different flags required for different sources (opt/opt2/opt3) - you can build all petsc/non

[petsc-users] complex number issues

2010-05-19 Thread Satish Balay
I don't understand the code change - but if you want to use c99 complex - you can build petsc with: --with-scalar-type=complex --with-clanguage=c [which is the default for clanguage] And no code-change needed for PETSc. However if your code uses c++ - you might have to address that.. Satish

[petsc-users] small dense solves

2010-06-21 Thread Satish Balay
PETSc just uses native lapack. There is nothing special here. petscblaslapack have c-prototypes for the native blas/lapack - so there is some error checking in petsc usage of blas/lapack. Satish On Mon, 21 Jun 2010, Andrew T Barker wrote: As one part of a larger parallel code, I want to do

[petsc-users] [petsc-dev] updating values in a DA Global array

2010-06-24 Thread Satish Balay
Actually gfortran-4.3+ should work [gfortran 4.2 is broken] Satish On Thu, 24 Jun 2010, Barry Smith wrote: Yes. On Jun 24, 2010, at 3:49 PM, Jed Brown wrote: On Thu, 24 Jun 2010 15:46:05 -0500, Barry Smith bsmith at mcs.anl.gov wrote: It turns out there was a bug in

[petsc-users] ublas sparse matrix bindings?

2010-06-29 Thread Satish Balay
On Tue, 29 Jun 2010, Luke Bloy wrote: Hi, In one of my projects I use the ublas compressed row major format for storing sparse matricies. One of the reasons i chose this format was interoperability with other packages. I'd like to use some petsc and slepc functionality in my project but I

[petsc-users] ublas sparse matrix bindings?

2010-06-29 Thread Satish Balay
On Tue, 29 Jun 2010, Luke Bloy wrote: Thanks for the quick responses. my matrices are fairly large 80,000 x 80,000 and not too sparse so I'd like to avoid any serialization if possible. I've worked out how to use MatCreateSeqAIJWithArrays so long as i control the types used in my ublas

[petsc-users] MatView for huge matrix output

2010-06-29 Thread Satish Balay
Perhaps you can build PETSc with the additional flag '--with-large-file-io=1' and see if this helps.. Satish On Tue, 29 Jun 2010, Yujie wrote: Dear PETSc developers, I want to output an about 36K*36K dense matrix using MatView in binary format. I use RedHat Enterprise 5 64bits system.

[petsc-users] Petsc parallel processing from a library, without mpiexec -np ...

2010-07-06 Thread Satish Balay
One additional point: If using 'parallel processing' - then you shouldn't be trying to circumvent this part of MPI standard [mpiexec -n]. And if you are primarily doing sequential stuff - you can build PETSc with '--with-mpi=0' [so MPI is not used. However when working with multiple MPI software

[petsc-users] my code runs too slow

2010-07-07 Thread Satish Balay
used during MatSetValues calls =1830 not using I-node routines On Jul 7, 2010, at 12:51 PM, Satish Balay wrote: total: nonzeros=1830, allocated nonzeros=29280 total number of mallocs used during MatSetValues calls =1830 There is something wrong with your

[petsc-users] my code runs too slow

2010-07-07 Thread Satish Balay
,TSDefaultComputeJacobianColor,matfdcoloring); These are the Jacobian related codes. On Jul 7, 2010, at 1:51 PM, Satish Balay wrote: total: nonzeros=1830 mallocs used during MatSetValues calls =1830 Looks like you are zero-ing out the non-zero structure - before assembling the matrix. Are you

[petsc-users] my code runs too slow

2010-07-07 Thread Satish Balay
at psu.edu wrote: On Jul 7, 2010, at 2:06 PM, Satish Balay wrote: On Wed, 7 Jul 2010, Xuan YU wrote: ierr = MatCreateSeqAIJ(PETSC_COMM_SELF,N,N,10,PETSC_NULL,J);CHKERRQ(ierr); ierr = MatAssemblyBegin(J,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr); ierr = MatAssemblyEnd(J,MAT_FINAL_ASSEMBLY

[petsc-users] my code runs too slow

2010-07-07 Thread Satish Balay
: [0]PETSC ERROR: MatGetColoring() line 491 in src/mat/color/color.c [0]PETSC ERROR: VecRestoreArrayPrivate3() line 201 in /work/petsc/pihm.c On Jul 7, 2010, at 2:28 PM, Satish Balay wrote: ierr = SNESComputeJacobian(ts_snes,CV_Y,J,J,flag);CHKERRQ(ierr); Perhaps

[petsc-users] How to get Version

2010-07-15 Thread Satish Balay
For eg: #if (PETSC_VERSION_RELEASE == 1) write(6,25) PETSC_VERSION_MAJOR,PETSC_VERSION_MINOR,PETSC_VERSION_SUBMINOR,PETSC_VERSION_PATCH 25 format('PETSc Release Version: ',i1,'.',i1,'.',i1,' Patch:',i2) #else write(6,25) 25 format('PETSc Development Version: ',

[petsc-users] configure problems

2010-07-21 Thread Satish Balay
Try: ./config/configure.py PETSC_DIR=$PWD PETSC_ARCH=win64-msvc-release --download-c-blas-lapack=1 --with-x=0 --with-debugging=0 --with-cc='win32fe cl' --with-fc=0 If it fails - send us configure.log at petsc-maint at mcs.anl.gov BTW: What do you have in your LIB env variable? [MS search this

[petsc-users] configure problems

2010-07-21 Thread Satish Balay
to get a usable MPICH2 installation, Dominik On Wed, Jul 21, 2010 at 6:09 PM, Satish Balay balay at mcs.anl.gov wrote: Try: ./config/configure.py PETSC_DIR=$PWD PETSC_ARCH=win64-msvc-release --download-c-blas-lapack=1 --with-x=0 --with-debugging=0 --with-cc='win32fe cl' --with-fc=0

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
sh: Warning: win32fe: File Not Found: /cygdrive/c/Program Files/MPICH2/lib/amd64/msmpi.lib Looks like - if win32fe doesn't find the file - it doesn't try to convert the path. This is not the cause of his configure failure. The primary issue is: [which is much ahead in configure.log] sh:

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
, Jul 22, 2010 at 3:45 PM, Satish Balay balay at mcs.anl.gov wrote: sh: Warning: win32fe: File Not Found: /cygdrive/c/Program Files/MPICH2/lib/amd64/msmpi.lib Looks like - if win32fe doesn't find the file - it doesn't try to convert the path. This is not the cause of his configure failure

[petsc-users] Missing symbols in libpetscts.h

2010-07-22 Thread Satish Balay
TSSetIFunction() doesn't exist in petsc-3.0.0. The same code will not work with both versions of petsc. You should stick with one version [i.e the latest ..] Satish On Thu, 22 Jul 2010, Filippo Spiga wrote: Dear all, I'm trying to compile PETSc 3.0.0-p12 and PETSc 3.1-p3. The

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
Completed test examples On Thu, Jul 22, 2010 at 4:13 PM, Satish Balay balay at mcs.anl.gov wrote: On Thu, 22 Jul 2010, Dominik Szczerba wrote: Many thanks for your precious feedback. My MPICH install is 64 bit but it seems that for some reason 32 bit VC compiler is invoked (I have both

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
Should have said: cd src/ksp/ksp/examples/tutorials make ex2 mpiexec -n 2 ./ex2.exe mpiexec -localonly -n 2 ./ex2.exe Satish On Thu, 22 Jul 2010, Satish Balay wrote: Try running manually. MPICH2 install by default insists on passwd - and this process isn't handled by petsc scripts properly

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
or directory make: *** No rule to make target `/conf/test'. Stop. On Thu, Jul 22, 2010 at 10:11 PM, Satish Balay balay at mcs.anl.gov wrote: Should have said: cd src/ksp/ksp/examples/tutorials make ex2 mpiexec -n 2 ./ex2.exe mpiexec -localonly -n 2 ./ex2.exe Satish On Thu

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
On Thu, Jul 22, 2010 at 10:23 PM, Satish Balay balay at mcs.anl.gov wrote: Sorry - you need PETSC_ARCH and PETSC_DIR with make. [normally these are set as env variables] make PETSC_DIR=/cygdrive/c/Users/Dominik/Programs/petsc-3.1-p3 PETSC_ARCH=win64-msvc-release ex2 Satish

[petsc-users] configure problems

2010-07-22 Thread Satish Balay
) failed, exhaused all end points Unable to connect to 'Tharsis:8676', sock error: Error = -1 ReadFile() failed, error 109 unable to start the local smpd manager. On Thu, Jul 22, 2010 at 10:52 PM, Satish Balay balay at mcs.anl.gov wrote: Have you tried 'mpiexec -localonly

[petsc-users] configure problems

2010-07-23 Thread Satish Balay
. On Thu, Jul 22, 2010 at 10:52 PM, Satish Balay balay at mcs.anl.gov wrote: Have you tried 'mpiexec -localonly'? Satish On Thu, 22 Jul 2010, Dominik Szczerba wrote: I tried: $ /cygdrive/c/Program\ Files/MPICH2/bin/mpiexec.exe -np 4 ./ex2

[petsc-users] Missing symbols in libpetscts.h

2010-07-23 Thread Satish Balay
On Fri, 23 Jul 2010, Jed Brown wrote: On Fri, 23 Jul 2010 08:47:39 -0400, Filippo Spiga filippo.spiga at disco.unimib.it wrote: Dear Satish, thanks for you reply. Anyway not in all case the libpetscts.{a, so} is compiled. Is there a specific option of the configure that build

[petsc-users] configure problems

2010-07-23 Thread Satish Balay
On Fri, 23 Jul 2010, Dominik Szczerba wrote: On Fri, Jul 23, 2010 at 4:06 PM, Satish Balay balay at mcs.anl.gov wrote: On Fri, 23 Jul 2010, Dominik Szczerba wrote: OK so to conclude this thread, on Windows I need to run the command line window in Administrator mode, and start

[petsc-users] configure problems

2010-07-23 Thread Satish Balay
On Fri, 23 Jul 2010, Dominik Szczerba wrote: Unfortunately, I need also a Windows version and I will leave the cygwin option as the last resort. For the moment I managed to successfully build petsc with --with-parmetis=1 --download-parmetis=1. When using --download-hypre=1 configure fails

[petsc-users] about speed

2010-07-26 Thread Satish Balay
Perhaps you missed a step here. After 'configure' you shoud do 'make all' to rebuild the PETSc libraries. Then you should try 'make test' to verify the build. After that - attempt to build your code. satish On Mon, 26 Jul 2010, Xuan YU wrote: On Jul 26, 2010, at 3:53 PM, Barry Smith wrote:

[petsc-users] run sundials from Petsc

2010-07-27 Thread Satish Balay
You list one compile with linux-gnu-c-debug - and the other with linux-gnu-c-opt. Looks like you rebuilt linux-gnu-c-debug with sundials - but not linux-gnu-c-opt build. Satish On Tue, 27 Jul 2010, Xuan YU wrote: On Jul 27, 2010, at 6:22 PM, Jed Brown wrote: On Tue, 27 Jul 2010 18:14:46

[petsc-users] run sundials from Petsc

2010-07-28 Thread Satish Balay
Perhaps your previous build is with --with-mpi=0 Just do: ./configure --download-sundials=1 --download-mpi=1 --download-f-blas-lapack=1 --with-debugging=0 PETSC_ARCH=linux-gnu-c-opt Satish On Wed, 28 Jul 2010, Xuan YU wrote: On Jul 28, 2010, at 5:07 PM, Jed Brown wrote: On Wed, 28 Jul

[petsc-users] PETSC with-debugging and petsc errors

2010-08-02 Thread Satish Balay
On Mon, 2 Aug 2010, John-Michael Fischer wrote: I was wondering if when PETSC is compiled with-debugging=1 (or any other configure flag for that matter) I can get it to report the line and filename that errors occur on in my source. It's dutifully reporting the file and line # of the

[petsc-users] Compile and link correctly using conf/petscvariables and conf/variables

2010-08-05 Thread Satish Balay
The idea here is - you would: include conf/variables in your makefile and use PETSC_LIB in the target - so everything works correctly. The complexity you mention is to manage --with-single-lib=0/1. So Instead of trying to list libraries directly in your makefile [which becomes unportable] - I

[petsc-users] Compile and link correctly using conf/petscvariables and conf/variables

2010-08-05 Thread Satish Balay
You configure script would have to create this makefile and echo ${PETSC_LIB} - and other stuff you need form a given install of PETSc - to detect what you want. Anythings else is not portable. [And won't work across different PETSc installations - of the same version] If you want a simple

<    1   2   3   4   5   6   7   8   9   10   >