Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Doc : how to build custom operator #57

Merged
merged 9 commits into from
Mar 29, 2024
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
rename sketch
  • Loading branch information
mbouyges committed Mar 27, 2024
commit 1d62c236508a8f4aa5a17d7c59bf1b264671b0df
2 changes: 1 addition & 1 deletion docs/src/manual/geometry.md
Original file line number Diff line number Diff line change
Expand Up @@ -10,4 +10,4 @@ F(\xi) = \sum_{i=1}^n \hat{\lambda}_i(\xi)M_i

where $(\lambda)_i$ are the Lagrange polynomials whose order matches the element order.

![sketch](../assets/sketch/sketch-integration.png)
![sketch](../assets/sketch/sketch-mapping.png)
4 changes: 2 additions & 2 deletions docs/src/manual/integration.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,10 +8,10 @@ To compute an integral on a geometrical element, for instance a curved element,

where we recall that $$F$$ is the reference to physical mapping and $$J$$ is the determinant of the jacobian matrix of this mapping. Depending on the shape and element order, this determinant is either hard-coded or computed with `ForwardDiff`.

Now, to compute the right side, i.e the integral on the reference shape, quadrature rules are applied:
Now, to compute the right side, i.e the integral on the reference shape, quadrature rules are applied to $\hat{g} = g \circ F$:

```math
\int_{\hat{\Omega}} g(\hat{x}) \mathrm{\,d} \hat{\Omega} = \sum_{i =1}^{N_q} \omega_i g(\hat{x}_i)
\int_{\hat{\Omega}} \hat{g}(\hat{x}) \mathrm{\,d} \hat{\Omega} = \sum_{i =1}^{N_q} \omega_i \hat{g}(\hat{x}_i)
```

A specific procedure is applied to compute integrals on a face of a cell (i.e a surfacic integral on a face of a volumic element).
2 changes: 1 addition & 1 deletion src/mesh/entity.jl
Original file line number Diff line number Diff line change
Expand Up @@ -512,7 +512,7 @@ struct isVolumic <: TopologyStyle end
"""
topology_style(::AbstractEntityType{topoDim}, ::Node{spaceDim, T}) where {topoDim, spaceDim, T}
topology_style(::AbstractEntityType{topoDim}, ::AbstractArray{Node{spaceDim, T}, N}) where {spaceDim, T, N, topoDim}
)


Indicate the `TopologyStyle` of an entity of topology `topoDim` living in space of dimension `spaceDim`.
"""
Expand Down
Loading