Derivative of JacobianInverse fails

Chiusa
  • Qualcosa non va?
#445 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

I maintainer di solito rispondono entro 1 giorno

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
4/5
Tempo stimato
3-5 giorni
Idoneità per principianti
25/100
Tipo di issue
Bug
Chiarezza
Da chiarire
Stato di attività
Ferma
Stack tecnologico
python
Ambito
compilers

Direzione di ricerca

Riprodurre il fallimento con l'esempio Python, quindi esaminare ffcx/ir/elementtables.py, in particolare get_ffcx_table_values e build_optimized_tables, insieme a ffcx/ir/integral.py. Tracciare come viene prodotto derivative_counts per ufl.geometry.JacobianInverse e confrontarlo con le dimensioni delle tabelle disponibili. Il lavoro è completato quando la forma fornita viene compilata e assemblata senza l'IndexError.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

import numpy as np

import ufl
import dolfinx.fem
from dolfinx.mesh import create_mesh
from mpi4py import MPI


domain = ufl.Mesh(ufl.VectorElement("Lagrange", ufl.Cell("interval", 3), 2))
mesh = create_mesh(MPI.COMM_WORLD, np.array([[0, 2, 1], [2, 4, 3]]), np.array(
    [[0.0, 0.0, 0.0], [0.25, 0.0, 0.0], [0.5, 0.0, 0.0], [0.75, 0.0, 0.0], [1.0, 0.0, 0.0]]), domain)

B = dolfinx.fem.TensorFunctionSpace(mesh, ("DG", 2), (mesh.topology.dim, mesh.topology.dim))
K = ufl.geometry.JacobianInverse(mesh).T
S = dolfinx.fem.form(ufl.inner(ufl.grad(K), ufl.grad(K)) * ufl.dx(mesh))

s = dolfinx.fem.assemble_scalar(S)

fails with

...
  File "/usr/local/lib/python3.9/dist-packages/dolfinx/jit.py", line 209, in ffcx_jit
    r = ffcx.codegeneration.jit.compile_forms([ufl_object], parameters=p_ffcx, **p_jit)
  File "/home/ffcx/ffcx/codegeneration/jit.py", line 168, in compile_forms
    impl = _compile_objects(decl, forms, form_names, module_name, p, cache_dir,
  File "/home/ffcx/ffcx/codegeneration/jit.py", line 232, in _compile_objects
    _, code_body = ffcx.compiler.compile_ufl_objects(ufl_objects, prefix=module_name, parameters=parameters)
  File "/home/ffcx/ffcx/compiler.py", line 103, in compile_ufl_objects
    ir = compute_ir(analysis, object_names, prefix, parameters, visualise)
  File "/home/ffcx/ffcx/ir/representation.py", line 94, in compute_ir
    irs = [
  File "/home/ffcx/ffcx/ir/representation.py", line 95, in <listcomp>
    _compute_integral_ir(fd, i, analysis.element_numbers, integral_names, finite_element_names,
  File "/home/ffcx/ffcx/ir/representation.py", line 360, in _compute_integral_ir
    integral_ir = compute_integral_ir(itg_data.domain.ufl_cell(), itg_data.integral_type,
  File "/home/ffcx/ffcx/ir/integral.py", line 84, in compute_integral_ir
    mt_table_reference = build_optimized_tables(
  File "/home/ffcx/ffcx/ir/elementtables.py", line 346, in build_optimized_tables
    t = get_ffcx_table_values(quadrature_rule.points, cell,
  File "/home/ffcx/ffcx/ir/elementtables.py", line 108, in get_ffcx_table_values
    tbl = tbl[basix_index(*derivative_counts)]
IndexError: index 4 is out of bounds for axis 0 with size 3

@mscroggs Any idea why most likely derivative_counts is wrong for this (a bit exotic) ufl symbol?

Lingua principale
Python
Stelle
196
Fork
47
Merge medio
2g 9h
PR unite (30g)
26

Preparare l'ambiente

Questo progetto non fornisce container di sviluppo, Dockerfile né guida per i contributori, quindi l'ambiente è a tuo carico: parti dal suo README e consulta la nostra guida al primo contributo per i passaggi generali.

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di FEniCS/ffcx

Tutte le issue di FEniCS/ffcx

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.