> double cuentoquads2=0;
> for (unsigned int cell=cuentohexs; cell<(cuentoquads+cuentohexs);
> ++cell) {
> if ((int)cells[cell].material_id ==boundary_value)
> cuentoquads2++;
> }
> std::cout <<"NUMBER OF QUADS: "<<cuentoquads2 <<std::endl; //where
> cuentoquads2= 0, i.e., this is different from cuentoquads, This is
> *WRONG*.
I believe that 'cells' is an array that has only 'cuentohexs' elements. You
access invalid elements here. If you want to investigate faces (quads
here) you have to look into the 'subcelldata' array, specifically the
boundary_quads member. Take a look at the queue of if-statements
immediately above the GridTools::delete_unused_vertices call in
grid_in.cc.
Best
W.
--
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
