<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Ferrand,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
We do not have support for PtAP with P in aij format and A in dense format yet. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">MatProductReplaceMats() is used only when a new matrix has the SAME non-zero pattern as the one to be replace, which is not your case. </span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important"><br>
</span></div>
<div style="font-size: 12pt; color: rgb(0, 0, 0);"><font face="Times New Roman, Times, serif">You may do:</font></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">B = A*P (MatMatMult)</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">C = Pt*B (MatTransposeMatMult)</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">in which, you simply set P as aij format and A as dense. These option should work. </span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">Let me know if this works or not. If it works, I can add this support to petsc.</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family:"Times New Roman", Times, serif;background-color:rgb(255, 255, 255);display:inline !important">Hong</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> petsc-users <petsc-users-bounces@mcs.anl.gov> on behalf of Ferrand, Jesus A. <FERRANJ2@my.erau.edu><br>
<b>Sent:</b> Wednesday, March 2, 2022 2:10 PM<br>
<b>To:</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> [petsc-users] Missing Function MatSeqAIJGetArray_C</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
Greetings:</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
I need to carry out a matrix product of the form P^tAP, where matrix P is always sparse, but matrix A can either be sparse or dense. To handle both cases, I define two (2) separate matrix products using the sequence.<br>
</div>
<div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
MatProductCreate()<br>
<div>MatProductSetType()</div>
<div>MatProductSetAlgorithm()<br>
</div>
<div>MatProductSetFill()</div>
<div>MatProductSetFromOptions()</div>
MatProductSymbolic()</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
And then when I need to carry out the multiplication I call.<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
MatProductReplaceMats()</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
MatProductNumeric()</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
When matrix A is sparse, the code runs flawlessly. When matrix A is dense, however, I get this error that says that there is a missing function "MatSeqAIJGetArray_C."</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
I have no clue as to why there could be a dependency on the type of my matrix A. I also tried MatPtAP() which essentially does the same as the "MatProduct" callouts I mentioned and the error is the same.<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
<div>[0]PETSC ERROR: No support for this operation for this object type</div>
<div>[0]PETSC ERROR: Cannot locate function MatSeqAIJGetArray_C in object</div>
<div>[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.</div>
<div>[0]PETSC ERROR: Petsc Release Version 3.16.0, unknown </div>
<div>[0]PETSC ERROR: ./par on a linux-c-dbg named F86 by jesus Wed Mar 2 14:55:05 2022</div>
<div>[0]PETSC ERROR: Configure options --with-32bits-pci-domain=1 --with-debugging =1 --download-ptscotch --download-metis --download-parmetis --download-chaco --download-hdf5</div>
<div>[0]PETSC ERROR: #1 MatSeqAIJGetArray() at /home/jesus/SAND/PETSc_install/petsc/src/mat/impls/aij/seq/aij.c:4550</div>
<div>[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[1]PETSC ERROR: No support for this operation for this object type</div>
<div>[1]PETSC ERROR: Cannot locate function MatSeqAIJGetArray_C in object</div>
<div>[1]PETSC ERROR: [0]PETSC ERROR: #2 KAFormKE2D() at /home/jesus/SAND/FEA/3D/PARALLEL.c:1845</div>
<div>[0]PETSC ERROR: #3 MESHTraverseDepth() at /home/jesus/SAND/FEA/3D/PARALLEL.c:411</div>
<div>[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.</div>
<div>[1]PETSC ERROR: Petsc Release Version 3.16.0, unknown </div>
<div>[1]PETSC ERROR: #4 KAFormK() at /home/jesus/SAND/FEA/3D/PARALLEL.c:1881</div>
<div>[0]PETSC ERROR: #5 main() at /home/jesus/SAND/FEA/3D/PARALLEL.c:2394</div>
<div>./par on a linux-c-dbg named F86 by jesus Wed Mar 2 14:55:05 2022</div>
<div>[1]PETSC ERROR: [0]PETSC ERROR: PETSc Option Table entries:</div>
<div>[0]PETSC ERROR: -benchmark_iter 2</div>
<div>Configure options --with-32bits-pci-domain=1 --with-debugging =1 --download-ptscotch --download-metis --download-parmetis --download-chaco --download-hdf5</div>
<div>[1]PETSC ERROR: #1 MatSeqAIJGetArray() at /home/jesus/SAND/PETSc_install/petsc/src/mat/impls/aij/seq/aij.c:4550</div>
<div>[0]PETSC ERROR: -mesh_filename L.msh</div>
<div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div>
<div>[1]PETSC ERROR: #2 KAFormKE2D() at /home/jesus/SAND/FEA/3D/PARALLEL.c:1845</div>
<div>[1]PETSC ERROR: #3 MESHTraverseDepth() at /home/jesus/SAND/FEA/3D/PARALLEL.c:411</div>
<div>application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0</div>
<div>[1]PETSC ERROR: #4 KAFormK() at /home/jesus/SAND/FEA/3D/PARALLEL.c:1881</div>
<div>[1]PETSC ERROR: #5 main() at /home/jesus/SAND/FEA/3D/PARALLEL.c:2394</div>
<div>[1]PETSC ERROR: PETSc Option Table entries:</div>
<div>[1]PETSC ERROR: -benchmark_iter 2</div>
<div>[1]PETSC ERROR: -mesh_filename L.msh</div>
<div>[1]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div>
<span>application called MPI_Abort(MPI_COMM_WORLD, 56) - process 1</span><br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div id="x_Signature">
<div>
<div style="font-family:"Times New Roman",Times,serif; font-size:12pt; color:rgb(0,0,0)">
</div>
<div></div>
<div></div>
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0px; margin-bottom:0px; margin-top:0; margin-bottom:0"></p>
<div>
<p style="margin-top:0px; margin-bottom:0px; margin:0px 0px 10.66px"></p>
<div>
<p style="background:white none repeat scroll 0% 0%; margin:0px 0px 11px; line-height:normal">
<span style="margin:0px; color:black; font-family:"Times New Roman",Times,serif; font-size:12pt">Sincerely:</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:red; font-family:"Courier New",monospace; font-size:14pt"><b style=""><span style="font-family:"Courier New",monospace">J.A. Ferrand</span></b></span><u><span style="margin:0px; color:red; font-family:"Swis721 BlkCn BT",sans-serif; font-size:12pt"></span></u></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt">Embry-Riddle Aeronautical University - Daytona Beach FL</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt">M.Sc. Aerospace Engineering | May 2022</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt">B.Sc. Aerospace Engineering</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt">B.Sc. Computational Mathematics</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:rgb(0,0,0); font-family:Courant; font-size:12pt"> </span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt"><span style="color:rgb(0,0,0); font-family:Courant; font-size:12pt; display:inline!important; font-size-adjust:none; background-color:rgb(255,255,255)">Sigma Gamma Tau</span><span style="color:rgb(0,0,0); font-family:Courant; font-size:12pt; display:inline!important; font-size-adjust:none; background-color:rgb(255,255,255)"> </span><br>
</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:rgb(0,0,0); font-family:Courant; font-size:12pt">Tau Beta Pi</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt"><span style="margin:0px; color:rgb(0,0,0); font-family:Courant; font-size:12pt">Honors Program</span></span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt"> </span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<u><span style="margin:0px; color:black; font-family:Courant; font-size:12pt">Phone:</span></u><span style="margin:0px; color:black; font-family:Courant; font-size:12pt"> (386)-843-1829</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal">
<u><span style="margin:0px; color:black; font-family:Courant; font-size:12pt">Email(s):</span></u><span style="margin:0px; color:black; font-family:Courant; font-size:12pt"> ferranj2@my.erau.edu</span></p>
<p style="background:white none repeat scroll 0% 0%; margin:0px; line-height:normal; text-indent:0.5in">
<span style="margin:0px; color:black; font-family:Courant; font-size:12pt"><span style="margin:0px">
</span>jesus.ferrand@gmail.com</span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>