Re: [firedrake] Solve a Variational problem in a part of the domain
Hi Anna, Not as written, no. I think that we would need to design a bespoke solver for that. all the best --cjc On 17 August 2015 at 12:18, Anna Kalogirou <a.kalogirou@leeds.ac.uk> wrote:
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> <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 listfiredrake@imperial.ac.ukhttps://mailman.ic.ac.uk/mailman/listinfo/firedrake
--
Dr Anna Kalogirou Research Fellow School of Mathematics University of Leeds
-- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916
Yes, obviously I wouldn't be able to write assemble(u*dx) for the trial function u, but I just wanted to show you what I have. Is it too difficult to design such a non-conventional solver? Thanks, Anna. On 17/08/15 13:04, Colin Cotter wrote:
Hi Anna, Not as written, no. I think that we would need to design a bespoke solver for that. all the best --cjc
On 17 August 2015 at 12:18, Anna Kalogirou <a.kalogirou@leeds.ac.uk <mailto:a.kalogirou@leeds.ac.uk>> wrote:
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> <mailto: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 <mailto: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/ <http://www1.maths.leeds.ac.uk/%7Ematak/>
-- http://www.imperial.ac.uk/people/colin.cotter
www.cambridge.org/9781107663916 <http://www.cambridge.org/9781107663916>
_______________________________________________ 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/
participants (2)
-
Anna Kalogirou
-
Colin Cotter