On Thu, 24 Jun 2010 12:16:53 -0500 (CDT), Roy Stogner 
<royst...@ices.utexas.edu> wrote:
> I suspect there are only two reasons why we didn't support mismatched
> types:
> 
> when that assertion was written we only had serial and parallel
> vectors, a more extreme mismatch

Yes, these are different in a deep way, they have different
communicators and likely different local sizes.

> As long as we don't copy the type over with operator=, this shouldn't
> be a concern for PetscVector.

No, just recognize that

  LocalFunction(xghost,fghost);
  yplain = xghost;
  LocalFunction(yplain,fghost);

would fail in LocalFunction instead of operator=.  I think this fine,
but it is a slightly later error message.

> I'd say that we should at least enable 
> "parallelvec = ghostedvec", if everyone agrees.

And ghostedvec=parallelvec.

> "parallelvec = serialvec" and "ghostedvec = serialvec"

Would the semantics of the latter be

  VecGetLocalForm(g,&lf);
  VecCopy(s,lf);
  VecRestoreLocalForm(g,&lf);

or something else?  Would it include a VecGhostUpdateBegin/End?  If so,
would that use ADD_VALUES or INSERT_VALUES?

Jed

------------------------------------------------------------------------------
ThinkGeek and WIRED's GeekDad team up for the Ultimate 
GeekDad Father's Day Giveaway. ONE MASSIVE PRIZE to the 
lucky parental unit.  See the prize list and enter to win: 
http://p.sf.net/sfu/thinkgeek-promo
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to