Hi,
I've managed to compile and run the examples in MSDev when i use MPICH.
However, if I compile without MPI, I 'll get errors during linking. How
should I modify my configuration file in MSDev so that I can run without
mpi?
Thank you very much!
-- next part --
An HTML
, 25 Nov 2006, Ben Tay wrote:
Hi,
I'm trying to use MPICH2 for the mpi during compilation of my Petsc. It
works OK if I use MPICH (mpich.nt.1.2.5). I try to specify the
directory,lib
(mpich2.lib),include etc but it keeps on saying that mpi_init cannot
work or
something like that. May I
.
On 11/25/06, Matthew Knepley knepley at gmail.com wrote:
You need to configure using --wih-mpi=0.
Matt
On 11/24/06, Ben Tay zonexo at gmail.com wrote:
Hi,
I've managed to compile and run the examples in MSDev when i use MPICH.
However, if I compile without MPI, I 'll get errors during
Hi,
I've some problems compiling PETSC on one of my school's server. I've
already compiled a working library on another linux system. Is it possible
for me to just copy the library from one to another. They are both redhat
linux system and I think they're running on xeon processors.
What are the
Hi,
I'm new to mpi. If I'm coding my code to use PETSC to solve my linear
eqn in parallel, must I compile in mpif90 (for fortran) or can I
compile in ifort and then instruct to send the job to n processors?
Thank you.
Hi,
I have a very simple fortran code. It compiles on a 32bit system with mkl
with no errors but on em64t, it gives undefined reference to error.
It works when I compiled with the supplied blas/lapack. However if I use
Intel mkl, it gives the error as stated above.
My code is
global.F
Yes it ran successfully. I've attached the output.
thank you very much.
On 1/11/07, Satish Balay balay at mcs.anl.gov wrote:
Do PETSc examples work?
Send us the output from
make test
Staish
On Wed, 10 Jan 2007, Ben Tay wrote:
Hi,
I have a very simple fortran code. It compiles
main.f90
Also - looks like you are compiling linking at the same step.. [i.e
using .o and .f90 in the same command]. Sugest crearing .o files before
linking..
i.e
${FLINKER} -o a.out global.o main.o ${PETSC_LIB}
Satish
On Thu, 11 Jan 2007, Ben Tay wrote:
Yes it ran successfully. I've
Hi,
I've the sequential version of PETSc working and I have some questions about
the parallel implementation.
My code is to solve the NS eqn. Mainly it has to solve 2 linear eqns -
momentum poisson. The poisson eqn solving takes the most % of time so it
going to be solved in parallel. Most
/examples/tutorials/ex29.c and ex22.c
Good luck,
Barry
On Sat, 13 Jan 2007, Ben Tay wrote:
Hi,
I've the sequential version of PETSc working and I have some questions
about
the parallel implementation.
My code is to solve the NS eqn. Mainly it has to solve 2 linear eqns
.
PS: Has the make clean on MPICH been updated on
BuildSystem/config/packages/MPI.py?
On 1/13/07, Ben Tay zonexo at gmail.com wrote:
Thanks Barry Aron.
I've tried to install mpich2 on a scratch directory and it finished in a
short while.
On 1/13/07, Barry Smith bsmith at mcs.anl.gov
hi,
i install PETSc using the following command:
./config/configure.py --with-vendor-compilers=intel --with-gnu-compilers=0
--with-blas-lapack-dir=/lsftmp/g0306332/inter/mkl/lib/32
--with-mpi-dir=/opt/mpich/intel/ --with-x=0 --with-shared
then i got:
Compilers:
C Compiler:
pan
--- Ben Tay zonexo at gmail.com wrote:
hi,
i install PETSc using the following command:
./config/configure.py --with-vendor-compilers=intel
--with-gnu-compilers=0
--with-blas-lapack-dir=/lsftmp/g0306332/inter/mkl/lib/32
--with-mpi-dir=/opt/mpich/intel/ --with-x=0
MPI job. So please try using the script.
Barry
On Tue, 16 Jan 2007, Ben Tay wrote:
Hi Pan,
I also got very big library files if I use PETSc with mpich2.
Btw, I have tried several options but I still don't understand why I
can't
get mpi to work with PETSc.
The 4 processors
Hi,
I can't get it to work with mpich2 (pre-built) although someone else told me
he got mpich2 working with PETSc. however his mpich2 was compiled from
source.
on the other hand, i got mpich1 (pre-built) to work and built with PETSc.
hope it helps.
On 1/12/07, Tan Meng YUE tmyue at
Hi,
I've encountered some MPI programming problems. Does anyone know of a
mailing list sort of discussion forum similar to this PETSc one which deals
with MPI?
Btw, my question is as follow, it will be great if someone can help too ;-)
It's in fortran but I guess it should be similar in C. Put
Hi,
In my MPI CFD code, I have to solve 2 linear eqns - momentum and poisson
eqn. There's some error in the code and I don't know which one causes the
error. May I know if it is possible to solve 1 eqn in sequential mode while
the other in MPI to pinpoint the error since my code works ok in a
Hi,
I was discussing with another user in another forum (cfd-online.com) about
using PETSc in my cfd code. I am now using KSP to solve my momentum and
poisson eqn by inserting values into the matrix. I was told that using PETSc
this way is only for unstructured grids. It is very inefficient and
Hi,
My school's server only has MKL installed on some nodes. Hence I was told
that if I need to use it, I'll have to link them statically. So how does
PETSc use MKL?
Btw, I'm currently using shared PETSc library.
Thank you.
-- next part --
An HTML attachment was
.
Is there some other directory where they are usually install?
I'm wondering if I can use lapack + gotoblas together? I tried to specify
the libraries but ./configure failed.
Thank you.
On 1/27/07, Satish Balay balay at mcs.anl.gov wrote:
On Sat, 27 Jan 2007, Ben Tay wrote:
My school's server has
Hi,
May I know what's the difference between the optimized and debugging
version?
Have the optimization flags such as -O3 been used in the compilcation of the
libraries for the optimized one? Whereas for the debugging version, there is
no optimization by the compiler?
If my program works fine
to functions etc in the debug version.
Barry
On Mon, 29 Jan 2007, Matthew Knepley wrote:
On 1/29/07, Ben Tay zonexo at gmail.com wrote:
Hi,
May I know what's the difference between the optimized and debugging
version?
Have the optimization flags such as -O3 been used
Hi,
somone suggested that I treat that face as a dirichlet boundary condition.
after 1 or a few iterations, the face value will be updated and it will be
repeated until covergerence. I wonder if that is possible as well?
It'll make the job much easier, although the iteration may take longer...
Hi,
I'm still relatively new to PETSC but I've finally managed to solve my own
poisson eqn using it in about a week's time.Thank you for writing such a
good solver.
However, I have some questions still about PETSC:
1. Are petscscaler the same as real nos.? I am writing in double precision
or
Hi,
regarding qn 4 ie the requirement of VecAssemblyBegin,VecAssemblyEnd, I
used VecPlaceArray, VecResetArray to change the values. I did not use
Assembly after the 1st time step and it works ok. So that is fine too,
right?
Thank you!
On 12/4/06, Satish Balay balay at mcs.anl.gov wrote:
On
Its normally better to change the values in place
So you mean I should use VecSetValue/s instead?
Thank you.
Regards.
On 12/4/06, Matthew Knepley knepley at gmail.com wrote:
On 12/4/06, Ben Tay zonexo at gmail.com wrote:
Hi,
regarding qn 4 ie the requirement of
VecAssemblyBegin
Hi,
I just tried but it's the same error. Any other suggestion?
thanks!
On 12/5/06, Matthew Knepley knepley at gmail.com wrote:
You have a broken Python installation. You can work around this using
the options --useThreads=0.
Thanks,
Matt
On 12/5/06, Ben Tay zonexo at gmail.com
own copy - and try with it.
configure python with:
./configure --prefix=${HOME}/python-2.4
make
make install
And configure PETSc with:
${HOME}/python-2.4/bin/python config/configure.py
Satish
On Wed, 6 Dec 2006, Ben Tay wrote:
Hi,
I just tried but it's the same error. Any other
wrote:
This is wierd. I don't think we've seen this error before. What kind
of linux is this running? I'm guessing there is a problem on the
system-software side.. [with threads]
Satish
On Wed, 6 Dec 2006, Ben Tay wrote:
Hi,
sorry, it still can't work. Installation of Python is ok
Hi,
I have been using a few different Krylov linear solver package such as
nspcg,sparsekit and now petsc to solve the linear eqns for my NS solver
momentum and poisson eqn. For moving bodies simulation, while solving the
poisson eqn, sometimes the result suddenly changes drastically, although the
Hi,
I tried to build the libraries on my school's em64t machines which is
supposed to be faster. The compilers are intel cc,fc ver 9. mkl is 7.21.I'm
trying to build a shared library. My command is
./config/configure.py --with-vendor-compilers=intel
--with-blas-lapack-dir=/opt/intel/mkl721/
Hi Satish,
Can you explain a bit how to change the RDict.py? I don't understand what
you mean.
Thank you.
On 12/7/06, Satish Balay balay at mcs.anl.gov wrote:
On Thu, 7 Dec 2006, Satish Balay wrote:
I'll sugest using a different machine for install.
Can you try the following change to
Hi,
I have some problems with using mpi. My code works with one server but when
I test it on another system, I get the wrong answers. In that system, the
mpi is located at /usr/local/topspin/mpi/mpich/. I tried to specify it with
--with-mpi-dir but it wasn't able to find during config. I then add
if the correct lib file is used.
Thanks.
On 12/13/06, Barry Smith bsmith at mcs.anl.gov wrote:
On Wed, 13 Dec 2006, Ben Tay wrote:
Hi,
I have some problems with using mpi. My code works with one server but
when
I test it on another system, I get the wrong answers. In that system
Hi,
I have some questions regarding the use of PETsc...
1. Is Intel MKL much faster than the downloaded BLAS/LAPACK? Or is it true
only for really large problem? Is ATLAS a good alternative too?
2. Can shared library be use in win32? I tried using --with-shared but after
./configure it still
Hi,
I'm using the downloaded mpich2 in the linux64 server. However, it is taking
up quite a lot of space. May I know which are the really important files
which I must keep? Which are the redundent files? I'm running out of disk
space
Thank you very much.
Regards
-- next part
Hi,
I'm now trying to modify my source code to enable parallel computing. I
tried to use the tutorial example and I managed to get it working. However,
in the example, there are a lot of ifs statement which is not efficient.
In my CFD code, I've subroutine meant for corner cells, north, south,
Hi,
I've used --download-hypre=1 to integrate hypre into PETSc. I've found the
hypre directory in externalpackages. hypre takes up about 136mb. Can I
delete the files other than the library files (such as src) after I've
compile them?
Thanks alot!
-- next part --
An HTML
Hi,
I tried to compile PETSc with compaq visual fortran in cygwin. My command is
$ ./config/configure.py --with-cc='win32fe cl' --with-fc='win32fe f90'
--download-f-blas-lapack=1 --with-f90-interface=win32 --with-x=0
--with-mpi-dir=/cygdrive/c/program\ files/MPICH/SDK
I then got this error
Hi,
I've managed to use PETSc in Compaq visual fortran using
--download-blas-lapack=1 with mpi installed.
Now I tried to recompile PETSc with --with-mpi=0. I then tried to remove
mpich.lib from the linking options. compiling is ok but during linking, I
got these error msg:
linking...
is a stub
library for one process) libmpiuni.lib
Barry
On Mon, 5 Mar 2007, Ben Tay wrote:
Hi,
I've managed to use PETSc in Compaq visual fortran using
--download-blas-lapack=1 with mpi installed.
Now I tried to recompile PETSc with --with-mpi=0. I then tried to remove
mpich.lib from
Hi,
I've compiled PETSc on my school's server and I wonder if I can simply copy
the compiled shared library/include files to my home comp's linux. Both uses
redhat.
Is there anything else I have to take care?
Thank you
-- next part --
An HTML attachment was scrubbed...
Hi,
i'm trying to solve my cfd code using PETSc in parallel. Besides the linear
eqns for PETSc, other parts of the code has also been parallelized using
MPI.
however i find that the parallel version of the code running on 4 processors
is even slower than the sequential version.
in order to find
Hi,
I've tried to use log_summary but nothing came out? Did I miss out
something? It worked when I used -info...
On 2/9/07, Lisandro Dalcin dalcinl at gmail.com wrote:
On 2/8/07, Ben Tay zonexo at gmail.com wrote:
i'm trying to solve my cfd code using PETSc in parallel. Besides the
linear
:
Impossible, please check the spelling, and make sure your
command line was not truncated.
Matt
On 2/9/07, Ben Tay zonexo at gmail.com wrote:
ya, i did use -log_summary. but no output.
On 2/9/07, Barry Smith bsmith at mcs.anl.gov wrote:
-log_summary
On Fri, 9 Feb
Ya, that's the mistake. I changed part of the code resulting in
PetscFinalize not being called.
Here's the output:
-- PETSc Performance Summary:
--
/home/enduser/g0306332/ns2d/a.out on a linux-mpi named
--with-mpi-dir=/opt/mpich/myrinet/intel/
-
On 2/10/07, Ben Tay zonexo at gmail.com wrote:
Hi,
I tried to use ex2f.F as a test code. I've changed the number n,m from 3
to 500 each. I ran the code using 1 processor and then with 4 processor. I
then repeat
eqn... wonder if hypre/BoomerAMG etc would help...
On 2/11/07, Lisandro Dalcin dalcinl at gmail.com wrote:
On 2/10/07, Ben Tay zonexo at gmail.com wrote:
In other words, for my CFD code, it is not possible to parallelize it
effectively because the problem is too small?
Is these true
Hi,
I'm trying to experiment with using external solvers. I have some questions:
1. Is there any difference in speed with calling the external software from
PETSc or directly using them?
2. I tried to install MUMPS using --download-mumps. I was prompted to
include --with-scalapack. After
.
Thank you.
On 2/12/07, Satish Balay balay at mcs.anl.gov wrote:
On Mon, 12 Feb 2007, Ben Tay wrote:
Hi,
I'm trying to experiment with using external solvers. I have some
questions:
1. Is there any difference in speed with calling the external software
from
PETSc or directly using
problem with hypre - remove
externalpackages/hypre-1.11.1b and retry.
Satish
On Mon, 12 Feb 2007, Ben Tay wrote:
Hi,
I tried to compile PETSc again and using --download-hypre=1. My command
given is
./config/configure.py --with-vendor-compilers=intel
--with-blas-lapack-dir=/lsftmp
Hi,
I have been using PETSc with visual fortran/intel mkl/mpich installed. This
has the same configuration as the configuration file .dsw supplied by PETSc.
However, now using another of my school's computer, MKL and MPICH are not
installed.
Is there anyway I can still use compaq visual fortran
Hi,
I've built static PETSc library with mpi/hypre w/o problems. The test e.g.
also worked.
When I tried to built my own a.out, using -static (fortran), the error is
lgcc_s not found. removing it in the make file resulted in a lot of error
for mpich2.
btw, i do not have root access. the
Hi,
I tried to compile PETSc and there's no problem. I also used the
*--with-batch=1
*option since my server is using a job scheduler. There's no shared library
and I'm installing hypre too.
However, when after sending the job for ex1f, I got the error msg:
Can't read MPIRUN_HOST.
So what's
, Ben Tay wrote:
Hi,
I tried to compile PETSc and there's no problem. I also used the
*--with-batch=1
*option since my server is using a job scheduler. There's no shared
library
and I'm installing hypre too.
However, when after sending the job for ex1f, I got the error msg:
Can't
It's ex1f. Thanks
On 5/22/07, Satish Balay balay at mcs.anl.gov wrote:
Is this a PETSc example? Which one is it?
Satish
On Tue, 22 May 2007, Ben Tay wrote:
ok. I sorted out that problem. Now I got this error:
*** glibc detected *** double free or corruption (!prev
:
call PCSetType(pc,'hypre',ierr)
call PCHYPRESetType(pc,'boomeramg',ierr)
Satish
On Tue, 22 May 2007, Ben Tay wrote:
Hi,
presently I am using ./a.out -pc_type hypre -pc_hypre_type boomeramg to
use
hypre.
However, I have 2 matrix to solve and now I only want one of them to use
]
Satish
On Tue, 22 May 2007, Ben Tay wrote:
It's ex1f. Thanks
On 5/22/07, Satish Balay balay at mcs.anl.gov wrote:
Is this a PETSc example? Which one is it?
Satish
On Tue, 22 May 2007, Ben Tay wrote:
ok. I sorted out that problem. Now I got this error
Ops how careless. Thanks again.
On 5/22/07, Barry Smith bsmith at mcs.anl.gov wrote:
pcsettype
On Tue, 22 May 2007, Ben Tay wrote:
Sorry I got this error during linking:
petsc_sub.o(.text+0x20c): In function `petsc_sub_mp_petsc_solver_pois_':
/nfs/home/enduser/g0306332/ns2d_c
have used :
call PCHYPRESetType(pc,'boomeramg',ierr)
call PCSetType(pc,'hypre',ierr)
I changed the 2 order but it still can't work.
Thanks again
On 5/22/07, Ben Tay zonexo at gmail.com wrote:
Ops how careless. Thanks again.
On 5/22/07, Barry Smith bsmith at mcs.anl.gov wrote
?
Thanks
On 5/22/07, Satish Balay balay at mcs.anl.gov wrote:
you are probably doing this inside a loop. check the code.
satish
On Tue, 22 May 2007, Ben Tay wrote:
Sorry to trouble you ppl again.
The 1st step worked but subsequently I got the error msg
Ok thank you very much
On 4/9/07, Barry Smith bsmith at mcs.anl.gov wrote:
On Mon, 9 Apr 2007, Ben Tay wrote:
Hi,
I am solving a poisson eqn and using hypre as a preconditioner has
greatly
increased the speed. I am trying to do some more tuning and I have some
questions:
1
Hi,
I've finished compiling PETSc and it use the fblaslapack external package.
May I know if I can delete the directories blas1-3, lapack1-6 and the
lapack-extra1-4?
Thanks
-- next part --
An HTML attachment was scrubbed...
URL:
Hi,
someone was talking abt core 2 duo performance on os x in some previous
email. it seems that due to memory issues, it's not possible to get 2x the
performance. there's also some mention of amd vs intel dual core.
for computation using PETSc, is there any reason to buy one instead of the
to do this is to manualy edit bmake/$PETSC_ARCH/petscconf and
change PCC_LINKER to use 'win32fe f90'
Satish
On Fri, 1 Jun 2007, Ben Tay wrote:
Hi,
I'm using ms vc++ and vf in windows xp. I tried to compile PETSc 2.3.3p0
.
Running configure has no problem. Compiling seems ok
Hi,
I found this link when doing a PETSc search.It is 30mins online lecture
which introduces and explains the use of PETSc on solving the laplace eqn.
Hope it's useful to new users
http://infostream.rus.uni-stuttgart.de/lec/307/537/real/
-- next part --
An HTML attachment
Hi,
I got this error during configure: Fortran could not successfully link C++
objects
My command is :
./config/configure.py --with-vendor-compilers=intel
--with-blas-lapack-dir=/usr/local/intel/mkl721cluster/ --with-debugging=0
--download-hypre=1 --with-x=0
Hi,
I need to insert values into a matrix and then solve this poisson eqn as
part of my cfd code. For efficiency sake, I simply use
do j=1,size_y
do i=1,size_x
matsetvalue
end do
end do
to insert values into all the cells
Then for specific cells, I need to enter some other
Smith bsmith at mcs.anl.gov wrote:
On Sun, 10 Jun 2007, Ben Tay wrote:
Hi,
I need to insert values into a matrix and then solve this poisson eqn as
part of my cfd code. For efficiency sake, I simply use
do j=1,size_y
do i=1,size_x
matsetvalue
end do
end do
Hi,
I am doing CFD and at each time step, I need to interpolate values at
different points. I'm doing bi-quadratic interpolation and hence to get each
interpolated value, I need to solve a 9x9 matrices. There's about 100+ pts
and hence 100+ 9x9 matrices to solve per time step.
I'm wondering if
= 1.0e-8
Thanks again
On 6/17/07, Barry Smith bsmith at mcs.anl.gov wrote:
Ben,
On Sat, 16 Jun 2007, Ben Tay wrote:
Hi,
I am trying to use PETSc to solve my poisson eqn. In the past, when I'm
not
using MPI, I set KSPSetTolerances to 1.0e-5. This seems to be good
enuff.
Now, I'm
Is the matrix symmetric? Do you expect/want it to be symmetric?
Barry
On Mon, 18 Jun 2007, Ben Tay wrote:
Oh btw, here's the output when KSPSetTolerances = 1.0e-5.
1:
0 KSP preconditioned resid norm 5.509130713952e+01 true resid norm
3.411770553001e+01 ||Ae||/||Ax
Hi,
I'm using ms vc++ and vf in windows xp. I tried to compile PETSc 2.3.3p0.
Running configure has no problem. Compiling seems ok but in the last part
when it tried to test examples, the error is
--Error detected during compile or link!---
See
Hi,
I am also trying to solve the NS eqns in 2d, finite volume mtd, using
fractional step formulation.
I am currently using KSPBCGS and PCILU to solve the momentum eqn matrix
and HYPRE's AMG and the default ksptype to solve the poisson eqn. Do you
ppl think this is a good choice? Or is there
Hi,
I'm formulating a NS solver which computes the flow past an airfoil. The
airfoil has a c-grid and it overlaps with the background cartesian grid. In
other words, it is a chimera or overset grid application.
The boundary cells of the c-grid are connected to the cartesian grid thru
Hi,
I'm thinking of ways to estimate and compare the performance of PETSc on
different CPUs. I think it will also enable one to make the wise choice
of whether to upgrade or not.
Of cos, the best way is to run your own code on the new machine to see
how much increase there is. However, most
Hi,
I am trying to solve a pressure poisson eqn resulting from the NS eqn.
How can I fixed the pressure at some nodes while solving the poisson eqn?
I tried to make the coefficient of the diagonal location of the node
very big e.g. 1e100. I also use
call
Hi,
I used 2 packages to solve my poisson eqn, which is part of my NS
unsteady solver. One is the NSPCG solver package which uses the
compressed row format to store the A matrix. The other is PETSc. I found
that using both solvers gave me similar answers for a number of time
step. However,
converging. Run with
-ksp_converged_reason and or -ksp_monitor. My guess is that the matrix is
becoming very ill-conditioned so the solvers with the default options are
failing.
Barry
On Fri, 3 Aug 2007, Ben Tay wrote:
Hi,
I used 2 packages to solve my poisson eqn, which is part
2007, Ben Tay wrote:
Hi,
I realised that the matrix storage format of NSPCG is actually ITPACK storage
format. Anyway, I tried to insert the values into a PETSc matrix, use MatAXPY
with the scalar a = -1 and then use MatNorm on the output matrix.
Using NORM_1, NORM_FROBENIUS
right?
Any other suggestion?
Thanks.
Barry Smith wrote:
Run with -ksp_type gmres -pc_type lu -ksp_monitor_true_residual
-ksp_truemonitor
At the bad iteration verify if the true residual is very small.
Barry
On Sun, 5 Aug 2007, Ben Tay wrote:
Well, ya, they're the same
Hi,
I've no problem writing out codes in fortran fixed format with PETSc.
However, is it possible to do it in fortran free format as well?
I'm using visual fortran and there's error.
original :
test.F
module global_data
implicit none
save
#include include/finclude/petsc.h
#include
subroutine xyz()
use foobar
implicit none
#include include/finclude/petsc.h
other includes
code
end subroutine
Satish
On Mon, 6 Aug 2007, Ben Tay wrote:
Hi,
I've no problem writing out codes in fortran fixed format with PETSc.
However,
is it possible to do it in fortran free
Hi,
I am trying to insert values into a matrix. I used to use INSERT_VALUES
with MatSetValues or MatSetValue. However, since I need to add values to
the same II,JJ location more than 1 time, I have to use ADD_VALUES with
MatSetValues.
However I got this error message:
Object is in wrong
be removed
from here..
other includes
#undef PETSC_AVOID_DECLARATIONS
moudle foobar
other module stuff
end module
Satish
On Mon, 6 Aug 2007, Ben Tay wrote:
Hi,
I tried to use
#define PETSC_AVOID_DECLARATIONS
#include include/finclude/petsc.h
other includes
#undef
variable declarations (in case you use a
PETSc type variable).
Good Luck,
Paul
Ben Tay wrote:
Hi,
I can compile and run with visual fortran by putting all the
declarations inside a subroutine e.g.
module ada
...
subroutine ddd
#define PETSC_AVOID_DECLARATIONS
#include
Hi,
when I changed my global.F to global.f90, I was told that
PETSC_COMM_SELF, PETSC_NULL_INTEGER and PETSC_NULL_CHARACTER does not
have a type.
I declare
MPI_Comm PETSC_COMM_SELF and it works
So, what about the other 2? You mention to use
integer NORM_MAX
parameter NORM_MAX=3
What values
Hi,
My original matrix is stored in Ellpack-Itpack format, which is used by
the NSPCG solver. I have problems when I insert the values of the matrix
into a PETSc matrix. I guessed I made some mistakes but I can't find
where during debugging. So I intent to just convert from the
Ellpack-Itpack
a type, and must have an explicit
type. [PETSC_NULL_CHARACTER]. Same for PETSC_NULL_INTEGER and
PETSC_COMM_SELF.
It was ok when I'm programming in fixed format.
Thanks
Satish Balay wrote:
On Thu, 9 Aug 2007, Ben Tay wrote:
Hi,
when I changed my global.F to global.f90,
Any
, everything's ok. You can give it a try too.
Lastly, I understand that you have been helping a lot of users and
there's really a lot of stupid users around (I know I'm one of them ;-)
). Thank you for your patience.
Satish Balay wrote:
On Thu, 9 Aug 2007, Ben Tay wrote:
Hi,
Guess I'm
can comiile locally.
Satish
On Thu, 9 Aug 2007, Ben Tay wrote:
Hi Satish,
I am sorry if I have made you angry. I may have missed out some important
points which you tried to highlight. I was converting several of my .F fixed
format files to .f90 (or .F90) free format and some files
2007, Ben Tay wrote:
Thank you Satish. I followed what you recommended exactly and used
sed -e 's/^C/\!/g'
With the new mpif.h, I got the following error:
global.i90
d:\cygwin\codes\MPICH\SDK\include\mpif.h(174) : Error: DEC$ ATTRIBUTES
DLLIMPORT attribute can appear only
Hi,
I am currently using the following way to solve my poisson eqn
in global.F (global variables)
implicit none
save
MatA_mat
Vecxx,b_rhs
KSP ksp
PCpc
PCType ptype
KSPTypeksptype
then they are initialized
call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
call
to the constructor are using 0-offset numbering?
I suspect that it would be 1-offset for both ia_spar and ja_spar which would
cause problems.
Matt
On 8/14/07, Ben Tay zonexo at gmail.com wrote:
Hi,
I am currently using the following way to solve my poisson eqn
in global.F (global variables
Well, I did call MatDestroy, but that's at the very last step when
everything's over. It's ok. I'll give to figure out where the problem
occurs.
Thanks
Matthew Knepley wrote:
On 8/14/07, Ben Tay zonexo at gmail.com wrote:
Oh, I will check up valgrind since I've not used it b4
Hi,
I found that there's a multigrid option in the manual which can be used
by adding -pctype mg -pcmgtype PC_MG_MULTIPLICATIVE.
However, it seems that I am actually just using LU since with or without
this option, I got exactly the same answer. Am I using multigrid the
wrong way?
Thanks
. Or if you have a
structured grid, consider using DMMG.
Thanks,
Matt
On 8/15/07, Ben Tay zonexo at gmail.com wrote:
Hi,
I found that there's a multigrid option in the manual which can be used
by adding -pctype mg -pcmgtype PC_MG_MULTIPLICATIVE.
However, it seems that I am
using DMMG.
Thanks,
Matt
On 8/15/07, Ben Tay zonexo at gmail.com wrote:
Hi,
I found that there's a multigrid option in the manual which can be used
by adding -pctype mg -pcmgtype PC_MG_MULTIPLICATIVE.
However, it seems that I am actually just using LU since with or without
it
lower?
Thank you.
Hong Zhang wrote:
On Thu, 16 Aug 2007, Ben Tay wrote:
Btw, is the file to download
ftp://ftp.mcs.anl.gov/pub/petsc/externalpackages/ml-5.0.tar.gz?
Yes. You can simply use '--download-ml=1' during configuration.
Run your program with '-help' to see all
also get the error msg:
*ML*WRN* This function has been compiled without the configure
*ML*WRN* option --with-ml_metis
*ML*WRN* I will put all the nodes in the same aggregate, this time...
*ML*WRN* (file ./Coarsen/ml_agg_METIS.c, line 954)
Thanks
Barry Smith wrote:
On Fri, 17 Aug 2007, Ben Tay
1 - 100 of 132 matches
Mail list logo