Subject: Re: [Fiat] [Branch ~fiat-core/fiat/fiat-new] Rev 140: Try to fix bug for GaussJacobiQuadratureLineRule.

From: Kristian Oelgaard <[email-masked]>
Date: Tue, 12 Jan 2010 14:07:57 +0100 (CET)

2010/1/12  <[email-masked]>:
> ------------------------------------------------------------
> revno: 140
> committer: Kristian B. Ølgaard <[email-masked]>
> branch nick: fiat-new
> timestamp: Tue 2010-01-12 13:55:32 +0100
> message:
>  Try to fix bug for GaussJacobiQuadratureLineRule.
> modified:
>  FIAT_NEW/quadrature.py

Could someone please check if I got this right? (and delete the commented lines)

Kristian

>
> --
> lp:~fiat-core/fiat/fiat-new
> https://code.launchpad.net/~fiat-core/fiat/fiat-new
>
> You are subscribed to branch lp:~fiat-core/fiat/fiat-new.
> To unsubscribe from this branch go to https://code.launchpad.net/~fiat-core/fiat/fiat-new/+edit-subscription.
>
> === modified file 'FIAT_NEW/quadrature.py'
> --- FIAT_NEW/quadrature.py      2010-01-08 13:01:22 +0000
> +++ FIAT_NEW/quadrature.py      2010-01-12 12:55:32 +0000
> @@ -21,13 +21,15 @@
>  class GaussJacobiQuadratureLineRule( QuadratureRule ):
>     """Gauss-Jacobi quadature rule determined by Jacobi weights a and b
>     using m roots of m:th order Jacobi polynomial."""
> -    def __init__( self , ref_el , a , b , m ):
> +#    def __init__( self , ref_el , a , b , m ):
> +    def __init__( self , ref_el , m ):
>         # this gives roots on the default (-1,1) reference element
> -        (xs_ref,ws_ref) = compute_gauss_jacobi_rule( a , b , m )
> +#        (xs_ref,ws_ref) = compute_gauss_jacobi_rule( a , b , m )
> +        (xs_ref,ws_ref) = compute_gauss_jacobi_rule( 0. , 0. , m )
>
>         Ref1 = reference_element.DefaultLine()
> -        A,b = reference_element.make_affine_mapping( Ref1.get_vertices , \
> -                                                     ref_el.get_vertices )
> +        A,b = reference_element.make_affine_mapping( Ref1.get_vertices() , \
> +                                                     ref_el.get_vertices() )
>
>         mapping = lambda x: numpy.dot( A , x ) + b
>
>
>
>