On Thu, Feb 27, 2014 at 07:12:09AM +0000, Atgeirr Rasmussen wrote: > Your observation is correct. The normals in the UnstructuredGrid are required > to have length equal to the > area of the corresponding face (see the Doxygen doc for UnstructuredGrid). > It was just a surprise to see such a difference and took time to realize it. I am probably not the last one to stumble over this. As currently the only place where this caused trouble was opm/core/pressure/tpfa/..., I worked around the problem there. The problem is that there are now faceNormal(const UG&, int) and faceNormal(const CpGrid, int) that actually return different things. This should probably be fixed.
One could adapt Cpgrid to store the scaled face normals. Or is there and argument against this? If this not an option, then it is a better idea to return a copy of a FieldVector from faceNormal(const CpGrid, int) that returns CpGrid::faceNormal(int) scaled with the face area. Markus -- Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany, USt-Id: DE279960836 Tel.: +49 (0) 160 97590858
signature.asc
Description: Digital signature
_______________________________________________ Opm mailing list [email protected] http://www.opm-project.org/mailman/listinfo/opm
