ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, May 25, 2007 at 09:01:04AM +0200, Garth N. Wells wrote:
> Anders Logg wrote:
> > 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 find an old changeset, but it still didn't work. I think that 
> the example has an extra factor in it (which I didn't have before) which 
> breaks things. It does work with quadrature.
> 
> Garth

Could it just be that numpy is generating the error since FFC tries to
create an array that is too large? (Product of all dimensions.)

/Anders


> > 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
> > _______________________________________________
> > 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