Go to the documentation of this file.
47 scalar magDir =
mag(direction_);
52 <<
"(almost)zero vector as direction for cell " << cellNo_
55 else if (
mag(magDir - 1) > SMALL)
58 direction_ /=
mag(direction_);
66 cellNo_(readLabel(is)),
69 scalar magDir =
mag(direction_);
74 <<
"(almost)zero vector as direction for cell " << cellNo_
77 else if (
mag(magDir - 1) > SMALL)
80 direction_ /=
mag(direction_);