Skip to content

Set NumberOfComponents to 3 for a VFS. - #260

Merged
ctjacobs merged 1 commit into
masterfrom
fix-parallel-vfs-write
Apr 14, 2014
Merged

Set NumberOfComponents to 3 for a VFS.#260
ctjacobs merged 1 commit into
masterfrom
fix-parallel-vfs-write

Conversation

@ctjacobs

Copy link
Copy Markdown
Contributor

Set NumberOfComponents to 3 in the pvtu file when writing a Function whose function space is a VectorFunctionSpace.

Fixes #259.

@kynan

kynan commented Apr 14, 2014

Copy link
Copy Markdown
Contributor

Have squashed the fixup commit and used isinstance instead of a type check. I think this is good to go.

ctjacobs added a commit that referenced this pull request Apr 14, 2014
@ctjacobs
ctjacobs merged commit 27c32e3 into master Apr 14, 2014
@ctjacobs
ctjacobs deleted the fix-parallel-vfs-write branch April 14, 2014 11:37
connorjward pushed a commit that referenced this pull request Dec 4, 2024
…pe-inference

Faster type inference in loopy codegen
hardik-corintis added a commit to Corintis/firedrake that referenced this pull request Aug 7, 2026
Threads dtype=RealType (or the numpy-derived real dtype from scalar_type
in TSFC) through every as_fiat_cell/create_element call site in
firedrake and tsfc, using FIAT PR firedrakeproject#260 and FInAT's new dtype support
(firedrakeproject/fiat#263). Removes the finat.element_factory.ufc_cell
monkeypatch in firedrake/utils.py, which is no longer needed now that
callers declare their working precision explicitly instead of relying
on a globally-overridden default.

Also fixes petsc_sparse's hardcoded rtol=1E-10 drop tolerance
(firedrake/preconditioners/fdm.py), which is too tight for single
precision: reference-cell geometry feeding into these matrices is only
accurate to float32 round-off (~1e-7), so entries that should be
exactly zero (e.g. in a Nedelec discrete gradient matrix) were leaking
into the sparsity pattern above the fp64-tuned threshold, corrupting
the exact combinatorial structure PETSc's PCBDDCNedelecSupport expects
and causing test_bddc_aij_simplex[N1curl-3-False] to fail.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Parallel writing of a Function in a VFS is broken

2 participants