Compatible PRanges between FESpaces and linear systems
@JordiManyer ya está trabajando en esto.
Desde el 30/11/2023.
Evaluación
Este issue todavía no se ha evaluado.
Descripción
Here are some notes on the issue:
- In some specific scenarios (which are not as rare), owned rows of the matrix contain column identifiers from dofs which are NOT available in the local FESpaces. This means that we will never be able to 100% match the PRanges from the PSparseMatrix and the DistributedFESpace. (See Alberto's note)
- Moreover, when we assemble ghost rows we also include some column ids which are NOT needed to multiply the matrix by a column vector. We need this extra rows as cache for in-place assembly, but these extra ghosts should not be included when communicating during matrix-vector product.
- On top of this, it would be advantageous to reorder the dofs in the local FESpaces so that the dof layout contains owned ids first, then ghost ids. This is how we do things in the linear system, and would bring us on par with other softwares like PETSc.
Given these observations, we concluded the following:
- We need to re-order the ids in the local FESpaces. Hopefully, after doing so,
LocalViewtype wont be needed anymore. - We need to create a new type
FECompatiblePSparseMatrix, that contains- The usual matrix
- Caches for the ghosts values that the matrix column PRange needs
- Separate PRanges for matrix-vector multiplication and for assembly
- Lenguaje dominante
- Julia
- Estrellas
- 150
- Forks
- 28
- Métricas de merge de PR
- Sin PR fusionados en 30 d
Preparar el entorno
Este proyecto no incluye contenedor de desarrollo, Dockerfile ni guía de contribución, así que la configuración corre por tu cuenta: empieza por su README y consulta nuestra guía para la primera contribución para los pasos generales.
Primeros pasos
- Lee el issue completo y luego la guía de contribución del proyecto.
- Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
- Haz un fork del repositorio y trabaja en una rama.
- Abre un pull request que haga referencia al número del issue.
Más de gridap/GridapDistributed.jl
-
A generic method of the abstract type DistributedDiscreteModel{Dc,Dp} directly accesses to a field name .face_gidsQuizá libre de nuevo @amartinhuertas la tomó hace 121 días y no hay ningún pull request abierto. Abierto
gridap/GridapDistributed.jl#211 · 1 asignado ·
-
Automatically add ghost cells to triangulation whenever a CellField is created out of a DistributedTriangulationQuizá libre de nuevo @JordiManyer la tomó hace 154 días y no hay ningún pull request abierto. Abierto
gridap/GridapDistributed.jl#209 · 2 asignados ·
-
Dificultad 4/5 3-5 días Aptitud para principiantes 35/100
gridap/GridapDistributed.jl#189 · 1 comentario ·
-
Issue with DG-SIP and FullyAssembledRowsQuizá libre de nuevo @amartinhuertas la tomó hace 359 días y no hay ningún pull request abierto. Abierto
gridap/GridapDistributed.jl#186 · 20 comentarios · 1 asignado ·
-
Potential issue when interpolating DistributedCelllFields that are only defined on own cellsQuizá libre de nuevo @JordiManyer la tomó hace 373 días y no hay ningún pull request abierto. Abierto
gridap/GridapDistributed.jl#184 · 1 comentario · 2 asignados ·
Todos los issues de gridap/GridapDistributed.jl
Issues similares
-
documentation invalid
Dificultad 1/5 1-3 horas Aptitud para principiantes 88/100
JuliaPhysics/BeamletOptics.jl#190 ·
Los mantenedores suelen responder en 2 días
-
mul! with a Symmetric or Hermitian factor throws for a matrix view whose columns are not contiguousAbierto
Dificultad 2/5 1-3 horas Aptitud para principiantes 72/100
JuliaLang/LinearAlgebra.jl#1829 ·
Los mantenedores suelen responder en 1 día
-
`EnsembleSummary` `ci_type = :SEM` ribbon is `sqrt(N)` too wide for scalar statesPosiblemente ocupada @isaacsas la tomó hoy. Abierto
Dificultad 2/5 1-3 horas Aptitud para principiantes 88/100
SciML/SciMLBase.jl#1654 ·
Los mantenedores suelen responder en 3 días
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
MagneticResonanceImaging/MRIReco.jl#320 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
rfordatascience/tidytuesday#1098 ·