Changeset 3266

Timestamp:

10/15/2009 10:17:06 AM (4 months ago)

Author:

guyer

Message:

Use existing L2norm() function

Files:

• trunk/fipy/terms/term.py (modified) (2 diffs)

trunk/fipy/terms/term.py

@@ -81 +81 @@
     def _calcResidual(self, var, matrix, RHSvector):
 
-        L2norm = numerix.sqrt(numerix.sum(self._calcResidualVector(var, matrix, RHSvector)**2))
-        RHSL2norm = numerix.sqrt(numerix.sum(RHSvector**2))
+        L2norm = numerix.L2norm(self._calcResidualVector(var, matrix, RHSvector))
+        RHSL2norm = numerix.L2norm(RHSvector)
         
         if RHSL2norm == 0:
@@ -193 +193 @@
     def justResidualVector(self, var, solver=None, boundaryConditions=(), dt=1., underRelaxation=None, residualFn=None):
         r"""
-        Builds and the `Term`'s linear system once. This method
+        Builds the `Term`'s linear system once. This method
         also recalculates and returns the residual as well as applying
         under-relaxation.