3
3
from builtins import object
4
4
__docformat__ = 'restructuredtext'
5
5
6
+ import logging
6
7
import os
7
8
8
9
from fipy import input
@@ -41,6 +42,9 @@ def __init__(self, coeff=1., var=None):
41
42
self ._RHSvector = None
42
43
self .var = var
43
44
45
+ self ._log = logging .getLogger (self .__class__ .__module__
46
+ + "." + self .__class__ .__name__ )
47
+
44
48
def _calcVars (self ):
45
49
raise NotImplementedError
46
50
@@ -105,6 +109,9 @@ def _getMatrixClass(self, solver, var):
105
109
return SparseMatrix
106
110
107
111
def _prepareLinearSystem (self , var , solver , boundaryConditions , dt ):
112
+
113
+ self ._log .debug ("BEGIN _prepareLinearSystem" )
114
+
108
115
solver = self .getDefaultSolver (var , solver )
109
116
110
117
var = self ._verifyVar (var )
@@ -151,6 +158,8 @@ def _prepareLinearSystem(self, var, solver, boundaryConditions, dt):
151
158
from fipy import input
152
159
input ()
153
160
161
+ self ._log .debug ("END _prepareLinearSystem" )
162
+
154
163
return solver
155
164
156
165
def solve (self , var = None , solver = None , boundaryConditions = (), dt = None ):
@@ -173,10 +182,14 @@ def solve(self, var=None, solver=None, boundaryConditions=(), dt=None):
173
182
Timestep size.
174
183
"""
175
184
185
+ self ._log .debug ("BEGIN solve" )
186
+
176
187
solver = self ._prepareLinearSystem (var , solver , boundaryConditions , dt )
177
188
178
189
solver ._solve ()
179
190
191
+ self ._log .debug ("END solve" )
192
+
180
193
def sweep (self , var = None , solver = None , boundaryConditions = (), dt = None , underRelaxation = None , residualFn = None , cacheResidual = False , cacheError = False ):
181
194
r"""
182
195
Builds and solves the `Term`'s linear system once. This method
0 commit comments