`allocate_matrix` with only a topology gives a bad error message

Repository metrics

Stars

439 stars (439 stars)

PR merge metrics

PR metrics pending (PR metrics pending)

Description

const grid = generate_grid(Hexahedron, (n, n, n))
const topology = ExclusiveTopology(grid)

const dh = DofHandler(grid)
const dh = DofHandler(grid)
add!(dh, :u, Lagrange{RefHexahedron, 2}()^3)
add!(dh, :p, Lagrange{RefHexahedron, 1}())
add!(dh, :q, DiscontinuousLagrange{RefHexahedron, 1}())
close!(dh)

allocate_matrix(dh; topology)

# ERROR: TypeError: in keyword argument interface_coupling, expected AbstractMatrix{Bool}, got a value of type Nothing
# Stacktrace:
#  [1] add_sparsity_entries!(sp::SparsityPattern, dh::DofHandler{…}, ch::Nothing; keep_constrained::Bool, coupling::Nothing, interface_coupling::Nothing, topology::ExclusiveTopology)

Contributor guide

Research direction

Trace the code path for `allocate matrix` with a topology keyword to identify where `Nothing` is passed as `interface coupling` and add a check for a valid value.

Tech stack

julia

Domain

backend

Issue type

Bug

DifficultyEstimated implementation difficulty for a new contributor, from 1 for very small changes to 5 for expert-level work.

2

Estimated timeA rough time range for an experienced contributor to investigate, implement, test, and prepare a pull request.

1-3 hours

Activity statusHow available the issue appears right now: fresh, active, stale, blocked, or waiting on maintainer input.

Fresh

ClarityHow clearly the issue explains the expected change, acceptance criteria, and next step.

Clear

Prerequisites

Julia

Newbie friendlinessA 1-100 score estimating how approachable this issue is for first-time contributors.

75