Dear Patrick,
Thanks for your help!
It essentially worked, however, for some reason the nonlinear solver does not allow me to use direct methods as the inner solver, i.e.,
('ksp_type': 'preonly',
'pc_type': 'lu',
'pc_factor_mat_solver_type': 'mumps')
produces some sort of error, possibly related to local_to_global_map, I do not understand the error actually.
Thanks anyway!
Send firedrake mailing list submissions to
firedrake@imperial.ac.uk
To subscribe or unsubscribe via the World Wide Web, visit
https://mailman.ic.ac.uk/mailman/listinfo/firedrake
or, via email, send a message with subject or body 'help' to
firedrake-request@imperial.ac.uk
You can reach the person managing the list at
firedrake-owner@imperial.ac.uk
When replying, please edit your Subject line so it is more specific
than "Re: Contents of firedrake digest..."
Today's Topics:
1. BC's - 1D Hermite elements (Roberto Federico Ausas)
2. Re: BC's - 1D Hermite elements (Patrick Farrell)
3. Problem installing Firedrake (Smail Merabet)
----------------------------------------------------------------------
Message: 1
Date: Mon, 6 Apr 2020 16:49:25 -0300
From: Roberto Federico Ausas <rfausas@icmc.usp.br>
To: firedrake@imperial.ac.uk
Subject: [firedrake] BC's - 1D Hermite elements
Message-ID:
<CAH6Ec4nWnhmM=HeMtd8uU9VQ9NF0_hVdorMFdK=p8Xc3Gus4xQ@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hello,
My name is Roberto. I'm currently a professor in the Applied Mathematics
Department of Univ. of S?o Paulo at S?o Carlos.
I'm working on some fluid-structure interaction problems in which the solid
will be modeled as a 1D beam. Our reference integration domain will be the
interval [0,L].
We are planning to use firedrake and Hermite elements. The problem is
nonlinear.
I already have a working example. Values degrees of freedom are easy to
impose by using DirichletBC(), however I'm not figuring out how to enforce
derivatives degrees of freedom as boundary conditions (let say at x=0),
which I need for a clamped beam.
Any help or trick would be appreciated
Regards
--
Roberto F. Ausas
Instituto de Ci?ncias Matem?ticas e de Computa??o - ICMC/USP
S?o Carlos, SP - Brasil
*http://www.lmacc.icmc.usp.br/~ausas* <http://www.lmacc.icmc.usp.br/~ausas>
-------------- next part --------------
HTML attachment scrubbed and removed
------------------------------
Message: 2
Date: Mon, 6 Apr 2020 20:57:06 +0100
From: Patrick Farrell <patrick.farrell@maths.ox.ac.uk>
To: firedrake@imperial.ac.uk
Subject: Re: [firedrake] BC's - 1D Hermite elements
Message-ID: <4a52e3a8-1bca-249c-b9ff-304237289370@maths.ox.ac.uk>
Content-Type: text/plain; charset=utf-8; format=flowed
Hi Roberto,
I did this before. I projected some functions onto 1D Hermite spaces to see
which dofs are the derivative value at the endpoints. In the end I implemented
class FixHermiteDerivativeValue(DirichletBC):
@utils.cached_property
def nodes(self):
if V.mesh().mpi_comm().size > 1: raise NotImplementedError
return [1] # the derivative component of the left endpoint
bcs = [FixHermiteDerivativeValue(V, 0, "on_boundary")]
This only fixes the derivative at the left; I guess the derivative on the right will be [-1],
but am not sure. The "on_boundary" is ignored; I just kept it to not have to override the
constructor of DirichletBC.
Hope this helps,
Patrick
On 06/04/2020 20:49, Roberto Federico Ausas wrote:
> Hello,
>
> My name is Roberto. I'm currently a professor in the Applied Mathematics Department of Univ. of S?o Paulo at S?o Carlos.
>
> I'm working on some fluid-structure interaction problems in which the solid will be modeled as a 1D beam. Our reference integration domain will be the interval [0,L].
>
> We are planning to use firedrake and Hermite elements. The problem is nonlinear.
>
> I already have a working example. Values degrees of freedom are easy to impose by using DirichletBC(), however I'm not figuring out how to enforce derivatives degrees?of freedom as boundary conditions (let say at x=0), which I need for a clamped beam.
>
> Any help or trick would be appreciated
> Regards
>
> --
> Roberto F. Ausas
> Instituto de Ci?ncias Matem?ticas e de Computa??o- ICMC/USP
> S?o Carlos, SP - Brasil
> _http://www.lmacc.icmc.usp.br/~ausas_
>
> _______________________________________________
> firedrake mailing list
> firedrake@imperial.ac.uk
> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
>
------------------------------
Message: 3
Date: Tue, 7 Apr 2020 08:10:54 +0100
From: Smail Merabet <merabetsmail@gmail.com>
To: firedrake@imperial.ac.uk
Subject: [firedrake] Problem installing Firedrake
Message-ID:
<CAEo0a_H1WQ=jFvkjeW5RiT4KkKM7PhF3-ya2X2pE2mjggKopeA@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hi firedrake users,
I have a problem with the installation of firedrake,
I work on Mojave 10.14.6
I did the following steps:
url -O
https://raw.githubusercontent.com/firedrakeproject/firedrake/master/scripts/firedrake-install
unset PYTHONPATH
unset PETSC_DIR
unset PETSC_ARCH;
python3 firedrake-install
the installation see to be done, I attached the log file to this email, but
when I excuse any script I get the following error message
(firedrake) (base) Merabets-MacBook-Pro:examples merabet$ python3 fi.py
Traceback (most recent call last):
File "/Users/merabet/firedrake/src/h5py/h5py/__init__.py", line 26, in
<module>
from . import _errors
ImportError: dlopen(/Users/merabet/firedrake/src/h5py/h5py/_
errors.cpython-37m-darwin.so, 2): Library not loaded:
@rpath/libhdf5.103.dylib
Referenced from: /Users/merabet/firedrake/src/h5py/h5py/_
errors.cpython-37m-darwin.so
Reason: image not found
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "fi.py", line 1, in <module>
import firedrake
File "/Users/merabet/firedrake/src/firedrake/firedrake/__init__.py", line
53, in <module>
from firedrake.checkpointing import *
File "/Users/merabet/firedrake/src/firedrake/firedrake/checkpointing.py",
line 3, in <module>
from firedrake.cython import hdf5interface as h5i
File "firedrake/cython/hdf5interface.pyx", line 2, in init
firedrake.cython.hdf5interface
File "/Users/merabet/firedrake/src/h5py/h5py/__init__.py", line 30, in
<module>
raise ImportError("You cannot import h5py from inside the install
directory.\nChange to another directory first.")
ImportError: You cannot import h5py from inside the install directory.
Change to another directory first.
(firedrake) (base) Merabets-MacBook-Pro:examples merabet$
Please, could any one help me to solve this problem ?
Best regards
-------------- next part --------------
HTML attachment scrubbed and removed
-------------- next part --------------
A non-text attachment was scrubbed...
Name: firedrake-install.log
Type: application/octet-stream
Size: 6649021 bytes
Desc: not available
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20200407/8cfd0625/attachment-0001.obj>
------------------------------
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake
End of firedrake Digest, Vol 40, Issue 4
****************************************
--
Roberto F. Ausas
Instituto de Ciências Matemáticas e de Computação - ICMC/USP
São Carlos, SP - Brasil