Je cherche à résoudre un problème qui est une PDE, et le domaine discrétisé 3D peut avoir une condition aux limites différente sur chacune des 6 limites (ou tout de même).Comment stockez-vous un domaine 3D discrétisé (pour résoudre PDE) sous un format clairsemé, lorsque les conditions aux limites peuvent changer?
Quelle est la meilleure façon de mettre cette matrice clairsemée dans un format compressé? Est-ce que la RSE sera ma seule option ici? J'ai pensé utiliser ellpack, mais je ne suis pas sûr de savoir comment cela fonctionne pour changer les conditions aux limites. En pensant à la représentation matricielle 2D de l'espace 3d ... il sera principalement diagonalement dominant avec 7 diagonales, mais ces diagonales changeront potentiellement le long des limites. Ne semble pas que je peux utiliser un format où je stocke les valeurs, et un décalage de la diagonale qui est la même à chaque fois.
Il est évident que je suis en train de régler ce problème pour être plus convivial pour mon cache solveur de CG qui fait beaucoup de vecteur-matrice multiplier