[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Poisson Eqn in a concentric cube
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Poisson Eqn in a concentric cube
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Tue, 24 Jul 2007 04:05:10 -0500
- Dkim-signature: a=rsa-sha1; c=relaxed/relaxed; d=gmail.com; s=beta; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=uvItlqdDKbqAItAZHr5F433BFhJha9FSfua4Dn7KnxDdEqt+ZrO9wZgrWQ6mexCo5gmMihcL0lDoD6x9lw6nrjrbPyc1D2lNRI3NoYFuey5KpC2QvA9eCeHc9800G4Z2XHnRpHcRTx4Vb+rBqivKmQsHHkURWMUkPVZpTRFeKRw=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=fPjh/xRpKjpIKmka6+jynSM8O7j0KDotdMjcOALwC4WynZO39Y/fVugEBZxZZk914MYjJS7HkGD0TaUZLpWguIKKImp9+aEF51juv9bB0FhQKnH41uw+tfdWO7SfmGEatmAIRP9goLLnTXEIF+iP0D41wLzO+qJx36K1w4hWg9k=
- In-reply-to: <A970687040E72044810535C74B5F0E4901071FB5@MAIL22.student.main.ntu.edu.sg>
- References: <A970687040E72044810535C74B5F0E4901071FB5@MAIL22.student.main.ntu.edu.sg>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
Currently, the DA does not directly support cutout grids, as you have
shown above. However, you could just substitute the identity for all
points in the cutout, and only suffer from an excessive preallocation.
Does this make sense?
Thanks,
Matt
On 7/23/07, #DOMINIC DENVER JOHN CHANDAR# <DOMI0002@xxxxxxxxxx> wrote:
Hi,
I would like to solve the Poisson equation (3D) in a concentric type of
cube as in the attached figure using PETSc. The solution domain would be the
shaded region only and it extends inside the cube. I would like to have
suggestions on ways to form and assemble the matrix associated with the
discretization. Also, is it a good practise to discretize the entire domain
and blank the grid points which are within the cube not in the domain ?
Regards,
Dominic
--
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