[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
question on MatMatMult
- To: petsc-users@xxxxxxxxxxx
- Subject: question on MatMatMult
- From: Randall Mackie <rlmackie862@xxxxxxxxx>
- Date: Mon, 17 Mar 2008 15:28:35 -0700
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:user-agent:mime-version:to:subject:content-type:content-transfer-encoding; bh=HfvE849FcJA/kcE6nu1AWb++RXSXqkCKeC6LpitVpzA=; b=oFW+oMLjw5Bl77h5RYxID7ox6ow77WJa2AI0gA2W/tDE6asZzTFD8sP3eV2qym79aY7Mllyhp0lLxC7dzTgaG44WKf90v15oIRUkNIBab7jBsFa2PIdD68mIcS3ZD0Z0w1gP3fiZEUPeRIqyKkiw5W/2zEDp+wRodkGCMm41PM4=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:user-agent:mime-version:to:subject:content-type:content-transfer-encoding; b=Mi9G0YLdSelGeB3TQJGLMkyFKZa21hXVxVDlgIVzTzAhPOsbNzPCQuevZvmL1k99MCLXAyNJdZ1RalLjvc9QN+6JuJ0JqDSkI9mHx9oO/jUBNGzQir8smB72SEC7TXdQFw1V/CZmobYMzZAJ3dUDJKFkfmfwX4zKbN0MA4RQmlo=
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
- User-agent: Thunderbird 2.0.0.12 (Windows/20080213)
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?
Thanks, Randy