Hi dune-pdelab,<br><br>Yesterday I was trying to send a message to <a href="mailto:dune-pdelab@dune-project.org">dune-pdelab@dune-project.org</a> mailing list <br>but, it seems, due to lack of permissions my message was not delivered.<br>
So I am mailing to common dune mailing list. <br><br>I have a question about "set_trivial_row" function in new _gridoperatorspace.hh at line 885  (precisely I am using Jacobian<br>function routine to construct matrices) and defined in gridoperatorspaceutilities.hh at line 971. <br>
If one uses <b>non overlapping grid(</b>for instance UG grid) and more than one solution component then this function routine<br>adds a zero entries for some diagonal entries of the ghost nodes in the global matrix.<br>
This function routine works well in case of one solution component presents in the problem.<br><br><br>right answer should looks like as follows (up to my knowledge):<br><br>f<b>or j=1 to no.of.comp</b><br>{<br>    B::clear_row(i<b>+j</b>,globalcontainer)<br>
   
B::access(globalcontainer,i+j,i+j) = 1,<br>}<br clear="all"><br><br>Is it right? or  didn't I understand correctly?<br><br>Thank you in advance.<br><br>Regards<br>Nagaiah<br>