MatMultHermitianTranspose could have a default implementation in terms of the others. Then it would work in all cases, right?<br><br><div class="gmail_quote">On Wed, Jul 18, 2012 at 7:59 AM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Now the SVD in SLEPc correctly applies complex conjugation in complex scalar builds. We have two options:<br>
- Implicit: uses MatMultHermitianTranspose<br>
- Explicit: uses MatTranspose plus MatConjugate (the default)<br>
<br>
The problem is that MatMultHermitianTranspose is defined only for SEQBAIJ (not for parallel BAIJ nor for AIJ sequential or parallel). So the implicit option is quite limited now for complex.<br>
<br>
Conversely, the explicit option works in all cases except for BAIJ matrices since MatConjugate is not defined for BAIJ.<br>
<br>
How difficult is to implement these operations?<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br>