On Thu, Feb 27, 2014 at 03:01:32PM +0100, Andreas Lauser wrote: > On Thursday 27 February 2014 14:47:50 Markus Blatt wrote: > > On Thu, Feb 27, 2014 at 07:12:09AM +0000, Atgeirr Rasmussen wrote: > > 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. > > can't the length of the returned vector just be undefined for faceNormal()? > the vector can easily be normalized thereafter...
Actually, the use in opm-core should mandate what is returned here. This is why my interface was introduced in the first place. There the normal scaled with the face area is needed. If it is undefined one has to do additional computations that are actually not needed. > > > One could adapt Cpgrid to store the scaled face normals. Or is there > > and argument against this? > > By virtue of Gauss' integral theorem, I encounter more unit normals than > scaled ones in equations. I'd prefer unit normals in CpGrid for this > reason... What to st -- 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
