Banyak pustaka FEM adaptif menggunakan struktur data mesh yang lebih canggih untuk menangani penambahan / penghapusan node, edge, triangles, tetrahedra, dll. Misalnya, p4est library menggunakan struktur data octree untuk perbaikan mesh adaptif; Anda tidak akan sering menemukan oktri yang digunakan untuk perhitungan pada jaring statis.
Apa perubahan pada sisi aljabar linier untuk FEM adaptif?
Cara paling tumpul yang bisa saya bayangkan adalah membangun kembali semua matriks sistem sepenuhnya setiap kali mesh tersebut disempurnakan atau kasar. Jika adaptasi mesh adalah operasi yang cukup jarang, maka biaya melakukannya pada akhirnya diamortisasi selama sisa perhitungan. Seseorang dapat dengan mudah memanfaatkan perangkat lunak aljabar linier jarang yang ada (PETSc, Trilinos, dll.) Dengan pendekatan ini.
Apakah metode blunt ini paling umum digunakan, atau apakah ada perpustakaan yang berhasil menggunakan kembali atau memodifikasi matriks lama selama penyempurnaan? Bagaimanapun, sebagian besar mesh dan matriks yang sesuai tidak berubah selama adaptasi mesh.
sumber