Dear All,
I hope you're having a nice day.
Compared the user manual of PETSc 3.1 and 3.4, I found that there are
several 'new' KSP types are added. Then, I have some confusions on which
KSP solver I need to pick. To my knowledge, if the matrix is symmetric,
positive definite, the CG can be used.
Dear All,
I hope you're having a nice day.
I'm trying to use BDDC as a preconditioner. The executable command is:
/mpiexec -np 8 ./ex45 -pc_type bddc -ksp_type cg -ksp_rtol 1.0e-7
-log_summary//
///However, it say:
/[//0]PETSC ERROR: - Error Message
I understand now. Thank for your answer.
best, :)
Alan
Alan zhenglun.wei at gmail.com writes:
Thanks, Dr. Knepley,
Can I treat AMG and MG like this way. To solve problems with
non-symmetric matrices, MG needs users to provide the coarse mesh
information. This take more complexity
Dear folks,
I hope you're having a nice day.
For the Poisson solver in /src/ksp/ksp/example/tutorial/ex45.c, I used
the ksp_type = CG to solve it before; it converges very fast with
pc_type = gamg. However, I was trying to check if the matrix generated
by the 'ComputeMatrix' is symmetric by using
grid
??: Sun, 31 Mar 2013 18:21:40 -0500
???:Zhenglun (Alan) Wei zhenglun.wei at gmail.com
???:Jed Brown jedbrown at mcs.anl.gov, Mark F. Adams
mark.adams at columbia.edu
Dear Dr. Brown,
Thank you so much for your answers. Here are my reply.
I reduced my code and attached here
Dear All,
I hope you're having a nice day.
Based on ksp ex45, a 3D Poisson solver with non-uniform grid is coded.
The PETSc options I used is:
/mpiexec -np 32 ./ex45 -pc_type gamg -ksp_type cg -pc_gamg_type agg
-pc_gamg_agg_nsmooths 1 -mg_levels_ksp_max_it 1 -mg_levels_ksp_type
richardson
Dear folks,
Here are several questions following up my previous ones:
1, I modified my code based on ex49.c. However, I found that if I don't
call 'DMDASetUniformCoordinates' before calling 'DMDAVecGetArray', then
an error comes up:
[0]PETSC ERROR: Null argument, when expecting valid
On 2/28/2013 6:31 PM, Matthew Knepley wrote:
On Thu, Feb 28, 2013 at 7:19 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear folks,
Here are several questions following up my previous ones:
1, I modified my code based on ex49.c
= dz = 0.2), the
mesh is shown in the picture 'doubleAll' in the attachment. It seems
that the mesh generated by the 'DMDASetUniformCoordinate' is always there.
best,
Alan
On 2/28/2013 7:25 PM, Matthew Knepley wrote:
On Thu, Feb 28, 2013 at 8:05 PM, Zhenglun (Alan) Wei
zhenglun.wei
Dear all,
I hope you're having a nice day.
I wonder if anyone could point out a example for me about the
DMDASetCoordinates(). I found this
(http://www.mcs.anl.gov/petsc/petsc-current/src/dm/impls/da/dacorn.c.html#DMDASetCoordinates)
online; however it does not really demonstrate how to create the
Dear folks,
Here I came across a problem.
[0]PETSC ERROR: Petsc has generated inconsistent data!
[0]PETSC ERROR: Have un-symmetric graph (apparently). Use
'-pc_gamg_sym_graph true' to symetrize the graph or '-pc_gamg_threshold
0.0' if the matrix is structurally symmetric.!
My code
than 4 processes (e.g. 6 or 8 processes).
Hope this helps to detect the problem.
thanks,
Alan
Original Message
Subject:A quick question on 'un-symmetric graph'
Date: Wed, 02 Jan 2013 15:41:47 -0600
From: Zhenglun (Alan) Wei zhenglun@gmail.com
To: PETSc
Dear folks,
I have a quick question on the DMDACreate3d.
In the manual, it says that the input format of this function is:
PetscErrorCode DMDACreate3d(MPI_Comm comm,DMDABoundaryType bx,DMDABoundaryType
by,DMDABoundaryType bz,DMDAStencilType stencil_type,PetscInt M,
are:
hg pull -u: under /petsc-dev
hg update
make all test
what else should I do? how can I check the PETSc version in my
computer?
thanks,
Alan
On Mon, Nov 26, 2012 at 2:16 AM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com
Thank you Sean,
I changed the
$ cd config/buildsystem
$ hg showconfig paths.default
to
https://bitbucket.org/petsc/buildsystem
It gives me 2 extra changes.
Although, it always giev me the warning of 'certificate not verified', it
updated lots of files. Should I just ignore
. It, then, shows that the 'petscdmmg.h' can not be found. Any
idea on this? :)
thanks,
Alan
On Sun, Nov 25, 2012 at 4:16 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com wrote:
Thank you Sean,
I changed the
$ cd config/buildsystem
$ hg showconfig paths.default
to
https://bitbucket.org
On Mon, Nov 26, 2012 at 12:59 AM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear Sean,
Thank you so much for your help. Here I have another question.
After updating the PETSc, I tried to 'make ex45' in
/src/ksp/ksp
changes. Since
you have many stale files around, you may want to use hg clean (this
will delete all untracked files!) to get back to a state like if you
had a fresh clone.
On Mon, Nov 26, 2012 at 1:25 AM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote
Dear Dr. Brown,
You're totally right. Thank you so much,
Alan
It tries to find the executable in PATH, maybe you need mpiexec -n
./app instead of mpiexec -n app.
On Mon, Nov 26, 2012 at 1:35 AM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote
Dear folks,
Here are some problems when I'm using the new PETSc and tried to
run src/ksp/ksp/example/tutorial/ex45.c
1, when I compile it, it says that 'DMSetInitialGuess' can not be found.
So, I commented it out; however, this command still appears in
Dear folks,
I came across an error while making a code:
makefile:7: *** missing separator. Stop.
This happens after I was tying to update the petsc-dev with some
error message; the configure.log, error.log and make.log are attached.
I've already sent it to petsc-maint list. Is this
something else?
thanks,
Alan
|
On Thu, 8 Nov 2012, Zhenglun (Alan) Wei wrote:
Dear folks,
I came across an error while making a code:
makefile:7: *** missing separator. Stop.
This happens after I was tying to update the petsc-dev with some error
message; the configure.log
output
please copy/paste the output of the following commands
hg parent
hg tip
hg status
hg pull
hg update
hg status
cd config/BuildSystem
hg parent
hg tip
hg status
hg pull
hg update
hg status
Satish
On Thu, 8 Nov 2012, Zhenglun (Alan) Wei wrote:
Thanks for your reply, Dr
'. If you
still have issues - send the new logs to petsc-maint.
Satish
On Thu, 8 Nov 2012, Zhenglun (Alan) Wei wrote:
*_hg parent_*
changeset: 24846:6e281c3f1f72
tag: tip
user:Matthew G Knepley knepley at gmail.com
date:Thu Nov 08 11:41:01 2012 -0500
summary
PETSC_ARCH
and retry [configure/make] - and then send us the logs.
Satish
On Thu, 8 Nov 2012, Zhenglun (Alan) Wei wrote:
Yes, it has some problem when I run
./$PETSC_ARCH/conf/reconfigure-$PETSC_ARCH.py
and similar problem when I do 'make'
I will send the log file to petsc-maint in 10 mins
Dear all,
I hope you're having a nice day.
I need some help to figure out what is wrong in my code. I had very
good performance after I change the preconditioner. The out_1, which is
attached, shows that. Last week, in order to use PETSc in out new
cluster, I reconfigure the PETSc
Dear folks,
I hope you're having a nice day.
I've been reading the thread of Enquiry regarding log summary
results and tried the similar things with -log_summary with my code.
However, my output seems totally different from Dr. TAY Wee-beng's. I
attached it here. My code adopts the
norm 1.800293604638e+01 ||r(i)||/||b|| 1.e+00
The finial question: what are the choice for the
/'/-mg_levels_ksp_type'. Are they as the same as 'ksp_type'?
thank you so much,
Alan
On 10/5/2012 12:00 PM, Jed Brown wrote:
On Fri, Oct 5, 2012 at 11:51 AM, Zhenglun (Alan) Wei
Dear Dr. Brown,
Thank you so much for your help. Your reply is really helpful and
knowledgeable.
sincerely,
Alan
On 10/5/2012 2:42 PM, Jed Brown wrote:
On Fri, Oct 5, 2012 at 2:33 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear Dr
Dear folks,
I was trying to check my current PETSc configuration. Is there any
command or file I can refer to? Specially, I want to know if my PETSc is
configured for distributed-memory or shared memory machine.
thanks,
Alan
Dear folks,
I did some test with -pc_type gamg with
/src/ksp/ksp/example/tutorial/ex45.c. It is not as good as default
-pc_type when my mesh (Cartisian) is 100*50*50; while it is a little bit
better than the default one when the mesh is 200*100*100. Therefore, I
guess this type of pc is
at gmail.com wrote:
On Fri, Sep 14, 2012 at 5:40 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear folks,
I did some test with -pc_type gamg with
/src/ksp/ksp/example/tutorial/ex45.c. It is not as good
on all of them.
It is generally better for performance to post the receives first, then
post the sends, then MPI_Waitall() on all the requests.
On Fri, Sep 7, 2012 at 4:50 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com wrote:
Dear folks,
I did more tests, since I want to make sure
:56 PM, Jed Brown wrote:
Numeric data that the solver sees should be stored in Vecs. You can
put other scalars in Vecs if you like.
On Sep 6, 2012 5:48 PM, Zhenglun (Alan) Wei zhenglun.wei at gmail.com
mailto:zhenglun.wei at gmail.com wrote:
Dear Dr. Brown,
I'm not quite
Dear All,
I hope you're having a nice day.
I met a memory problem for MPI data communication. I guess here is
a good place to ask this question since you guys are experts and may
experienced the same problem before.
I used the MPI derived data type (MPI_Type_contiguous,
of my code?
thanks,
Alan
On 9/6/2012 5:44 PM, Jed Brown wrote:
Are you familiar with VecScatter?
On Sep 6, 2012 5:38 PM, Zhenglun (Alan) Wei zhenglun.wei at gmail.com
mailto:zhenglun.wei at gmail.com wrote:
Dear All,
I hope you're having a nice day.
I met a memory
? Then likely it is an issue with your code, if
only one has problems then it is an MPI implementation issue.
Barry
On Sep 6, 2012, at 5:38 PM, Zhenglun (Alan) Wei zhenglun.wei at gmail.com
wrote:
Dear All,
I hope you're having a nice day.
I met a memory problem for MPI
Dear All,
I hope you're having a nice day.
I coded my program based on the
/src/ksp/ksp/example/tutorial/ex45.c. However, I met a problem that the
memory is accumulating while the program is running. I have a doubt on
the DMDAVecGetArray(da, b, bb) and DMDAVecRestoreArray(da,b,bb).
Thanks Barry. It helps a lot.
Alan
On 8/30/2012 6:18 PM, Barry Smith wrote:
On Aug 30, 2012, at 6:08 PM, Zhenglun (Alan) Wei zhenglun.wei at
gmail.com wrote:
Dear All,
I hope you're having a nice day.
I coded my program based on the /src/ksp/ksp/example/tutorial/ex45.c.
However
Dear All,
I hope you're having a nice day.
I'm trying to use PETSc to program a Poisson Solver with local
refinement grid (or adaptive mesh). Is that possible? I checked DMDA
functions and found DMDASetRefinementFactor. Is this what I need to look
into? or it is only for overall
Dear All,
I hope you're having a nice day.
I was trying to solve the Poisson equation with
/src/ksp/ksp/exmaples/tutorial/ex45.c; however, it always give me the
same result. I recalled the PetscLogStages and resolved this problem.
However, the process of reconstructing DM costs a lot
On 5/17/2012 1:11 PM, Matthew Knepley wrote:
On Thu, May 17, 2012 at 1:47 PM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear All,
I hope you're having a nice day.
I was trying to solve the Poisson equation with
/src/ksp/ksp
Dear All,
I hope you're having a nice day.
I wonder if there is any example in PETSc to solve 3D Poisson
equation with the periodic boundary condition in one direction. 2D is
fine if there is no 3D example.
thanks,
Alan
Thank you so much, Dr. Knepley.
Alan
On 4/19/2012 9:51 AM, Matthew Knepley wrote:
On Thu, Apr 19, 2012 at 10:49 AM, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Dear Dr. Knepley,
It is very nice to hear that. I will read the manual. Do
PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Fri, Apr 13, 2012 at 17:38, Zhenglun (Alan) Wei zhenglun.wei at gmail.com
wrote:
I have a final question on it. Is it taken a lot of memory for doing
this? As I understand, pointers won't occupy many memories and it works
like an alias
Dear All,
Any suggestions on a 3D Poisson Solver with both Dirichlet and
Neumann BC's?
Also, does PETSc have any example of Poisson solver (either 2D or
3D) with periodic boundary conditions.
thanks
Zhenglun
by doing this.
Is that true?
thanks,
Alan
On 4/12/2012 1:22 PM, Jed Brown wrote:
On Thu, Apr 12, 2012 at 13:16, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Thank you for your reply. It gives me another way to solve
this problem. However
Dear All,
I hope you're having a nice day.
I came across a problem while using MPI derived data type.
Although it does not related to PETSc, I hope you guys can give me some
suggestions or comments.
The reason I want to do this is that, in my PETSc code, I want to
send a set
On 4/12/2012 11:52 AM, Jed Brown wrote:
On Thu, Apr 12, 2012 at 11:37, Zhenglun (Alan) Wei
zhenglun.wei at gmail.com mailto:zhenglun.wei at gmail.com wrote:
Everything is fine in AlanFloatStruct.c. *However, when I was
trying to use 'calloc' to allocate memory to 'a[][]' rather
() line 51 in src/dm/impls/da/dagetarray.c
My code is attached.
thanks in advance,
Alan
On 4/7/2012 5:30 PM, Jed Brown wrote:
On Sat, Apr 7, 2012 at 17:28, Jed Brown jedbrown at mcs.anl.gov
mailto:jedbrown at mcs.anl.gov wrote:
On Sat, Apr 7, 2012 at 17:18, Alan Wei zhenglun.wei
Dear All,
I hope you're having a nice day.
I'm working on ex45. It does not have any output subroutine inside.
Therefore, I was trying to write out mine.
I was trying to use DMDAVecGetArray; but I found that the DM object
has been destroyed at very beginning of the main code.
Dear All,
I hope you're having a nice weekend.
I got a problem with /ksp/ksp/example/tutorial. I was trying to execute
this example with petsc-dev in a share memory dell workstation with 8
cores. The parameter for DMDACreate2d is 500, 500, which means it generates
a 1997*1997 2d grid. It
Thanks,
Alan
On Sat, Apr 7, 2012 at 3:49 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Apr 7, 2012 at 15:41, Alan Wei zhenglun.wei at gmail.com wrote:
I got a problem with /ksp/ksp/example/tutorial. I was trying to
execute this example with petsc-dev in a share memory dell
Dear All,
I have another questions on /src/ksp/ksp/example/tutorial/ex29.c. I
just install PETSc-dev in a distributed-memory. Everything is fine and it
finished 'make test'.
However, when I was trying to 'make' ex29.c, it comes this error
message:
makefile:18: /conf/variables: No such file
A little problem comes up. I commented out VTK part, the code still has the
same error message popping up.
Alan
On Sat, Apr 7, 2012 at 4:34 PM, Alan Wei zhenglun.wei at gmail.com wrote:
Thanks,
Alan
On Sat, Apr 7, 2012 at 3:49 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Apr 7
Equation solver in PETSc code?
thanks,
Alan
On Sat, Apr 7, 2012 at 5:12 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sat, Apr 7, 2012 at 17:10, Alan Wei zhenglun.wei at gmail.com wrote:
A little problem comes up. I commented out VTK part, the code still has
the same error message popping up
and ComputeMatrix (I assert some print-out
probes inside these two functions). However, If I use it again, it just does
not call the new one; rather, it calls the old one.
My code is attached. Could you please help me to fix this problem ^_^.
best,
Alan
On Sun, Oct 2, 2011 at 5:52 PM, Alan Wei
That's great, Matt. Thanks to point out the example for me.
best,
Alan
On Mon, Oct 3, 2011 at 12:06 PM, Matthew Knepley knepley at gmail.com wrote:
On Mon, Oct 3, 2011 at 12:00 PM, Alan Wei zhenglun.wei at gmail.com wrote:
Dear all,
I still have the problem. _
I can not find
Dear all,
I hope you're having a nice weekend. I'm still using this program to
solve Poisson Equations. After several discussions before, I know how to use
it in general. However, I met a problem when I want to use it twice to solve
a same Poisson Equations with two different RHS.(same
side. _ any suggestions on that?
best,
Alan
On Sun, Oct 2, 2011 at 3:56 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Oct 2, 2011, at 12:04 PM, Alan Wei wrote:
Dear all,
I hope you're having a nice weekend. I'm still using this program to
solve Poisson Equations. After several
In this example, i recalled the ComputeRHS.
Alan
On Sun, Oct 2, 2011 at 5:22 PM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Sun, Oct 2, 2011 at 16:31, Alan Wei zhenglun.wei at gmail.com wrote:
BTW, what is the DMMG interface you are talking about? is it DMMGSetKSP()?
Every object
Dear folks,
I just update petsc-dev following
PETSc can now be configured in the usual way, specified on the Installation
page http://www.mcs.anl.gov/petsc/petsc-as/documentation/installation.html.
To update your copy of petsc-dev
- cd petsc-dev
- hg pull -u
- cd config/BuildSystem
][i] - u[j-1][i] - u[j][i+1] - u[j][i-1] + rhs * dx*dy =
0
thanks,
Alan
On Wed, Sep 21, 2011 at 8:22 AM, Barry Smith bsmith at mcs.anl.gov wrote:
On Sep 20, 2011, at 10:47 PM, Alan Wei wrote:
Dear Dr. Smith,
I figure out this problem. Actually, I made my own RHS, but I did
, therefore -\Delta u = -1, which, then, will give
me a good result for Delta p = 1, is that true?
thanks in advance,
Alan
On Wed, Sep 21, 2011 at 11:19 AM, Matthew Knepley knepley at gmail.com wrote:
On Wed, Sep 21, 2011 at 4:16 PM, Alan Wei zhenglun.wei at gmail.com wrote:
However, why signs
bsmith at mcs.anl.gov wrote:
Alan,
It is very possible that the comment in the example code is wrong and has
an incorrect sign. If you tell us what example this is we'll check it and
fix the comment it if is wrong.
Barry
On Sep 21, 2011, at 11:28 AM, Alan Wei wrote:
Oh, Thanks
No problem, thank you guys to help me figure out the issue. ^_^
Alan
On Wed, Sep 21, 2011 at 1:26 PM, Barry Smith bsmith at mcs.anl.gov wrote:
The comment will be fixed in the next patch
Thanks for letting us know about the error
Barry
On Sep 21, 2011, at 1:19 PM, Alan Wei wrote
Dear folks,
I hope you guys are having a nice day.
I'm reading the /src/ksp/ksp/examples/tutorials/ex29.c.html and wonder
how to set up the convergence criteria? Currently I use it as a poisson
solver to solve a Poisson Equation with three direction Neumann BC's and one
direction
Thanks, Jed. Let me finish this part first and I will upgrade to version 3.2
to try new things.
best,
Alan
On Tue, Sep 13, 2011 at 1:10 AM, Jed Brown jedbrown at mcs.anl.gov wrote:
On Tue, Sep 13, 2011 at 05:59, Alan Wei zhenglun.wei at gmail.com wrote:
I'm still working on this poisson
Dear guys,
I hope you're having a nice day.
I'm still working on this poisson solver, and find a small problem.
I found from the original code that DMMGGetx(dmmg) can get the solution
vector after calculation. However, I want to have them as a format of
2-dimensional array (since I
Thanks, Matt. I thought there is a ghosted grid around Dirichlet and Neumann
BC's. ^_^
best,
Alan
On Fri, Sep 2, 2011 at 7:14 AM, Matthew Knepley knepley at gmail.com wrote:
On Thu, Sep 1, 2011 at 9:47 PM, Alan Wei zhenglun.wei at gmail.com wrote:
Dear All,
I hope you're having a nice
:20, Alan Wei zhenglun.wei at gmail.com wrote:
I thought there is a ghosted grid around Dirichlet and Neumann BC's.
Only if you ask for it with DMDA_BOUNDARY_GHOSTED (not available in
petsc-3.1). This example does not.
-- next part --
An HTML attachment was scrubbed
Dear All,
I hope you're having a nice day.
I'm still working on this Poisson solver example, and I still have a
little confusion on the boundary condition. I wonder if I have a 17*17 grid.
(index i = 0~16 and j = 0~16, solving for variable 'u')
1) If a Dirichlet BC applied on the left
Dear all,
I hope you're having a nice day.
I met a problem on makefile of PETSc. Original PETSc has the makefile in
the same directory of the .c files. However, I want to put all .c files to
another directory. I changed the makefile also; however it always has a
error saying the PETSc
${PETSC_SNES_LIB}
Also it generally is much easier to have make files that site in the same
directory as the source code, less hassle and less likely for things to go
wrong.
Barry
On Aug 25, 2011, at 1:56 PM, Alan Wei wrote:
Dear all,
I hope you're having a nice day.
I met
Dear all,
I hope you're having a nice day.
I'm using PETSc as a beginner. I found it is very easy to find a
definition of a function like 'DACreate2d' online and it explains very
clearly. However, I think it is really confusing for those explanation of
data structures like DMMG, DA, DM
at mcs.anl.gov wrote:
On Aug 18, 2011, at 9:40 PM, Alan Wei wrote:
Dear Barry,
That is really a great idea, and it is something that I want. I've
never noticed that the structure in C is so marvelous. However, I have a
further question. In that example, a variable x is defined like Field **x
Thanks Barry,
^_^
Alan
On Fri, Aug 19, 2011 at 3:51 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Aug 19, 2011, at 2:55 PM, Alan Wei wrote:
Dear Barry,
In ex19, the X is a kinda of global vec, while the x is a local vec.
Since DAVecGetArray makes the x and X sharing the same
Dear Barry,
Thanks for your reply.
best,
Alan
On Thu, Aug 18, 2011 at 5:49 PM, Barry Smith bsmith at mcs.anl.gov wrote:
On Aug 18, 2011, at 5:02 PM, Alan Wei wrote:
Dear all,
I hope you're having a nice day.
I'm trying to get the x, y-value of a coordinate created
Dear All,
I'm trying to 'draw' a 2-D cylinder in PETSc. It simply restore the x
and y coordinate values in two arrays. I wonder if I can use Vec in PETSc to
implement it (i.e. defining a vec circle) so that I can call its x or y
coordinate by circle.x or circle.y. Is this possible? also, is
].u = x[j][i].u;
f[j][i].v = x[j][i].v;
f[j][i].omega = x[j][i].omega - (x[j][i+1].v - x[j][i].v)*dhx;
f[j][i].temp = x[j][i].temp;
note that this means the values are stored interlaced in the vector.
Barry
On Aug 18, 2011, at 9:03 PM, Alan Wei wrote:
Dear All
at 11:22 AM, Alan Wei zhenglun.wei at gmail.com wrote:
Thanks, John. However, on my opinion, one of them has been commented out.
^_^
best,
Alan
On Thu, Jul 28, 2011 at 12:55 AM, John Chludzinski
jchludzinski at gmail.com wrote:
You have:
*typedef enum {DIRICHLET, NEUMANN} BCType
at 3:46 PM, Alan Wei zhenglun.wei at gmail.com wrote:
Dear John,
Thanks for your quick reply. I think the error message says that the
BCType has some problem. However, it is originally there and I have not
changed it at all. Is that because I define it in ComputeRHS again and
let
ex29
Dear Sir/Madam,
I hope you're having a nice day.
I asked several questions about src/ksp/ksp/example/tutorial/ex29.c
before, and learned some beginning staff of PETSc. Right now, I want to
start to use it to my own program. Firstly, I want to isolate the ComputeRHS
function in ex29, which
,
Alan
On Wed, Jul 27, 2011 at 2:04 PM, John Chludzinski jchludzinski at
gmail.comwrote:
Without wading through the attached zip file, it looks like you're missing
some typedef's that come from a missing include file. ---John
On Wed, Jul 27, 2011 at 2:38 PM, Alan Wei zhenglun.wei at gmail.com
Dear Sir/Madam,
Recently, I studied the /src/ksp/ksp/example/tutorial/ex29.c of solving
a 2D Poisson equation. It employs the Multigrid Method on a uniform
coordinate. I wonder if there is a way to solve the Poisson equation with a
nested grid in PETSc. Any suggestions on any example or
column file with, say MPI rank number
as the first column, can I do that? Could you please suggest a function to
call or I need to loop it by myself.
thanks in advance, ^_^
Alan
On Fri, Jun 10, 2011 at 5:03 PM, Alan Wei zhenglun.wei at gmail.com wrote:
Oh, I got it. Indeed I found the original
ooops, sorry, here it is.
^_^
Alan
On Fri, Jun 10, 2011 at 3:09 PM, Matthew Knepley knepley at gmail.com wrote:
On Fri, Jun 10, 2011 at 3:05 PM, zhenglun.wei at gmail.com wrote:
Dear Matt,
Following your suggestions, I modified the VecView_VTK in
/src/ksp/ksp/example/tutorial/ex29.c to
Dear Sir/Madam,
Lately, I read there is a optimized version of PETSc by using
'--with-debugging=0'. I have two questions here:
1, Do I need to use '/config/configure.py ... --with-debugging=0' when I
configure the PETSc in order to use the non-debugging version of PETSc?
Therefore, if I want
Hi guys,
I hope you're having a nice day.
I was trying to run ex29.c in src/ksp/ksp/example/tutorial. It is normal
to go through without MPI. However, with the same grid size, if I run it
with 4 processes with MPI, a error message of Caught signal number 11 SEGV:
Segmentation Violation,
89 matches
Mail list logo