Use different mesh for disc without interaction domain

Signed-off-by: Christian Glusa <[email protected]>

Author: cgcgcg

base/PyNucleus_base/DenseLinearOperator_{SCALAR}.pxi

@@ -196,6 +196,30 @@ cdef class {SCALAR_label}Dense_VectorLinearOperator({SCALAR_label}VectorLinearOp
                     y[i, k] += self.data[i, j, k]*x[j]
         return 0
 
+    cdef INDEX_t matvecTrans(self,
+                             {SCALAR}_t[::1] x,
+                             {SCALAR}_t[:, ::1] y) except -1:
+        cdef:
+            INDEX_t i, j, k
+        for i in range(self.num_rows):
+            for k in range(self.vectorSize):
+                y[i, k] = 0.
+            for j in range(self.num_columns):
+                for k in range(self.vectorSize):
+                    y[i, k] += self.data[j, i, k]*x[j]
+        return 0
+
+    cdef INDEX_t matvecTrans_no_overwrite(self,
+                                          {SCALAR}_t[::1] x,
+                                          {SCALAR}_t[:, ::1] y) except -1:
+        cdef:
+            INDEX_t i, j, k
+        for i in range(self.num_rows):
+            for j in range(self.num_columns):
+                for k in range(self.vectorSize):
+                    y[i, k] += self.data[j, i, k]*x[j]
+        return 0
+
     def isSparse(self):
         return False
 

base/PyNucleus_base/utilsFem.py

@@ -1073,10 +1073,8 @@ def process(self, override={}):
         if params['test']:
             import psutil
             p = psutil.Process()
-            try:
+            if hasattr(p, 'cpu_affinity'):
                 p.cpu_affinity(list(range(psutil.cpu_count())))
-            except AttributeError:
-                pass
         self._timer = TimerManager(self.logger, comm=self.comm, memoryProfiling=params['showMemory'])
 
         for fun in self.processHook:
@@ -1450,12 +1448,12 @@ def __init__(self, name, params=[]):
             elif p == int:
                 fields.append('[+-]?[0-9]+')
             elif p == float:
-                fields.append('[+-]?[0-9]+\.[0-9]*')
+                fields.append(r'[+-]?[0-9]+\.[0-9]*')
             elif p == bool:
                 fields.append('True|False')
             else:
                 raise NotImplementedError()
-        self.regexp = re.compile(name+'\(?'+','.join(['\s*(' + f + ')\s*' for f in fields])+'\)?')
+        self.regexp = re.compile(name+r'\(?'+','.join([r'\s*(' + f + r')\s*' for f in fields])+r'\)?')
 
     def match(self, s):
         return self.regexp.match(s) is not None
@@ -1516,6 +1514,8 @@ def __call__(self):
                             needToBuild = True
                         else:
                             dependencyLogger.log(self.logLevel, 'Values for {} are identical: \'{}\' == \'{}\''.format(prop, oldValue, newValue))
+                    elif newValue is None and oldValue is None:
+                        dependencyLogger.log(self.logLevel, 'Values for {} are identical: \'{}\' == \'{}\''.format(prop, oldValue, newValue))
                     elif newValue != oldValue:
                         cached_args[prop] = newValue
                         dependencyLogger.log(self.logLevel, 'Values for {} differ: \'{}\' != \'{}\', calling \'{}\''.format(prop, oldValue,

drivers/testDistOp.py

@@ -69,13 +69,13 @@
                                       h=h)
     elif nPP.domain == 'disc':
         if d.horizonToMeshSize <= 0. or nPP.kernel.horizon.value == np.inf:
-            h = 0.04/2**(nPP.noRef-3)
+            h = 0.04/2**(nPP.noRef-4)
         else:
             h = nPP.kernel.horizon.value/d.horizonToMeshSize/np.sqrt(2)
         mesh, _ = nonlocalMeshFactory(nPP.domain,
                                       kernel=nPP.kernel,
                                       boundaryCondition=HOMOGENEOUS_DIRICHLET,
-                                      h=h,
+                                      hTarget=h,
                                       max_volume=h**2/2,
                                       projectNodeToOrigin=False)
     elif d.domain == 'gradedDisc':

fem/PyNucleus_fem/factories.py

@@ -145,7 +145,7 @@ def __call__(self, mesh, *args, **kwargs):
 
 # meshes
 from . mesh import (simpleInterval, simpleSquare, simpleLshape, simpleBox, box,
-                    circle, graded_circle, cutoutCircle, twinDisc, dumbbell, wrench,
+                    uniform_disc, circle, graded_circle, cutoutCircle, twinDisc, dumbbell, wrench,
                     Hshape, ball, rectangle, crossSquare,
                     gradedSquare, gradedBox,
                     squareWithCircularCutout, boxWithBallCutout,
@@ -178,6 +178,7 @@ def __call__(self, mesh, *args, **kwargs):
 meshFactory.register('squareWithInteraction', squareWithInteractions, 2)
 meshFactory.register('simpleLshape', simpleLshape, 2, aliases=['Lshape', 'L-shape'])
 meshFactory.register('circle', circle, 2, aliases=['disc', 'unitDisc', 'ball2d', '2dball'])
+meshFactory.register('uniform_disc', uniform_disc, 2, aliases=['uniform_ball2d', '2dball_uniform'])
 meshFactory.register('graded_circle', graded_circle, 2, aliases=['gradedCircle'])
 meshFactory.register('discWithInteraction', discWithInteraction, 2)
 meshFactory.register('twinDisc', twinDisc, 2)

fem/PyNucleus_fem/mesh.py

@@ -720,10 +720,7 @@ def discWithInteraction(radius, horizon, h=0.25, max_volume=None, projectNodeToO
                                      innerRadius=radius,
                                      max_volume=max_volume)
     else:
-        return circle(n,
-                      radius=radius,
-                      max_volume=max_volume,
-                      projectNodeToOrigin=projectNodeToOrigin)
+        return uniform_disc(radius=radius)
 
 
 def gradedDiscWithInteraction(radius, horizon, mu=2., h=0.25, max_volume=None, projectNodeToOrigin=True):
@@ -946,6 +943,23 @@ def Lshape(d):
     return mesh2d(vertices, cells)
 
 
+def uniform_disc(radius=1., **kwargs):
+    points = [(0., 0.)]
+    cells = []
+    n = 6
+    for i in range(n):
+        points.append((radius*np.cos(i*2*np.pi/n), radius*np.sin(i*2*np.pi/n)))
+    for i in range(1, len(points)-1):
+        cells.append((0, i, i+1))
+    cells.append((0, len(points)-1, 1))
+
+    mesh = mesh2d(np.array(points, dtype=REAL),
+                  np.array(cells, dtype=INDEX))
+    from . meshCy import radialMeshTransformer
+    mesh.setMeshTransformation(radialMeshTransformer())
+    return mesh
+
+
 def circle(n, radius=1., returnFacets=False, projectNodeToOrigin=True, **kwargs):
     from meshpy.triangle import MeshInfo, build
     mesh_info = MeshInfo()

nl/PyNucleus_nl/clusterMethodCy.pyx

@@ -1193,7 +1193,7 @@ cdef class tree_node:
                 for i in range(self.coefficientsDownVec.shape[0]):
                     for l in range(self.coefficientsDownVec.shape[1]):
                         y[dof, l] += self.value[componentNo, k, i]*self.coefficientsDownVec[i, l]
-            k += 1
+                k += 1
         else:
             for c in self.children:
                 for i in range(c.transferOperator.shape[0]):

tests/cache_runFractional.py--domaindisc--sconst(0.25)--problemconstant--elementP0--solvercg-mg--matrixFormatH2

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.13400055486309811
-  L2 error: 0.0530257052232404
-  L2 error interpolated: 0.017212921206495174
-  Linf error interpolated: 0.04041162523460817
-  relative Hs error: 0.09110983276196648
-  relative L2 error: 0.04257104506006856
-  relative interpolated L2 error: 0.013812422648097445
-  relative interpolated Linf error: 0.0469573643243436
+  Hs error: 0.13190712640577038
+  L2 error: 0.05014710987068258
+  L2 error interpolated: 0.01711413239312614
+  Linf error interpolated: 0.042577264837976025
+  relative Hs error: 0.08968646614351379
+  relative L2 error: 0.04025999965392971
+  relative interpolated L2 error: 0.013732997326192832
+  relative interpolated Linf error: 0.04947323085350075
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconst(0.25)--problemconstant--elementP0--solvercg-mg--matrixFormatdense

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.1413324571924916
-  L2 error: 0.05318118990883034
-  L2 error interpolated: 0.017767471098211238
-  Linf error interpolated: 0.0402337391432151
-  relative Hs error: 0.0960949493963012
-  relative L2 error: 0.042695874056279615
-  relative interpolated L2 error: 0.014257418438872843
-  relative interpolated Linf error: 0.046750664842368186
+  Hs error: 0.1403179566911808
+  L2 error: 0.050326820499312745
+  L2 error interpolated: 0.01771990305453638
+  Linf error interpolated: 0.042613546494036036
+  relative Hs error: 0.09540516888676681
+  relative L2 error: 0.0404042781550261
+  relative interpolated L2 error: 0.014219089561682037
+  relative interpolated Linf error: 0.04951538881627348
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconst(0.25)--problemconstant--elementP1--solvercg-mg--matrixFormatH2

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.180573736324132
-  L2 error: 0.05758412182442717
-  L2 error interpolated: 0.017270441090162007
-  Linf error interpolated: 0.08537962523475567
-  relative Hs error: 0.12277593129747377
-  relative L2 error: 0.04623071460552205
-  relative interpolated L2 error: 0.013926906465042695
-  relative interpolated Linf error: 0.09919993999254309
+  Hs error: 0.18185981625380002
+  L2 error: 0.05782159835474865
+  L2 error interpolated: 0.017338240822480377
+  Linf error interpolated: 0.08846749274971044
+  relative Hs error: 0.12365036444762217
+  relative L2 error: 0.04642136976098792
+  relative interpolated L2 error: 0.013982271432707184
+  relative interpolated Linf error: 0.1027876375415333
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconst(0.25)--problemconstant--elementP1--solvercg-mg--matrixFormatdense

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.18352538491659642
-  L2 error: 0.05767242742488538
-  L2 error interpolated: 0.01741375173232188
-  Linf error interpolated: 0.08492371703495766
-  relative Hs error: 0.12478282007421258
-  relative L2 error: 0.04630160968707422
-  relative interpolated L2 error: 0.014042472355826032
-  relative interpolated Linf error: 0.09867023438728048
+  Hs error: 0.1839933908571473
+  L2 error: 0.057885119791182965
+  L2 error interpolated: 0.017437708550776834
+  Linf error interpolated: 0.08802653240222696
+  relative Hs error: 0.12510102728625555
+  relative L2 error: 0.046472367176700564
+  relative interpolated L2 error: 0.014062486305142948
+  relative interpolated Linf error: 0.10227529938252666
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconst(0.75)--problemconstant--elementP1--solvercg-mg--matrixFormatH2

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.05860698560008363
-  L2 error: 0.0021486188473719033
-  L2 error interpolated: 0.001290467850492567
-  Linf error interpolated: 0.0014304080053778728
-  relative Hs error: 0.06761001433404402
-  relative L2 error: 0.004579199453127269
-  relative interpolated L2 error: 0.0027517017386181047
-  relative interpolated Linf error: 0.003417393926565291
+  Hs error: 0.059725648882225826
+  L2 error: 0.0022274080583107514
+  L2 error interpolated: 0.0013188143633481523
+  Linf error interpolated: 0.0017403068951004769
+  relative Hs error: 0.06890052330266251
+  relative L2 error: 0.004747117328410178
+  relative interpolated L2 error: 0.002812374690576744
+  relative interpolated Linf error: 0.004157774698768523
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconst(0.75)--problemconstant--elementP1--solvercg-mg--matrixFormatdense

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.05957553465651644
-  L2 error: 0.0021947415915123092
-  L2 error interpolated: 0.0013546361081100044
-  Linf error interpolated: 0.001421361318103019
-  relative Hs error: 0.0687273489814099
-  relative L2 error: 0.004677497597073466
-  relative interpolated L2 error: 0.002888529561165253
-  relative interpolated Linf error: 0.0033957804470319086
+  Hs error: 0.060319591944560894
+  L2 error: 0.002256341047519089
+  L2 error interpolated: 0.0013588886928846355
+  Linf error interpolated: 0.0017321548640112927
+  relative Hs error: 0.06958570611059765
+  relative L2 error: 0.004808780162896719
+  relative interpolated L2 error: 0.002897833291318784
+  relative interpolated Linf error: 0.004138298646181701
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconstantNonSym(0.25)--problemconstant--elementP1--solvergmres-jacobi--matrixFormatH2

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.18057373621467981
-  L2 error: 0.057584121294463754
-  L2 error interpolated: 0.017270390898283507
-  Linf error interpolated: 0.08537914587854994
-  relative Hs error: 0.12277593122305488
-  relative L2 error: 0.046230714180047344
-  relative interpolated L2 error: 0.013926865990245698
-  relative interpolated Linf error: 0.09919938304342657
+  Hs error: 0.18185981616987204
+  L2 error: 0.057821602128125614
+  L2 error interpolated: 0.01733821843880583
+  Linf error interpolated: 0.08846702811244911
+  relative Hs error: 0.12365036439055775
+  relative L2 error: 0.04642137279039791
+  relative interpolated L2 error: 0.013982253381590213
+  relative interpolated Linf error: 0.10278709769390203
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconstantNonSym(0.25)--problemconstant--elementP1--solvergmres-jacobi--matrixFormatdense

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.1835253848079427
-  L2 error: 0.05767242809363617
-  L2 error interpolated: 0.01741370463401979
-  Linf error interpolated: 0.08492323905357663
-  relative Hs error: 0.12478282000033661
-  relative L2 error: 0.04630161022397269
-  relative interpolated L2 error: 0.014042434375693132
-  relative interpolated Linf error: 0.09866967903553016
+  Hs error: 0.18399339204392906
+  L2 error: 0.05788512423832981
+  L2 error interpolated: 0.017437686447626173
+  Linf error interpolated: 0.08802604582549106
+  relative Hs error: 0.12510102809317378
+  relative L2 error: 0.046472370747038534
+  relative interpolated L2 error: 0.01406246848025199
+  relative interpolated Linf error: 0.10227473404409995
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconstantNonSym(0.75)--problemconstant--elementP1--solvergmres-jacobi--matrixFormatH2

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.058607115732549314
-  L2 error: 0.0021486244053874164
-  L2 error interpolated: 0.0012904756531338202
-  Linf error interpolated: 0.0014305222217047692
-  relative Hs error: 0.06761016445706741
-  relative L2 error: 0.004579211298532808
-  relative interpolated L2 error: 0.002751718376414617
-  relative interpolated Linf error: 0.0034176668012838207
+  Hs error: 0.0597255551387594
+  L2 error: 0.00222740318228125
+  L2 error interpolated: 0.001318806267789667
+  Linf error interpolated: 0.0017403189548499884
+  relative Hs error: 0.06890041515860693
+  relative L2 error: 0.0047471069364732115
+  relative interpolated L2 error: 0.002812357426779486
+  relative interpolated Linf error: 0.004157803510768023
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--sconstantNonSym(0.75)--problemconstant--elementP1--solvergmres-jacobi--matrixFormatdense

@@ -1,13 +1,13 @@
 Timers: {}
 errors:
-  Hs error: 0.059575659109656745
-  L2 error: 0.002194748254364462
-  L2 error interpolated: 0.0013546456020427384
-  Linf error interpolated: 0.0014213792199474634
-  relative Hs error: 0.06872749255266705
-  relative L2 error: 0.00467751179713923
-  relative interpolated L2 error: 0.0028885498053512672
-  relative interpolated Linf error: 0.003395823216405576
+  Hs error: 0.06031949823523683
+  L2 error: 0.0022563364021103744
+  L2 error interpolated: 0.0013588811288227196
+  Linf error interpolated: 0.0017321724455461873
+  relative Hs error: 0.06958559800592927
+  relative L2 error: 0.0048087702624655315
+  relative interpolated L2 error: 0.0028978171609391934
+  relative interpolated Linf error: 0.004138340650302444
 meshes: {}
 results: {}
 sysInfo: {}

tests/cache_runFractional.py--domaindisc--stwoDomainNonSym(0.25,0.75)--problemknownSolution--elementP1--solvergmres-mg--matrixFormatH2

@@ -1,11 +1,11 @@
 Timers: {}
 errors:
-  L2 error: 0.005930197043816925
-  L2 error interpolated: 0.003413020217808587
-  Linf error interpolated: 0.01082136808655848
-  relative L2 error: 0.0051832219786591745
-  relative interpolated L2 error: 0.002984844935336751
-  relative interpolated Linf error: 0.01082136808655848
+  L2 error: 0.005826340789746348
+  L2 error interpolated: 0.0031194213241089513
+  Linf error interpolated: 0.011279648619087423
+  relative L2 error: 0.0050924475887455896
+  relative interpolated L2 error: 0.0027282922216046363
+  relative interpolated Linf error: 0.011279648619087423
 meshes: {}
 results: {}
 sysInfo: {}