← Back to team overview

ffc team mailing list archive

Re: [FIAT-dev] evaluate_basis and derivatives

 

Quoting Robert Kirby <robert.c.kirby@xxxxxxxxx>:

> Yes.  The columns of the coeff matrix correspond to the orthogonal basis
> functions.You need those to construct the 9 "non-bubble" bases.

OK, that's was how I understood things too, thanks.

Kristian

>
>
> On Tue, Aug 18, 2009 at 2:37 PM, Anders Logg <logg@xxxxxxxxx> wrote:
>
> > On Tue, Aug 18, 2009 at 09:12:49PM +0200, Kristian Oelgaard wrote:
> > >
> > >
> > > Hello,
> > >
> > > I want to eliminate a basis function value (i.e. not generate code for
> > it) for
> > > evaluate_basis() (and evaluate_basis_derivative()).
> > >
> > > Calling tabulate(x) on the FIAT element FiniteElement("CG", "triangle",
> > 3)
> > > will return the dot product of the expansion coefficients (a matrix) and
> > the
> > > bvals (an array) which is evaluated at the point x.
> > >
> > > In the current case the coefficients is a 10x10 matrix and the bvals is a
> > 10x1.
> > > If I just want to return the first 9 basis function values should I
> > remove the
> > > last row of the coefficient matrix only, or should I also remove the last
> > > column as well along with the last row in the bvals?
> > >
> > > Kristian
> >
> > I think you should remove the last row to get a 9x10 matrix which you
> > multiply with a 10x1 vector to get a 9x1 vector with the values.
> >
> > --
> > Anders
> >
> > -----BEGIN PGP SIGNATURE-----
> > Version: GnuPG v1.4.9 (GNU/Linux)
> >
> > iEYEARECAAYFAkqLAu4ACgkQTuwUCDsYZdEkfQCfYwIjbuiAeTS4gCUsLb4yqXuz
> > NF4An3HHeqXGfxRcwWOFsW6qBgahlV8x
> > =BxOj
> > -----END PGP SIGNATURE-----
> >
> > _______________________________________________
> > FFC-dev mailing list
> > FFC-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/ffc-dev
> >
> >
>




References