Skip to content

bug in create_interpolation_matrix on a complex environment #19

@pierricmora

Description

@pierricmora

Hi, and thanks again for this great package. I use complex numbers and I was getting crazy outputs until I found the following.

The _MPI.DOUBLE used to inform the block size in https://github.com/scientificcomputing/fenicsx_ii/blob/main/src/fenicsx_ii/interpolation.py#L124 is not suited to dtype=np.complex128. The effect is that recv_basis_functions is a shuffled version of basis_values_on_V.

I have proposed a fix on:
https://github.com/pierricmora/fenicsx_ii/blob/complex/src/fenicsx_ii/interpolation.py#L128

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions