ments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/5b805a6c/attachment.htm>
o do this.
Thanks.
Khalid
Here is the code:
PetscViewerBinaryOpen(appctx->comm,"Pxbin_10.out",FILE_MODE_READ,&viewer);
PetscViewerBinaryGetDescriptor(viewer,&fd);
PetscBinaryRead(fd,header,1,PETSC_SCALAR);
?for (k=zs; khttp://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/83123661/attachment.htm>
nualpages/Vec/VecScatterCreate.html
).
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/54f8a5cf/attachment.htm>
On Aug 29, 2011, at 8:32 PM, Tabrez Ali wrote:
> Hello
>
> If during Vecscatter (say from a global vector x to a local vector y) if 'is'
> and 'iy' are such that almost all values being scattered happen to be on the
> local proc then almost 0 or very little MPI calls would be made internally.
Hello
If during Vecscatter (say from a global vector x to a local vector y) if
'is' and 'iy' are such that almost all values being scattered happen to
be on the local proc then almost 0 or very little MPI calls would be
made internally. Is this correct?
Such a setup would not cause scalability
On Aug 29, 2011, at 5:49 PM, David Fuentes wrote:
> Is there any reason why petsc compiled to link with goto blas shared
> libraries would not run multi-threaded by default ?
We don't do anything to prevent it from using multi-threaded. Te first thing
I would suggest is make sure it truly is
THANK YOU! that was indeed the case!
I switched to defaults and it now converges with one iteration, and it
does so with my jacobian as well :)
I'll go play now :)
Here it is :)
-- next part --
A non-text attachment was scrubbed...
Name: milan-debug-session.log
Type: application/octet-stream
Size: 29689 bytes
Desc: not available
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/a2d10f86/attachment-0001.obj>
: not available
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/1b0ef54b/attachment-0001.obj>
Is there any reason why petsc compiled to link with goto blas shared
libraries would not run multi-threaded by default ?
I've set (OMP/GOTO)_NUM_THREADS=8
but when I call dgemm from PETSc I can't seem to get it to run on
multiple cores (<= 100% cpu usage from top).
I checked and the test routines
It is 27*27 with 343 non-zeros, i use LU for all the linear equations.
On Mon, August 29, 2011 5:46 pm, Jed Brown wrote:
> On Mon, Aug 29, 2011 at 16:36, Likun Tan wrote:
>
>
>> If i keep using the uncoupled small system and use ksp to solve the
>> linear equations sequentially, there is nowher
> Or your system is somewhat ill-conditioned. If you add '-mat_fd_type ds' to
> the last command line, are the errors similar or improve a bit?
Norm of matrix ratio 3.24879e-07 difference 3.24879e-05
So I guess it is the same...
> A different way to test for a correct Jacobian is to run with
> -
Thank you very much.
Since i have millions of such linear equations to solve, so efficiency is
my big concern.
If i keep using the uncoupled small system and use ksp to solve the linear
equations sequentially, there is nowhere i can take advantage of parallel
computing.
I think of the coupled sy
our right hand sides for that process in a matrix B, and
call MatMatSolve().
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/2d6b6f9a/attachment-0001.htm>
Instead of solving Ax=b with different right-hand-side sequentially, we
can also form a sparse block diagonal matrix A and a vector b composed of
all the elements. Then we can set values to each section of b concurrently
and solve the enlarged system in parallel, is this an efficient way?
And also
me matrix on
every processor and do a direct solve with all the right hand sides.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/b89854da/attachment.htm>
u want anyway.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/6ad212c0/attachment.htm>
>
>
--
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/2c6560b0/attachment-0001.htm>
TML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/d67ebe9c/attachment.htm>
Ok... I had a problem with indexing (-1 offset) in the local to global
mapping and now I get the following:
Testing hand-coded Jacobian, if the ratio is
O(1.e-8), the hand-coded Jacobian is probably correct.
Run with -snes_test_display to show difference
of hand-coded and finite difference Jacobia
ct type to KSPSetTolerances()? (Fortran doesn't warn
about this, sadly.)
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/b57166d6/attachment.htm>
hments/20110829/91d28087/attachment.htm>
log.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/92b9e931/attachment.htm>
Thanks,
I read ex16.c, it seems that the two linear equations are solved
sequentially?
On Mon, August 29, 2011 10:52 am, Satish Balay wrote:
> On Mon, 29 Aug 2011, Matthew Knepley wrote:
>
>
>> On Mon, Aug 29, 2011 at 2:40 PM, Likun Tan
>> wrote:
>>
>>
>>> Dear all,
>>>
>>>
>>> I have linear s
Dear all,
I have linear systems Ax=b with different right-hand-side vectors. These
vectors are independent. Could i generate b simultaneously and solve Ax=b
simultaneously?
best,
Likun
robustness is still a
problem and making efficient software for it is more tricky and AFAIK, has
not been done.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/05ebb9b6/attachment.htm>
On Mon, 29 Aug 2011, Matthew Knepley wrote:
> On Mon, Aug 29, 2011 at 2:40 PM, Likun Tan wrote:
>
> > Dear all,
> >
> > I have linear systems Ax=b with different right-hand-side vectors. These
> > vectors are independent. Could i generate b simultaneously and solve Ax=b
> > simultaneously?
> >
>
HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/62696564/attachment.htm>
Did you run with -snes_test_display it will display the matrices computed
both ways so you can see which entries are different and maybe track down the
problem.
Barry
On Aug 29, 2011, at 6:39 AM, Milan Mitrovic wrote:
> Ok... I had a problem with indexing (-1 offset) in the local to gl
Hello Gong Ding,
thank you very much for sharing your work on FVM support for libmesh.
I have not been able to find this part in libmesh. Is it available as a
separate patch to the libmesh sources?
Thanks
Marek
- Original Message -
From: Gong Ding
To: Marek Schmitt ; PETSc users list
30 matches
Mail list logo
