[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: question on MatMatMult
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: question on MatMatMult
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Mon, 17 Mar 2008 17:40:12 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=oD1/KTf5oil/VASKLtFOT+OVmsNUPH+0Xj8FAFrqj9g=; b=cbztk84LyzKl+m7+ZlBjO8cimozV8J23BO2FMVYgEP1+cKfMFxtg1v71tDOS5bH0L9XnMYtVhUIDWfmSHnp5x2jdLFDD03K3+gq3JlrhQkgxRPrqsyq72KTiFOsgw2mE2IPE/a6MuZNUWuas8kMJo8nG2zsSpZ947Dy29Uj34iU=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=YCNHbP3WPnnKhiGVVpHlI06lWOqDASD3zRZ8SS2EV1TET1Z4efU6Xl5GfdEjc5liRrstClTbLRAnqNbjaE7DRiJ/QuyVctVsB+fs7MPCDSQTSPontCe8LMNnGnYzoJJsYzftpkiOYEQV9/Bk8dDSqBR2i7eI7VD+jqhzyOBgnAA=
- In-reply-to: <47DEF093.5040605@gmail.com>
- References: <47DEF093.5040605@gmail.com>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
On Mon, Mar 17, 2008 at 5:28 PM, Randall Mackie <rlmackie862@xxxxxxxxx> wrote:
> If I want to compute H = L^T L where L is a sparse matrix which is an approximation to the
> laplacian, and hence H is the biharmonic (also sparse), and if I have L as an MPI matrix
> in PETSc, will MatMatMult work for this (assuming I create the transpose of L first).
>
> In other words, does MatMatMult look at the non-zero structure only that would result,
> or does it think the result is a dense matrix?
It builds the structure dynamically, which explains the "fill"
argument. This is not a great
thing to do unless you have no idea how to form it directly.
Matt
> Thanks, Randy
--
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