Dear Felix Heimann,<br>Thank you for your response. I am sorry to say that I did not give complete details.<br>you are right. I am referring vector valued solution which can be represented by PowerGridFunctionSpace in PDELab.<br>
I know only that we can set physical boundary conditions (Neumann or Dirichlet BC's)<br>using Dune::PDELab::constriants() function. I do not know that we can impose <br>constraints (artificial BC) for ghost elements as well. Also, I am not sure how to set the constraints <br>
for ghost elements in this constraints() function. <br><br>Is it a good idea to set artificial boundary conditions for ghost elements <br>in non overlapping grids case using Dune::PDELab::constraints(..)?<br>I would expect that these atrificial boundary conditions should be automatically taken care by<br>
gridoperator function routines. am I stupid here or did not understand parallel pdelab routines properly?<br><br>Thank you in advance.<br><br>Best Regards<br>Nagaiah<br><br><br><br><div class="gmail_quote">On Fri, Jan 14, 2011 at 4:46 PM, Felix Heimann <span dir="ltr"><<a href="mailto:felix.heimann@iwr.uni-heidelberg.de">felix.heimann@iwr.uni-heidelberg.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear Nagaiah,<br>
<br>
I am not sure, what exactly you mean by "more than one solution<br>
component". If you refer to e.g. a vector valued solution, then this<br>
should be represented by e.g. a PowerGridFunctionSpace. In such a case,<br>
the set_trivial_row() function will work correctly as long as the<br>
indices in the constraints container are consistent to the indices of<br>
the grid function space (this should be guaranteed if you set up the<br>
container with the Dune::PDELab::constraints(..) function).<br>
<br>
However, this functionality is the same for all trivially constrained<br>
dofs, independent of whether they result from a ghost cell in a domain<br>
decomposed non-overlapping application or a Dirichlet boundary in a<br>
standard application. Hence, I do not see the connection to the specific<br>
grid i.e. application and may not have grasped your problem. May be you<br>
can clarify...<br>
<br>
Best Regards,<br>
Felix Heimann<br>
<br>
<br>
<br>
Am Mittwoch, den 12.01.2011, 22:57 +0100 schrieb Nagaiah Chamakuri:<br>
<div><div></div><div class="h5">> Dear dune-pdelab,<br>
><br>
> I have a question about "set_trivial_row" function in<br>
> new_gridoperatorspace.hh at line 885  (precisely I am using Jacobian<br>
> function routine to construct matrices) and defined in<br>
> gridoperatorspaceutilities.hh at line 971.<br>
> Mainly I am using non overlapping grid (UG grid) and diagonal elements<br>
> of the ghost nodes are set to zero in the global matrix<br>
> in case of more than one solution component presents.<br>
> This function routine works well in case of one solution component<br>
> presents in the problem.<br>
><br>
><br>
> right answer should look like as follows (up to my knowledge):<br>
><br>
> for j=1 to no.of.comp<br>
> B::clear_row(i+j,globalcontainer)<br>
><br>
> B::access(globalcontainer,i+j,i+j) = 1,<br>
><br>
><br>
><br>
> Is it right? or  didn't I understand correctly?<br>
><br>
> Thanks & Regards<br>
> Nagaiah<br>
</div></div>> _______________________________________________<br>
> dune-pdelab mailing list<br>
> <a href="mailto:dune-pdelab@dune-project.org">dune-pdelab@dune-project.org</a><br>
> <a href="http://lists.dune-project.org/mailman/listinfo/dune-pdelab" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune-pdelab</a><br>
<br>
--<br>
Felix Heimann<br>
Universität Heidelberg<br>
Interdisziplinäres Zentrum für Wissenschaftliches Rechnen<br>
Arbeitsgruppe Paralleles Rechnen<br>
IWR 368, Raum 422<br>
Tel: 06221 / 54 8881<br>
<br>
<br>
_______________________________________________<br>
dune-pdelab mailing list<br>
<a href="mailto:dune-pdelab@dune-project.org">dune-pdelab@dune-project.org</a><br>
<a href="http://lists.dune-project.org/mailman/listinfo/dune-pdelab" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune-pdelab</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Dr. Chamakuri Nagaiah<br>Institute of Mathematics and Scientific Computing<br>Heinrichstrasse-36, University of Graz,<br>A-8010, Graz, Austria. Ph(office) : (0043) 316 380 5063<br>
<br>