<div dir="ltr">On Wed, Jan 9, 2013 at 11:39 AM, Eduardo <span dir="ltr"><<a href="mailto:erocha.ssa@gmail.com" target="_blank">erocha.ssa@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">So, does the v block (the logically two-dimensional input array of<br>
values) still have memory positions for the lower-triangular? I mean<br>
do I still have to allocate a full v block even if the<br>
lower-triangular is never touched?<br></blockquote><div><br></div><div style>Yes.</div><div style><br></div><div style>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks a lot,<br>
Eduardo<br>
<br>
On Wed, Jan 9, 2013 at 3:24 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
> When using SBAIJ, you only set the upper-triangular part. If you want to be<br>
> able to set all entries anyway, run with -mat_ignore_lower_triangular or<br>
> call MatSetOption(mat,MAT_IGNORE_LOWER_TRIANGULAR,PETSC_TRUE). There is no<br>
> automatic way to symmetrize, but you're welcome to compute half and<br>
> symmetrize before calling MatSetValues().<br>
><br>
><br>
> On Wed, Jan 9, 2013 at 11:18 AM, Eduardo <<a href="mailto:erocha.ssa@gmail.com">erocha.ssa@gmail.com</a>> wrote:<br>
>><br>
>> Hi all,<br>
>><br>
>> Is there any way to assemble a block that is symmetric to a matrix<br>
>> (also symmetric)? I mean, as far as I know, the MatSetValues assumes a<br>
>> full block, i.e the parameter v in:<br>
>><br>
>>   PetscErrorCode  MatSetValues(Mat mat,PetscInt m,const PetscInt<br>
>> idxm[],PetscInt n,const PetscInt idxn[],const PetscScalar<br>
>> v[],InsertMode addv)<br>
>><br>
>> is a full block (local matrix) that is assembled into the global matrix<br>
>> mat.<br>
>><br>
>> Thanks in advance,<br>
>> Eduardo<br>
><br>
><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <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
</div></div>