<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body 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="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="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>
</body>
</html>