Thanks Barry, it works fine! :)
Eric
On 2022-08-31 16:28, Barry Smith wrote:
You can use PetscObjectSetOptionsPrefix((PetscObject) partitioner,...
On Aug 31, 2022, at 10:21 AM, Eric Chamberland
wrote:
Anyone have something to say about this?
Thanks,
Eric
On 2022-07-06 22:32, Eric
Anyone have something to say about this?
Thanks,
Eric
On 2022-07-06 22:32, Eric Chamberland wrote:
Hi,
I was looking for PetscPartitionerSetOptionsPrefix function, but only
found this post mentioning it.
I was "blindly" coding with in mind the possibility to manage
different o
can
easily and low-error prone call C, any other viable candidates?
--
Eric Chamberland, ing., M. Ing
Professionnel de recherche
GIREF/Université Laval
(418) 656-2131 poste 41 22 42
oner changeable from options.
Do you have anything against?
Vaclav
--
Eric Chamberland, ing., M. Ing
Professionnel de recherche
GIREF/Université Laval
(418) 656-2131 poste 41 22 42
I find the repository containing the data for these tests?...
Or what is the right value for DATAFILESPATH?
Thanks,
Eric
--
Eric Chamberland, ing., M. Ing
Professionnel de recherche
GIREF/Université Laval
(418) 656-2131 poste 41 22 42
for taking news! :)
Eric
On 2021-03-30 1:47 p.m., Barry Smith wrote:
Eric,
How are things going on this OpenMP front? Any bug fixes from
hypre or SuperLU_DIST?
BTW: we have upgraded to OpenMPI 4.1 perhaps this resolves some
issues?
Barry
On Mar 22, 2021, at 2:07 PM, E
I added some information here:
https://github.com/xiaoyeli/superlu_dist/issues/69#issuecomment-804318719
Maybe someone can say more than I on what PETSc tries to do with the 2
mentioned tutorials that are timing out...
Thanks,
Eric
On 2021-03-15 11:31 a.m., Eric Chamberland wrote
ix it.
Barry
I got the problem to occur with ex56 with 2 MPI ranks and 4 OpenMP
threads, if I used less than 4 threads it did not generate an
indefinite preconditioner.
On Mar 14, 2021, at 1:18 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Hi again,
ok, just saw that some matrices have lines of "0" in case of 3D hermite
DOFs (ex: du/dz derivatives) when used into a 2D plane mesh...
So, my last problem about hypre smoother is "normal".
However, just to play with one of this matrix, I tried to do a "LU" with
mumps icntl_24 optio
would not have seen the bug.
Hopefully they can just fix it.
Barry
I got the problem to occur with ex56 with 2 MPI ranks and 4 OpenMP
threads, if I used less than 4 threads it did not generate an
indefinite preconditioner.
On Mar 14, 2021, at 1:18 PM, Eric Chamberland
<mailto:eric.cham
Reported timeout bugs to SuperLU_dist too:
https://github.com/xiaoyeli/superlu_dist/issues/69
Eric
On 2021-03-14 2:18 p.m., Eric Chamberland wrote:
Done:
https://github.com/hypre-space/hypre/issues/303
Maybe I will need some help about PETSc to answer their questions...
Eric
On 2021-03
tps://github.com/hypre-space/hypre/issues>
On Mar 14, 2021, at 3:44 AM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
For us it clearly creates problems in real computations...
I understand the need to have clean test for PETSc, but for me, it
reveals that hypre isn
test suite
for each test, you can add the extra flag in your command line
TIMEOUT=180 (default is 60, units are seconds).
Thanks, I’ll ping you on GitLab when I’ve got something ready
for you to try,
Pierre
On 12 Mar 2021, at 8:54 PM, Eric Chamberland
mailto
seconds).
Noted, I will add this to my scripts and the dockerfiles too...
Thanks, I’ll ping you on GitLab when I’ve got something ready for you
to try,
Thank you for all your work!
Eric
Pierre
--
Eric Chamberland, ing., M. Ing
Professionnel de recherche
GIREF/Université Laval
(418) 656-2131 poste 41 22 42
(sorry for the delay).
Thanks,
Eric
On 2021-03-11 9:03 a.m., Eric Chamberland wrote:
Hi Pierre,
ok, that's interesting!
I will try to build a docker image until tomorrow and give you the
exact recipe to reproduce the bugs.
Eric
On 2021-03-11 2:46 a.m., Pierre Jolivet wrote:
On 11 Mar 2021
counts/ksp_ksp_tutorials-ex5f_superlu_dist_2.counts
ok ksp_ksp_tutorials-ex5f_superlu_dist_2 # SKIP Fortran required for
this test
Barry
On Mar 10, 2021, at 11:03 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Barry,
to get a some follow up on --with-openmp=1 fail
Barry,
to get a some follow up on --with-openmp=1 failures, shall I open gitlab
issues for:
a) all hypre failures giving DIVERGED_INDEFINITE_PC
b) all superlu_dist failures giving different results with initia and
"Exceeded timeout limit of 60 s"
c) hpddm failures "free(): invalid next siz
ified MKL-Pardiso Intel-Windows32',
None,
[os.path.join(dir,libdir,'mkl_core.lib'),'mkl_intel_c.lib','mkl_intel_thread.lib']+i+['libiomp5md.lib'],'32','yes')
return
The assumption is that the link will fail unless the correct librar
linux2-c-opt-ompi/tests/counts/snes_tutorials-ex56_attach_mat_nearnullspace-0_bddc_approx_hypre.counts
ok snes_tutorials-ex56_attach_mat_nearnullspace-0_bddc_approx_hypre
ok diff-snes_tutorials-ex56_attach_mat_nearnullspace-0_bddc_approx_hypre
TEST
arch-linux2-c-opt-ompi/tests/counts/ksp_ks
Mar 2021, at 9:14 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Hi,
It all started when I wanted to test PETSC/CUDA compatibility for our
code.
I had to activate --with-openmp to configure with --with-cuda=1
successfully.
I then saw that PETSC_HAVE_OPENMP is
linking configuration for example, should it be a good thing to have
this (--with-openmp=1) tested into the pipelines (with external packages
of course)?
Thanks,
Eric
Thanks,
Pierre
On 2 Mar 2021, at 9:14 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Hi,
It a
e detailed information.
I have cc:ed Pierre so he can look at the HPDDM failures.
On Mar 2, 2021, at 2:14 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Hi,
It all started when I wanted to test PETSC/CUDA compatibility for
our code.
I had to activate --w
cc:ed Pierre so he can look at the HPDDM failures.
On Mar 2, 2021, at 2:14 PM, Eric Chamberland
<mailto:eric.chamberl...@giref.ulaval.ca>> wrote:
Hi,
It all started when I wanted to test PETSC/CUDA compatibility for our
code.
I had to activate --with-openmp to configure with
know):
https://giref.ulaval.ca/~cmpgiref/dernier_ompi/2021.03.02.02h16m01s_make_test.log
https://giref.ulaval.ca/~cmpgiref/dernier_ompi/2021.03.02.02h16m01s_configure.log
Thanks,
Eric
--
Eric Chamberland, ing., M. Ing
Professionnel de recherche
GIREF/Université Laval
(418) 656-2131 poste 41 22 42
Perfect!
I confirm the master branch is working again for me!
Thanks a lot! :)
Eric
On 2019-04-23 10:48 p.m., Stefano Zampini wrote:
Eric,
Is in maint and master now
On Apr 23, 2019, at 9:22 AM, Smith, Barry F. via petsc-dev
wrote:
Great, thanks Stefano for fixing them all
On Apr
Hi,
I just noticed there are warnings at compile time for some
mkl_cpardiso.c routines:
/pmi/cmpbib/compilation_BIB_gcc_redhat_petsc-master_debug/COMPILE_AUTO/petsc-master-debug/src/mat/impls/aij/mpi/mkl_cpardiso/mkl_cpardiso.c:
In function ‘MatDestroy_MKL_CPARDISO’:
/pmi/cmpbib/compilation_B
-clusters-and-hpc-technology/topic/797761
Eric
On 19/04/18 09:01 AM, Eric Chamberland wrote:
Hi,
this morning, mpich/master with PETSc is 100% working again for us.
Thanks to both commits:
https://github.com/pmodels/mpich/commit/c597c8d79deea220a42751fda0f01ce70764c260
https://github.com
Hi,
since we recently upgraded our nighlty valgrind validations to PETSc
3.9.2, valgrind found a bug in PetscOptionsFindPair_Private, called from:
PetscOptionsHasName("M_AssembleurGD_H_/pmi/cmpbib/compilation_BIB_gcc_valgrind/COMPILE_AUTO/BIBTestValidation/Ressources/opt/Test.assemblageDynamiq
master branch soon.
Thanks,
Min
On 2018/04/17 14:10, Eric Chamberland wrote:
Hi,
are we talking about the "tag" passed to MPI_Isend for example?
but does that mean there is something to change for any MPI call which
involves tags usage or is it only a PETSc "bad" tag usage?
y wrote:
On Tue, 13 Mar 2018, Eric Chamberland wrote:
Hi,
each night we are testing mpich/master with our petsc-based code. I don't
know if PETSc team is doing the same thing with mpich/master? (Maybe it is a
good idea?)
Everything was fine (except the issue
https://github.com/pmodels
ease.
Thanks,
Eric
On 14/03/18 03:35 AM, Eric Chamberland wrote:
Hi,
fwiw, the actual mpich/master branch doesn't passes the PETSc "make
test" after a fresh installation... It hangs just afer the 1 MPI
process test, meaning it is locked into the 2 process
orials/ex19 run successfully with 1
MPI process
^Cmakefile:151: recipe for target 'test' failed
make: [test] Interrupt (ignored)
thanks,
Eric
On 13/03/18 08:07 AM, Eric Chamberland wrote:
Hi,
each night we are testing mpich/master with our petsc-based code. I
don't know if PETS
Hi,
each night we are testing mpich/master with our petsc-based code. I
don't know if PETSc team is doing the same thing with mpich/master?
(Maybe it is a good idea?)
Everything was fine (except the issue
https://github.com/pmodels/mpich/issues/2892) up to commit 7b8d64debd,
but since co
Hi,
it is not a big issue, but I just wanted to know if it has been noticed
that super-lu-dist 5.3.0 libs are now in PETSC_DIR/lib64 instead of
PETSC_DIR/lib like all other libs?
For me, after installation, it is the only package that have libs
installed there... I just added it to my LD_LI
On 12/02/18 12:11 PM, Jed Brown wrote:
Eric Chamberland writes:
why does mpi stuff appears in the first one and not the second?
I would prefer to have "-lmpi" in the second one too... When is it
decided to "hide" external packages "under" PETSc? At configura
Hi,
just to let know that since yesterday or before yesterday, I have errors
in mkl_cpardiso.c on petsc/master:
http://www.giref.ulaval.ca/~cmpgiref/petsc-master-debug/2018.01.31.02h00m01s_make.log
other logs:
http://www.giref.ulaval.ca/~cmpgiref/petsc-master-debug/2018.01.31.02h00m01s_confi
if useful.
Satish
*From: *Eric Chamberland <mailto:eric.chamberl...@giref.ulaval.ca>
*Sent: *Wednesday, September 13, 2017 8:41 AM
*To: *petsc-dev <mailto:petsc-dev@mcs.anl.gov>
*Subject: *[petsc-dev] ParMETIS change
Hi,
having trouble to compile our code with PETSc-master new default
Hi,
having trouble to compile our code with PETSc-master new defaults
options for PARMEtis.
Parmetis was using double precision before...
By curiosity, is there a big performance gain to use float?
Thanks,
Eric
Hi,
thanks for these infos.
Isn't there a bug into the macros?
Let's say I retreived a copy of petsc-master 3.7.2 then I have:
#define PETSC_VERSION_RELEASE0
#define PETSC_VERSION_MAJOR 3
#define PETSC_VERSION_MINOR 7
#define PETSC_VERSION_SUBMINOR 2
Now an API changed in offi
Hi,
Juste compiled with last night tarball of petsc-master and found that
MatGetSubMatrix has changed to MatCreateSubMatrix.
Fine, but I can't use the petscversion.h to switch into my code with
#ifdef over any of:
#define PETSC_VERSION_MAJOR 3
#define PETSC_VERSION_MINOR 7
#define
On 18/01/17 02:32 PM, Barry Smith wrote:
Go for it.
Ok, here it is! :)
https://bitbucket.org/petsc/petsc/pull-requests/618/bug-fixe-the-debug_mode-for-abusive/diff
It also includes an optimisation...
Eric
On 18/01/17 11:26 AM, Stefano Zampini wrote:
You’re right. Maybe we should use MatFindZeroDiagonals if missing is true.
I may try to do my first pull request if no one is under the way to
correct this
Is that ok?
Thanks,
Eric
Hi Stefano,
On 18/01/17 11:07 AM, Stefano Zampini wrote:
Perhaps this is a valid fix
- if (missing)
SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_WRONGSTATE,"Matrix is missing
diagonal entry in row %D",d);
+ if (missing) {
+for (i=0; i
I think that d only contains the first row with
Hi,
we are calling MatZeroRows_SeqAIJ (petsc 3.7.5) with:
N=3
rows[0]=15
rows[1]=30
rows[2]=65
diag = 1
x=0
b=0
and (Mat_SeqAIJ*)A->data->keepnonzeropattern is true.
but MatMissingDiagonal_SeqAIJ finds a missing diagonal at line 1 and
then return on error here:
if (missing) SETERRQ1(PETSC
Hi,
shouldn't there be a MPI_Reduce in MatMissingDiagonal_MPIAIJ so it is
really collective as stated in documentation?
Thanks,
Eric
Excellent!
Thanks a lot! :)
Eric
Le 2017-01-06 à 19:37, Barry Smith a écrit :
Added in master
On Jan 6, 2017, at 2:46 PM, Eric Chamberland
wrote:
Hi,
we are trying to add a "petsc" developpement into our code but we need to call
MatSeqSBAIJGetArray_SeqSBAIJ for i
Thanks!!! :)
Eric
On 06/01/17 04:14 PM, Barry Smith wrote:
Run ./configure with the additional option --with-visibility=0
We'll add a MatSeqSBAIJGetArray()
Barry
On Jan 6, 2017, at 2:46 PM, Eric Chamberland
wrote:
Hi,
we are trying to add a "petsc" developpemen
Hi,
we are trying to add a "petsc" developpement into our code but we need
to call MatSeqSBAIJGetArray_SeqSBAIJ for it.
Unlike MatSeqAIJGetArray_SeqAIJ which is accessible with either
MatSeqAIJGetArray
or
ierr =
PetscUseMethod(A,"MatSeqAIJGetArray_C",(Mat,PetscScalar**),(A,array));CHKERRQ(
On 16/12/16 01:06 PM, Stefano Zampini wrote:
What about raising an error from the interface if mat_reuse_matrix is
requested with the same input and output matrix?
Me too, I was asking myslef if this shall be done!
Eric
the reuse flag on a matrix which is not
of type SeqDense. Can you provide a MWE to reproduce the problem?
Thanks,
Stefano
2016-12-16 16:39 GMT+03:00 Eric Chamberland :
Another precision: the problem may still be there since dec 15, not tonight...
Eric
On 16/12/16 08:30 AM, Eric Chamberland
Ok, the offending call made is:
ierr = (*conv)(mat,newtype,reuse,M);CHKERRQ(ierr);
with:
print newtype
$1 = (MatType) 0x7fdc7bf8414a "seqdense"
print reuse
$2 = MAT_REUSE_MATRIX
and mat is:
(gdb) print *mat
$4 = {hdr = {classid = 1211214, bops = {{getcomm = 0x7fdc6adda05b
, view = 0x7fdc6b
Another precision: the problem may still be there since dec 15, not
tonight...
Eric
On 16/12/16 08:30 AM, Eric Chamberland wrote:
Hi,
something seems broken with master branch since last night.
With our nighlty tests, we now have many times this PETSC ERROR:
[0]PETSC ERROR: Cannot reuse
Hi,
something seems broken with master branch since last night.
With our nighlty tests, we now have many times this PETSC ERROR:
[0]PETSC ERROR: Cannot reuse matrix of type (null)
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
for trouble shooting.
[0]PETSC ERROR: Pets
Hi,
is there an error in code documentation in pvecimpl.h:
typedef struct {
VECHEADER
PetscIntnghost; /* length of local portion
including ghost padding */
isn't "nghost" the length of ghost padding only?
Thanks,
eric
with m,n declared "int", for me 2 warnings arise from this line:
first, "write" takes a size_t argument:
write(fd,vals,(size_t)(m*n)*sizeof(double))
and returns a "ssize_t" value, so the rest of the line compiles fine
like this:
!= (ssize_t)((size_t)(m*n)*sizeof(double)
which gives:
write(
Le 2016-09-12 à 22:29, Barry Smith a écrit :
It's funny we've been doing this for almost 20 years and no one has noticed
before
but we are only 11 days after clang 3.9.0 release... :)
Eric
Le 2016-09-12 à 17:32, Satish Balay a écrit :
Do you get these warnings with PETSc library build aswell?
I can't tell since I didn't tried to build PETSc with clang...
The logging code tries to log all messages in library and in
application - and prints a summary with -info.
You can disable
Hi,
I just discovered that all my calls to MPI are replaced by PETSc #define
macros...
I don't really care, but I now have this warning with clang 3.9.0:
/pmi/cmpbib/compilation_BIB_gcc_redhat_squash/COMPILE_AUTO/GIREF/src/commun/Parallele/PABroadcast.icc:324:41:
warning: possible misuse of
On 28/08/16 06:24 PM, Jed Brown wrote:
Eric Chamberland writes:
Le 2016-08-27 à 15:22, Jed Brown a écrit :
Barry Smith writes:
Well in Eric's case he knows he needs a bug fix in sub minor 3 and wants to
generate an error if linked against an earlier one that does not have th
Le 2016-08-27 à 15:22, Jed Brown a écrit :
Barry Smith writes:
Well in Eric's case he knows he needs a bug fix in sub minor 3 and wants to
generate an error if linked against an earlier one that does not have the bug
fix. There needs to be a way for users to indicate if they need above
Hi,
Le 2016-08-27 à 15:07, Barry Smith a écrit :
On Aug 27, 2016, at 1:59 PM, Jed Brown wrote:
What about
#define PetscInitialize(argc,argv,file,help)
PetscInitialize_Private(argc,argv,file,help,PETSC_VERSION_MAJOR,PETSC_VERSION_MINOR,PETSC_VERSION_SUBMINOR,PETSC_VERSION_RELEASE)
and Petsc
Hi,
wanting to prevent a possible error: I compile and link (dynamically) my
code with petsc-3.7.2 BUT I may change my LD_LIBRARY_PATH to point to
3.7.0...
I would the like my code to stop and tell me I am using the wrong libs...
I could add something like:
if (PETSC_VERSION_MAJOR != PetscI
Hi guys,
something is broken for configuring superlu_dist with
petsc-master.tar.gz available at 02h00 am:
http://www.giref.ulaval.ca/~cmpgiref/petsc-master-debug/2016.03.22.02h00m01s_configure.log
http://www.giref.ulaval.ca/~cmpgiref/petsc-master-debug/2016.03.22.02h00m01s_RDict.log
I also t
Hi Barry,
the configuration of hypre looks broken since last night...?
Here is the log issued from our automatic compilation of
petsc-master.tar.gz at 2016-02-28 02h00:
http://www.giref.ulaval.ca/~cmpgiref/petsc-master-debug/configure_20160228_0200.log
Everything was fine until last night...
Le 2016-01-27 11:49, Barry Smith a écrit :
Sorry about that. Fixed in master.
Thanks for the quick fix!!! :)
Eric
Hi,
since dec 15, we compile and test our code with the nightly
"petsc-master" tarball
(http://ftp.mcs.anl.gov/pub/petsc/petsc-master.tar.gz)
and everything was fine until this morning.
I noticed 4 tests in error:
0]PETSC ERROR: - Error Message
--
On 21/10/15 03:36 PM, Barry Smith wrote:
Sure, but in that "case" the count would be identical on both processes
anyways and so would not detect a problem? How could the count be different?
hmmm, sorry, I think I am wrong and you are right.
It won't stop the code earlier neither detect a d
On 21/10/15 03:02 PM, Barry Smith wrote:
On Oct 21, 2015, at 1:54 PM, Eric Chamberland
wrote:
On 21/10/15 02:21 PM, Barry Smith wrote:
Thanks. I don't think I need the Icount but. I think I can just generate an
error if the FILE and LINE don't match exactly?
Barry, I highly
On 21/10/15 02:21 PM, Barry Smith wrote:
Thanks. I don't think I need the Icount but. I think I can just generate an
error if the FILE and LINE don't match exactly?
Barry, I highly recommend you to use the static counter too otherwise
you won't catch the following bug:
if all processes s
SomewhatPetscValidLogicalCollectiveInt(__FILE__,comm);
sorry, should be something like:
...
SomewhatPetscValidLogicalCollectiveInt(function_that_sum_chars_to_int(__FILE__),comm);
...
Eric
On 21/10/15 11:36 AM, Barry Smith wrote:
On Oct 21, 2015, at 10:28 AM, Eric Chamberland
wrote:
On 21/10/15 10:56 AM, Barry Smith wrote:
I similarly saw the same behavior in the debugger that you reported and was
mystified.
This is why I called it a nasty bug. The setting of
On 21/10/15 10:56 AM, Barry Smith wrote:
I similarly saw the same behavior in the debugger that you reported and was
mystified.
This is why I called it a nasty bug. The setting of different nonew flag
meant that one process called the Allreduce() at the end of
MatAssemblyEnd_MPIAIJ()
i
Hi Barry,
if I may suggest, is it possible to have this test case added as a
nightly test since, I don't know but... a matrix with no local lines
may be something not widely tested...???
Thanks!
Eric
On 20/10/15 10:47 PM, Barry Smith wrote:
Eric,
Thanks for the test case. I have d
Thanks Barry! :)
another question: while trying to understand this, I used the
"-on_error_attach_debugger ddd" which worked, BUT the line where it was
breaking all seemed ok to me. I mean, the value of "op" variable tested
by at line petsc-3.6.2/src/mat/interface/matrix.c:5264 :
5256 PetscE
SET_2nd_PROC_TO_HAVE_NO_LOCAL_LINES
to see that it runs fine if the 2nd process have local lines...
The "magic" happens only if I call MatConvert on the "original" (C)
matrix...
Hope this helps!
Thanks,
Eric
On 08/10/15 04:18 PM, Barry Smith wrote:
On Oct 8, 20
On 08/10/15 03:30 PM, Satish Balay wrote:
The commented code is sequential code - and shouldn't make a
difference.
...but it does!
Perhaps your application has other issues.
Can you verify if your code is valgrind clean?
http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
Yes
I think the same problem is in MatShift_MPIBAIJ.
I mean if I have a processor with no (local) lines at all in a matrix
(which is my bugging case here), it should not call Mat*SetPreallocation
at that point...
That sounds ok to you?
thanks,
Eric
On 08/10/15 02:45 PM, Eric Chamberland wrote
preallocated) {
ierr = MatMPIAIJSetPreallocation(Y,1,NULL,0,NULL);CHKERRQ(ierr);
} /*else if (!aij->nz) {
ierr = MatSeqAIJSetPreallocation(maij->A,1,NULL);CHKERRQ(ierr);
}*/
ierr = MatShift_Basic(Y,a);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
Eric
On 08/10/15 09:45 AM, Eric Chamberland
Hi!
hint: For petsc 3.6.1, if I apply the patches that Patrick gave (in the
head of this thread), it all works correctly,
Eric
On 07/10/15 10:11 PM, Eric Chamberland wrote:
Hi Jed,
If I recall correctly, that patch is only checking for incorrect use by
ensuring that the flags are set
Hi Jed,
If I recall correctly, that patch is only checking for incorrect use by
ensuring that the flags are set collectively. If your code is correct
then it doesn't affect you and if your code is incorrect, it will error
now where as it would behave incorrectly before. The implementation was
d
/petsc-3.6.2_opt_openmpi-1.10_mkl_mt/lib/libpetsc.so.3.6(MatShift_MPIAIJ+0xf6)
[0x7fa4eb90da45]
#012:
/opt/petsc-3.6.2_opt_openmpi-1.10_mkl_mt/lib/libpetsc.so.3.6(MatShift+0x9e)
[0x7fa4eb6dbf94]
while the other continued...
thanks again!
Eric
On 07/10/15 03:02 PM, Eric Chamberland wrote:
Hi
Hi Barry,
just compiled/tested with 3.6.2.
On a 2 processes example, the non-debug version is hanging silently and
indefinitely but in debug mode, it now abort when calling MatShift when
a process has no entries:
#0 0x7fffdbc6c065 in raise () from /lib64/libc.so.6
#1 0x7fffdbc6d4e8
On 06/08/2015 04:54 PM, Satish Balay wrote:
Please do check src/docs/website/documentation/changes/dev.html and
update it with any missing changes entries..
http://www.mcs.anl.gov/petsc/documentation/changes/dev.html
thanks,
Satish
On Mon, 8 Jun 2015, Barry Smith wrote:
If the nightly te
On 03/25/2015 11:10 PM, Barry Smith wrote:
On Mar 25, 2015, at 7:59 PM, Eric Chamberland
wrote:
Hi Barry,
On 03/25/2015 07:25 PM, Barry Smith wrote:
This is an issue with any use of "nested" matrices (of which MPIBAIJ is a
special, two level case). One way we could handle
Hi Barry,
On 03/25/2015 07:25 PM, Barry Smith wrote:
This is an issue with any use of "nested" matrices (of which MPIBAIJ is a
special, two level case). One way we could handle this is (I think universally) would be
to introduce a unique PETSC_ERR_MAT_ENTRY_NONZERO error flag and then whe
On 09/08/2014 06:08 PM, Jed Brown wrote:
Barry Smith writes:
When should we install external packages as shared, currently if
—with-PACKAGENAME-shared is set or if —with-shared-libraries is set
we do it for some packages (like MPICH). Other packages, such as
hypre don’t build p
On 11/04/2014 04:40 PM, Barry Smith wrote:
Unfortunately the test only checks if the link works, it doesn’t test
if the prefetch actually does anything or might generate an error.
Writing a proper (runtime) test for this, especially in a batch
environment seems very difficult.
Ok.
#2- How c
Hi,
it all looks like we have a bug with prefetch with PGI 14.7 + cpu
specific compilation options. We have a "signal 4 : illegal instruction
(not reset when caught) " within MatMult_SeqAIJ_Inode. The only thing
we see there than can bug is the prefetch thing...
We were asking ourself two
On 04/29/2014 03:14 PM, Jed Brown wrote:
Eric Chamberland writes:
> You're supposed to use MatGetLocalSubMatrix(). Translating from global
> indices to local indices is a disaster that we want to avoid. So we go
> the other way. Speak the language of "split l
On 04/27/2014 12:32 PM, Jed Brown wrote:
> Barry Smith writes:
>
>>I am totally confused and likely users are to. Please point me to an
>>example that actually uses Matnest the way it is “supposed to be
>>used”.
>
> src/snes/examples/tutorials/ex28.c
>
> Note that this does not contai
On 04/17/2014 03:47 PM, Jed Brown wrote:
Eric Chamberland writes:
Now I noticed another thing while measuring performances (notice: I used
a PETSc compiled with debugging=yes... please don't flame me...):
The first assembly for a 549252 square matrix takes me 4.505s but if i
sort the in
On 04/16/2014 06:00 PM, Barry Smith wrote:
Removed in barry/rm-null-values-to-matsetvalues and next
doh! ;-)
Ok, I have done it with my own vector of zeros...
Now I noticed another thing while measuring performances (notice: I used
a PETSc compiled with debugging=yes... please don't flam
On 04/16/2014 09:44 AM, Jed Brown wrote:
Eric Chamberland writes:
I was writing (new) code which do the firsts MatSetValues after the
MatXAIJSetPreallocation. We have to do this because the "real"
non-zeros will be added later by a mix of ADD_VALUES and
INSERT_VALUES... which would
On 04/12/2014 06:10 PM, Barry Smith wrote:
Patrick,
Currently it is supported for MatSetValues_SeqAIJ() and
MatSetValues_MPIAIJ(). It is not handled for BAIJ and SBAIJ, the code while
crash with segv in those cases (since we don’t check for a valid v pointer
before use :-). But su
Le 2013-03-12 16:57, Tim Tautges a ?crit :
>
> Also, I find that infographics here help a bunch, like
> http://mercurial.selenic.com/wiki/QuickReferenceCardsAndCheatSheets?action=AttachFile&do=get&target=Mercurial-Usage-v1.0bz-300dpi-color.png.
>
> I looked a few days ago and couldn't find a go
Hi,
On 10/04/2012 11:12 PM, Barry Smith wrote:
> supposedly in petsc-dev
I just want to know if it will be included in next release... the
problem still occurs with my freshly compiled 3-3.p5...
Thanks!
Eric
Jed, Hong,
thank you for your answers!
Eric
Le 12/06/2012 01:25 PM, Jed Brown a ?crit :
> Eric, we would welcome an implementation that batched together a
> run-time-specified number of RHS vectors and used MUMPS' native
> multiple-RHS routines. The problem is that there has to be enough
> me
Hi,
I just want to know if there is a known difficulty or unfeasible task if
someone would like to implement MatMatSolve_MUMPS?
Thanks,
Eric
Hi,
we have compile petsc-2.3.3-p15 with openMPI 1.3.4. To have it working
with our C++ code, I did a "#define OMPI_SKIP_MPICXX".
In petsc.h, this is already done but for MPICH, on line #137:
#define MPICH_SKIP_MPICXX 1
It would be nice to also have a:
#define OMPI_SKIP_MPICXX" 1
But I don't
Hi,
we are using PETSc in our code and we have a problem with, I think, the
ghosted values that we expect.
We developed our own pre-conditioner, and try to use it in a parallel
environment. With other precondtioners (PETSc built-in), everything
works fine. But with our home-made one, here is
1 - 100 of 101 matches
Mail list logo