<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Mangal;
panose-1:0 0 4 0 0 0 0 0 0 0;}
@font-face
{font-family:Mangal;
panose-1:0 0 4 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@SimSun";}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Texte de bulles Car";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.TextedebullesCar
{mso-style-name:"Texte de bulles Car";
mso-style-priority:99;
mso-style-link:"Texte de bulles";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hello everyone,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am trying to solve a A^T x = b system. For my applications, I had realized that the MatSolveTranspose does not converge properly. Therefore, I had implemented a MatTranspose followed by a MatSolve. This proved to converge perfectly (which
is strange since the transposed matrix has the same eigenvalues as the untransposed…). The problem is that for bigger matrices, the MatTranspose is very costly and thus cannot be used.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I tried using the superlu_dist package. Although it the package works perfectly for the MatSolve, I get the an “No support for this operation for this object type” error with MatSolveTransopse. I reproduced the error using the MatView
an ex10 tutorial. I can provide the matrix and rhs upon request. My command line was:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">ex10 -f0 A_and_rhs.bin -pc_type lu -pc_factor_mat_solver_package superlu_dist –trans<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">So it there an additional parameter I need to use for the transposed solve?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: No support for this operation for this object type<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Matrix type mpiaij<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Petsc Release Version 3.5.1, unknown <o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: /gpfs/fs1/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/ksp/examples/tutorials/ex10 on a ARGUS_impi_opt named hpc-user11 by ad007804 Thu Aug 28 16:41:15 2014<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: Configure options --CFLAGS="-xHost -axAVX" --download-hypre --download-metis --download-ml --download-parmetis --download-scalapack --download-superlu_dist --download-mumps --with-c2html=0 --with-cc=mpiicc --with-fc=mpiifort
--with-cxx=mpiicpc --with-debugging=yes --prefix=/gpfs/fs1/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/petsc-3.5.1 --with-cmake=/gpfs/fs1/aero/SOFTWARE/TOOLS/CMAKE/cmake-2.8.7/bin/cmake --with-valgrind=/gpfs/fs1/aero/SOFTWARE/TOOLS/PROGRAMMING/valgrind-3.9.0/bin/valgrind
--with-shared-libraries=0<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #1 MatSolveTranspose() line 3473 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/mat/interface/matrix.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #2 PCApplyTranspose_LU() line 214 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/pc/impls/factor/lu/lu.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #3 PCApplyTranspose() line 573 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/pc/interface/precon.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #4 KSP_PCApply() line 233 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/include/petsc-private/kspimpl.h<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #5 KSPInitialResidual() line 63 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/ksp/interface/itres.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #6 KSPSolve_GMRES() line 234 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/ksp/impls/gmres/gmres.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #7 KSPSolveTranspose() line 704 in /gpfs/fs2/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/ksp/interface/itfunc.c<o:p></o:p></p>
<p class="MsoNormal">[0]PETSC ERROR: #8 main() line 324 in /gpfs/fs1/aero/SOFTWARE/TOOLS/PROGRAMMING/petsc/src/ksp/ksp/examples/tutorials/ex10.c<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">FYI, the transpose solve is a typical application for adjoint optimization. There should be a big adjoint community of developers that try to solve the transposed matrix.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Any help is much appreciated,<o:p></o:p></p>
<p class="MsoNormal">Best,<o:p></o:p></p>
<p class="MsoNormal">Antoine<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F">Antoine DeBlois</span></b><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F"><br>
Specialiste ingenierie, MDO lead / Engineering Specialist, MDO lead<br>
Aéronautique / Aerospace</span><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#131E29"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F">514-855-5001, x 50862<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#131E29"><a href="mailto:antoine.deblois@aero.bombardier.com"><span style="color:blue">antoine.deblois@aero.bombardier.com</span></a><br>
<br>
</span><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F">2351 Blvd Alfred-Nobel</span></i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F"><br>
<i>Montreal, Qc<o:p></o:p></i></span></p>
<p class="MsoNormal"><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F">H4S 1A9</span></i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#7F7F7F"><br>
</span><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#131E29"><br>
<img border="0" width="162" height="34" id="Image_x0020_1" src="cid:image001.jpg@01CFC2DB.A38AE950" alt="Description : Description : http://signatures.ca.aero.bombardier.net/eom_logo_164x39_fr.jpg"><o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:8.5pt;font-family:"Arial","sans-serif";color:#555555;mso-fareast-language:EN-US">CONFIDENTIALITY NOTICE</span></b><span style="font-size:8.5pt;font-family:"Arial","sans-serif";color:#666666;mso-fareast-language:EN-US">
- This communication may contain privileged or confidential information.<br>
If you are not the intended recipient or received this communication by error, please notify the sender<br>
and delete the message without copying</span><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>