[petsc-users] MatCreateMPIAIJWithSplitArrays

[Matthew Knepley]

On Tue, Aug 2, 2011 at 5:09 PM, Kuhlemann, Verena <vkuhlem at emory.edu> wrote:

>  HI,
>
>  I am creating a matrix with MatCreateMPIAIJWithSplitArrays, but I would
> like the arrays that I am using
> to create the matrix to be freed once the matrix is destroyed. Is this the
> correct way to do this?
>
>
> MatCreateMPIAIJWithSplitArrays(PETSC_COMM_WORLD,ne,nf,PETSC_DETERMINE,PETSC_DETERMINE,ia_E_d,ja_E_d,v_E_d,ia_E_o,ja_E_o,v_E_o,E);
> //set flags so that arrays are freed once the matrix is destroyed
>  mpimat= (Mat_MPIAIJ*)(*E)->data;
>  a=(Mat_SeqAIJ*)(mpimat->A)->data;
>  b=(Mat_SeqAIJ*)(mpimat->B)->data;
>  a->free_a  = PETSC_TRUE;
>  a->free_ij = PETSC_TRUE;
>  b->free_a  = PETSC_TRUE;
>  b->free_ij = PETSC_TRUE;
>

Yes, this should work.

   Matt


> Thanks,
>
>  Verena
>
> ------------------------------
>
> This e-mail message (including any attachments) is for the sole use of
> the intended recipient(s) and may contain confidential and privileged
> information. If the reader of this message is not the intended
> recipient, you are hereby notified that any dissemination, distribution
> or copying of this message (including any attachments) is strictly
> prohibited.
>
> If you have received this message in error, please contact
> the sender by reply e-mail message and destroy all copies of the
> original message (including attachments).
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110802/17a7c475/attachment.htm>