> Hi,
>
> I'm trying to solve a two-dimensional partial differential equation using a
> DG FEM. The differential equation involves an integral term in only one
> direction of the mesh:
> \int_{-2}^{2} b(x,y,y') u(x,y') dy'
> with the integral kernel b. Is there a simple way to build the corresponding
> matrix of the double integral which arises during the FEM discretisation
> using the usual assembling procedures?
I do not see a simple way (that is a solution that can be expressed
purely as an argument to generic_assembly::set() method).
My guess would be to do the calculations in properly implemented
nonlinear term that would be inserted into assembly chain with
generic_assembly::push_nonlinear_term(pnonlinear_elem_term) but
probably you have figured it out yourself.
Regards,
Roman
--
Roman Putanowicz, PhD < [email protected] >
Institute for Computational Civil Engng (L-5)
Dept. of Civil Engng, Cracow Univ. of Technology
www.l5.pk.edu.pl, tel. +48 12 628 2569, fax 2034
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
