<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body>
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Excellent question. It may well be some quirk of my own setup that somehow makes VecView run rather sluggishly, making it hardly relevant to the general case.</p>
<p><br>
</p>
<p>- Peder</p>
</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>Fra:</b> Matthew Knepley <knepley@gmail.com><br>
<b>Sendt:</b> 4. juni 2021 16:48:05<br>
<b>Til:</b> Peder Jørgensgaard Olesen<br>
<b>Cc:</b> Jose E. Roman; petsc-users@mcs.anl.gov<br>
<b>Emne:</b> Re: [petsc-users] SLEPc: HDF5 support for SVD/EPSValuesView?</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">On Fri, Jun 4, 2021 at 10:28 AM Peder Jørgensgaard Olesen <<a href="mailto:pjool@mek.dtu.dk">pjool@mek.dtu.dk</a>> wrote:<br>
</div>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div dir="ltr">
<div id="gmail-m_8835617947864268106x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<p>What Matt suggests could perhaps be used as a workaround, though it appears neither elegant nor in my experience efficient. I may be wrong on this.</p>
</div>
</div>
</div>
</blockquote>
<div>Why is it not efficient?</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div dir="ltr">
<div id="gmail-m_8835617947864268106x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<p>My attempted solution was to use a binary viewer instead, but I can't find a clear way to retrieve values thus stored with PETSc, as these do not lend themselves to be read using VecLoad. Again one might go about that by wrapping them inside a vector and
send that to the viewer, which just brings us back around to the original problem.</p>
<p><br>
</p>
<p>- Peder</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_8835617947864268106x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Fra:</b> Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>><br>
<b>Sendt:</b> 4. juni 2021 16:01:04<br>
<b>Til:</b> Matthew Knepley<br>
<b>Cc:</b> Peder Jørgensgaard Olesen; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">
petsc-users@mcs.anl.gov</a><br>
<b>Emne:</b> Re: [petsc-users] SLEPc: HDF5 support for SVD/EPSValuesView?</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt">
<div>The problem is that here I am writing PetscReal's not PetscScalar's, and in EPS I am writing PetscComplex even in real scalars. That is why I did not implement it as you are suggesting.<br>
<br>
Jose<br>
<br>
<br>
> El 4 jun 2021, a las 15:56, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> escribió:<br>
> <br>
> On Fri, Jun 4, 2021 at 9:31 AM Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br>
> I could try and adapt the HDF5 code in PETSc for this, but I am no HDF5 expert. Furthermore I am busy at the moment, so it is faster if you could prepare a merge request with this addition.<br>
> <br>
> I think the easiest thing is to wrap a Vec around the values and just VecView() it into the same HDF5 file.<br>
> <br>
> THanks,<br>
> <br>
> Matt<br>
> <br>
> Jose<br>
> <br>
> <br>
> > El 4 jun 2021, a las 13:10, Peder Jørgensgaard Olesen via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> escribió:<br>
> > <br>
> > Hello<br>
> > <br>
> > In SLEPc one may write singular vectors to an HDF5 file using SVDVectorsView(), but a similar option doesn't seem to work for singular values using SVDValuesView(). The values are viewed correctly using other viewers, but nothing is seemingly produced when
using an HDF5 viewer. Looking at the source code seems to confirm this, and suggests that there is a similar situation with EPS.<br>
> > <br>
> > Is this correct, and if so, might there be plans to include such functionality in future releases?<br>
> > <br>
> > Med venlig hilsen / Best Regards<br>
> > <br>
> > Peder Jørgensgaard Olesen<br>
> > PhD Student, Turbulence Research Lab<br>
> > Dept. of Mechanical Engineering<br>
> > Technical University of Denmark<br>
> > Koppels Allé<br>
> > Bygning 403, Rum 105<br>
> > DK-2800 Kgs. Lyngby<br>
> <br>
> <br>
> <br>
> -- <br>
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
> -- Norbert Wiener<br>
> <br>
> <a href="https://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br>
<br>
</div>
</span></font></div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr" class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener</div>
<div><br>
</div>
<div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>