In issue #45 (not submitted by me) it was indicated that this should be
directed here: https://bitbucket.org/fenics-project/dolfin/issue/45/

P1_DG FunctionSpace initialisation seems to exhibit quadratic serial
scaling. e.g. with the following test code:

from dolfin import *
m = 128
for p in range(1, 15):
  n = 2 ** p
  mesh = UnitSquareMesh(m, n)
  timer = Timer("dg")
  space = FunctionSpace(mesh, "DG", 1)
  print "%ix%i: %.6f" % (m, n, timer.stop())

I get (FEniCS 1.2):

128x2: 0.017284
128x4: 0.007458
128x8: 0.022625
128x16: 0.077661
128x32: 0.295081
128x64: 1.119631
128x128: 4.450091
128x256: 17.725417
128x512: 70.331817

P2_DG exhibits a similar issue, whereas P0 seems fine. I haven't tested
higher degree.

James Maddison

_______________________________________________
fenics-support mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics-support

Reply via email to