Fixed a bit notation.

src/bvpy/bvp.py

@@ -96,7 +96,7 @@ class BVP(object):
         self._dirichlet_class: list = []
         self.neumannCondition: dict = {'constant': [],
                                        'variable': []}
-        self.neumann_domain_id: list = {'constant': [],
+        self._neumann_domain_id: list = {'constant': [],
                                         'variable': []}
 
         self.solver: Union[LinearSolver, NonLinearSolver] = None
@@ -303,11 +303,11 @@ class BVP(object):
 
         if isinstance(dir_class, (ConstantNeumann)):
             self.neumannCondition['constant'].append(dir_class.apply(self.domain.bdata))
-            self.neumann_domain_id['constant'].append(dir_class._id)
+            self._neumann_domain_id['constant'].append(dir_class._id)
         elif isinstance(dir_class, (VariableNeumann, NormalNeumann)):
             self.neumannCondition['variable'].append(dir_class.apply(self.functionSpace,
                                         self.domain.bdata))
-            self.neumann_domain_id['variable'].append(dir_class._id)
+            self._neumann_domain_id['variable'].append(dir_class._id)
 
     def add_boundary_condition(self, bnd_cond,
                                method='topological',
@@ -402,11 +402,11 @@ class BVP(object):
         if isinstance(self.vform.rhs, Zero):
 
             for neumann, domain_id in zip(self.neumannCondition['constant'],
-                                          self.neumann_domain_id['constant']):
+                                          self._neumann_domain_id['constant']):
                 lhs -= fe.dot(neumann, #.apply(self.domain.bdata),
                               self.testFunction)*self.domain.ds(domain_id)
             for neumann, domain_id in zip(self.neumannCondition['variable'],
-                                          self.neumann_domain_id['variable']):
+                                          self._neumann_domain_id['variable']):
                 lhs -= fe.dot(neumann, #.apply(self.functionSpace,
                                         #    self.domain.bdata),
                               self.testFunction)*self.domain.ds(domain_id)
@@ -422,11 +422,11 @@ class BVP(object):
         else:
 
             for neumann, domain_id in zip(self.neumannCondition['constant'],
-                                          self.neumann_domain_id['constant']):
+                                          self._neumann_domain_id['constant']):
                 rhs += fe.dot(neumann, #.apply(self.domain.bdata),
                               self.testFunction)*self.domain.ds(domain_id)
             for neumann, domain_id in zip(self.neumannCondition['variable'],
-                                          self.neumann_domain_id['variable']):
+                                          self._neumann_domain_id['variable']):
                 rhs += fe.dot(neumann, #.apply(self.functionSpace,
                                         #    self.domain.bdata),
                               self.testFunction)*self.domain.ds(domain_id)