Dear dune-pdelab,<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>
Mainly I am using non overlapping grid (UG grid) and diagonal elements 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 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 clear="all"><br><br>Is it right? or  didn't I understand correctly?<br><br>Thanks & Regards<br>Nagaiah<br>