ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

I'm pretty sure the problem is a few changes we made to get BDM
elements working. I was looking at it but ran out of time. I'll fix it
as soon as I can (which should be early next week). If you want to
take a look, try to find out at which changeset it broke.

I tried to simplify the form a bit and find the simplest possible form
where it still breaks, but it seems it only breaks if the form is
complicated enough (has enough factors in it).

/Anders


On Thu, May 24, 2007 at 09:55:15PM +0200, Garth N. Wells wrote:
> Any ideas what going on here? It's causing problem for me.
> 
> Garth
> 
> Anders Logg wrote:
> > I'll take a look.
> > 
> > /Anders
> > 
> > 
> > On Thu, May 17, 2007 at 03:45:10PM +0200, Garth N. Wells wrote:
> >> When trying to compile a form file, I'm getting the error
> >>
> >> Computing tensor representation for term 0... *** dimensions too large.
> >> *** To get more information about this error, rerun ffc with the option -d1.
> >>
> >> Running with -d1,
> >>
> >>   File
> >> "/usr/lib/python2.5/site-packages/ffc/compiler/representation/tensor/monomialintegration.py", 
> >>
> >> line 54, in integrate
> >>      A0 = __compute_product(psis, vscaling * monomial.numeric * weights)
> >>    File
> >> "/usr/lib/python2.5/site-packages/ffc/compiler/representation/tensor/monomialintegration.py", 
> >>
> >> line 268, in __compute_product
> >>      A0 = numpy.transpose(A0, rearrangement)
> >>    File "/usr/lib/python2.5/site-packages/numpy/core/fromnumeric.py",
> >> line 118, in transpose
> >>      return transpose(axes)
> >> ValueError: dimensions too large.
> >>
> >> I was able to compile this form last weeks. Below is the form.
> >>
> >> Garth
> >>
> >>
> >>
> >> scalar1 = FiniteElement("Lagrange", "triangle", 1)
> >>
> >> r  = TestFunction(scalar1)
> >> p1 = TrialFunction(scalar1)
> >>
> >> p0   = Function(scalar1)
> >> lam  = Function(scalar1)
> >> FFF  = Function(scalar1)
> >>
> >> stab = Function(scalar1)
> >> h    = Constant("triangle")
> >>
> >> k     = Constant("triangle")
> >> theta = Constant("triangle")
> >> dt    = Constant("triangle")
> >>
> >> beta = mult(h,stab)
> >> SUPG = - beta*lam*k*dot(grad(p0),grad(r))
> >>
> >> L_saturation_stab = SUPG*lam*k*dot(grad(p0),grad(FFF))*dx
> >>
> >> L = L_saturation_stab
> >>
> >> _______________________________________________
> >> FFC-dev mailing list
> >> FFC-dev@xxxxxxxxxx
> >> http://www.fenics.org/mailman/listinfo/ffc-dev
> > _______________________________________________
> > FFC-dev mailing list
> > FFC-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/ffc-dev
> > 
> 
> 
> _______________________________________________
> FFC-dev mailing list
> FFC-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/ffc-dev