Universite-Gustave-Eiffel / elastodynamicsx

Elastodynamics with FEniCSx/DOLFINx
https://universite-gustave-eiffel.github.io/elastodynamicsx/
MIT License
4 stars 0 forks source link

Installation issue #1

Closed mlehujeur closed 1 year ago

mlehujeur commented 1 year ago

Bonjour Pierric,

Merci pour le partage de ton code, J'ai essayé d'installer le programme selon l'option 2 (Inside a FEniCSx Docker image / For use with shell commands)

Et j'obtiens l'erreur suivante (depuis le container) :

root@PNAN-20-123:~# python3 examples/weq_2D-SH_FullSpace.py 
Number of points per wavelength at central frequency:  2.5
Number of time steps per period at central frequency:  62.5
CFL condition: Courant number =  0.04
Traceback (most recent call last):
  File "/root/examples/weq_2D-SH_FullSpace.py", line 130, in <module>
    tStepper = TimeStepper.build(pde.m, pde.c, pde.k, pde.L, dt, V, bcs=bcs, scheme='leapfrog')
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/solvers/timestepper/timestepper.py", line 44, in build
    if scheme.lower() in s_.labels: return s_(*args, **kwargs)
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/solvers/timestepper/leapfrog.py", line 53, in __init__
    self._L =-dt_*dt_*k_(self._u_nm1, v) + 2*m_(self._u_nm1,v) - m_(self._u_nm2,v)
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/pde/pde.py", line 42, in <lambda>
    return lambda u,v: sum([mat.k(u,v) for mat in self.materials])
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/pde/pde.py", line 42, in <listcomp>
    return lambda u,v: sum([mat.k(u,v) for mat in self.materials])
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/pde/elasticmaterial.py", line 36, in k
    return self.k_DG(self._DGvariant)
  File "/usr/local/lib/python3.10/dist-packages/elastodynamicsx/pde/elasticmaterial.py", line 134, in k_DG
    n = ufl.FacetNormal(V)
  File "/usr/local/lib/python3.10/dist-packages/ufl/geometry.py", line 85, in __init__
    self._domain = as_domain(domain)
  File "/usr/local/lib/python3.10/dist-packages/ufl/domain.py", line 272, in as_domain
    return extract_unique_domain(domain)
  File "/usr/local/lib/python3.10/dist-packages/ufl/domain.py", line 345, in extract_unique_domain
    domains = extract_domains(expr)
  File "/usr/local/lib/python3.10/dist-packages/ufl/domain.py", line 338, in extract_domains
    for t in traverse_unique_terminals(expr):
  File "/usr/local/lib/python3.10/dist-packages/ufl/corealg/traversal.py", line 134, in traverse_unique_terminals
    for op in unique_pre_traversal(expr, visited=visited):
  File "/usr/local/lib/python3.10/dist-packages/ufl/corealg/traversal.py", line 70, in unique_pre_traversal
    visited.add(expr)
TypeError: unhashable type: 'FunctionSpace'

Est ce que ça te parle?

Merci,

Maximilien

mlehujeur commented 1 year ago

Je suis reparti d'une image antérieure de dolfinx/dolfinx en remplaçant la première ligne de Dockerfile.shell de FROM dolfinx/dolfinx:stable vers FROM dolfinx/dolfinx:v0.5.2

et ça a résolu mon bug... j'ai pu lancer le teste 2D SH! ++