<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Hello All,</p>
<p><br>
</p>
<p>My apologies, I have closed this e-mail window and the first
thing I read on the manual is "ALL processes that share a matrix
MUST call <a
href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatAssemblyBegin.html#MatAssemblyBegin">MatAssemblyBegin</a>()
and <a
href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatAssemblyEnd.html#MatAssemblyEnd">MatAssemblyEnd</a>()
the SAME NUMBER of times". So I understand that petsc simply does
not support unequal number of assembly calls.</p>
<p><br>
</p>
<p>My question evolves then. I have a problem at hand where I do not
know how many calls each process will make to MatAssembly
routines. Any suggestions to make this work?<br>
</p>
<br>
<div class="moz-cite-prefix">On 28-01-2018 22:09, Ali Berk Kahraman
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:60ec2ddf-fe37-4946-1a00-afdde34ceb2c@yahoo.com">Hello
All,
<br>
<br>
<br>
The code takes a parallel matrix and calls a function using that
matrix. That function fills the specified row of that matrix with
the id of the process that part of the matrix belongs in. You can
see the short code in the attachment, it is about 80 lines.
<br>
<br>
<br>
The problem is that the code gets into a deadlock at some point,
usually at the last row of each process except for the last
process (greatest pid). I use petsc with configure options
"--with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich
--download-fblaslapack".
<br>
<br>
<br>
I am a beginner with MPI, so I do not know what may be causing
this. My apologies in advance if this is a very trivial problem.
<br>
<br>
<br>
Best Regards to All,
<br>
<br>
<br>
Ali Berk Kahraman
<br>
<br>
M.Sc. Student, Mechanical Eng.
<br>
<br>
Bogazici Uni., Istanbul, Turkey
<br>
<br>
</blockquote>
<br>
</body>
</html>