discretize.TreeMesh.average_edge_y_to_cell

TreeMesh.average_edge_y_to_cell

Averaging operator from y-edges to cell centers (scalar quantities).

This property constructs a 2nd order averaging operator that maps scalar quantities from y-edges to cell centers. This averaging operator is used when a discrete scalar quantity defined on y-edges must be projected to cell centers. Once constructed, the operator is stored permanently as a property of the mesh. See notes.

Returns:

(n_cells, n_edges_y) scipy.sparse.csr_matrix

The scalar averaging operator from y-edges to cell centers

Notes

Let \(\boldsymbol{\phi_y}\) be a discrete scalar quantity that lives on y-edges. average_edge_y_to_cell constructs a discrete linear operator \(\mathbf{A_{yc}}\) that projects \(\boldsymbol{\phi_y}\) to cell centers, i.e.:

\[\boldsymbol{\phi_c} = \mathbf{A_{yc}} \, \boldsymbol{\phi_y}\]

where \(\boldsymbol{\phi_c}\) approximates the value of the scalar quantity at cell centers. For each cell, we are simply averaging the values defined on its y-edges. The operation is implemented as a matrix vector product, i.e.:

phi_c = Ayc @ phi_y