introduce implementation of normSquared to avoid often unnecessary roots and squaring;  compute norm by sqrt(normSquared); double->field_type; harmonized argument names with base class (i.e. v->f)

[[Imported from SVN: r12076]]