Hi all,

I have to solve a problem with the bilinear form similar to the following:

a = ( v*u + dt*inner(grad(u),grad(v)) )*dx + ( v*dt*assemble(u*dx) )*dx

The last term is essentially the product of the integral of u and the integral of test function v.
Is it even possible to solve this in Firedrake?

Best,

Anna.

On 14/08/15 15:43, Lawrence Mitchell wrote:

      

        
On 14 Aug 2015, at 11:55, Anna Kalogirou <a.kalogirou@leeds.ac.uk> wrote:

Hi,

I was planning to define the heavyside function as a DG0 function.

Any ideas about solving equation (1b), which contains a time-update for both a function and a scalar?


      

      
I haven't been fully following along, but it looks like the formulation in (1f) introduces, effectively, a global coupling.  Is that right?  Assuming you didn't have any implementation constraints, how would you go about solving this?

Lawrence

      

      

      

      
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake

    

    

    
-- 
 
 Dr Anna Kalogirou
 Research Fellow
 School of Mathematics
 University of Leeds

 http://www1.maths.leeds.ac.uk/~matak/