Non-uniform 2D mesh questions

yaron at oak-research.com yaron at oak-research.com
Tue Jan 30 17:07:18 CST 2007


Barry-
So far I only thought of having a single large sparse matrix.

Yaron


	-------Original Message-------
  From: Barry Smith
  Subject: Re: Non-uniform 2D mesh questions
  Sent: 30 Jan '07 10:58


   Yaron,

   Do you want to end up generating a single large sparse matrix? Like a
MPIAIJ
  matrix? Or do you want to somehow not store the entire huge matrix but
still
  be able to solve with the composed matrix? Or both?

   Barry


  On Mon, 29 Jan 2007, [LINK:
http://webmail.oak-research.com/compose.php?to=yaron@oak-research.com]
yaron at oak-research.com wrote:

  > Barry-
  > Yes, each block is a rectangular portion of the domain. Not so small
  > though (more like 100 x 100 nodes)
  >
  > Yaron
  >
  >
  > 	 -------Original Message-------
  > From: Barry Smith
  > Subject: Re: Non-uniform 2D mesh questions
  > Sent: 29 Jan '07 19:40
  >
  >
  > Yaron,
  >
  > Is each one of these "blocks" a small rectangular part of the
  > domain (like a 4 by 5 set of nodes)? I don't understand what you
  > want to do.
  >
  > Barry
  >
  >
  > On Mon, 29 Jan 2007, [LINK:
  > [LINK:
http://webmail.oak-research.com/compose.php?to=yaron@oak-research.com]
http://webmail.oak-research.com/compose.php?to=yaron@oak-research.com]
  > [LINK:
http://webmail.oak-research.com/compose.php?to=yaron@oak-research.com]
yaron at oak-research.com wrote:
  >
  > > Hi all
  > > I have a laplace-type problem that's physically built from repeating
  > > instances of the same block.
  > > I'm creaing matrices for the individual blocks, and I'd like to
reuse
  > > the individual block matrices in order to compose the complete
  > problem.
  > > (i.e if there 10K instances of 20 blocks, I'd like to build 20
  > matrices,
  > > then use them to compose the large complete matrix)
  > > Is a 2D DA the right object to do that? And if so, where can I find
a
  > > small example of building the DA object in parallel, then using the
  > > different (for every instance) mappings of local nodes to global
nodes
  > in
  > > order to build the complete matrix?
  > >
  > >
  > > Thanks
  > > Yaron
  > >
  >
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20070130/b3d1ef28/attachment-0001.htm>


More information about the petsc-users mailing list