Apa yang dimaksud dengan prasyarat terukur untuk Helmholtz frekuensi tinggi?

15

Metode multigrid dan dekomposisi domain standar tidak berfungsi, tetapi saya memiliki masalah 3D besar dan pemecah langsung bukan merupakan pilihan. Metode apa yang harus saya coba?

Bagaimana pilihan saya dipengaruhi oleh pertimbangan berikut?

  • koefisien bervariasi pada beberapa urutan besarnya, atau
  • elemen hingga versus metode yang berbeda hingga digunakan
Jed Brown
sumber
1
Dalam 3D, pemecah iteratif biasanya berkinerja buruk, saya sarankan merujuk pada beberapa pemecah langsung HSS-matrix reordering dari Ming Gu, Xia, dan Chandrasekaran.
Shuhao Cao

Jawaban:

5

Saya pikir secara umum perlu diingat bahwa metode yang paling efisien yang kita miliki (geometri dan aljabar multigrid serta, sampai taraf tertentu, dekomposisi domain) bergantung pada kenyataan bahwa solusi PDE seringkali mulus dan penyelesaian masalah yang lebih kasar dapat menghasilkan perkiraan yang baik untuk masalah skala halus. Masalah dengan persamaan Helmholtz untuk frekuensi tinggi adalah bahwa asumsi ini tidak benar: Anda memang membutuhkan mesh yang relatif baik untuk mewakili solusinya, dan pemecah mesh kasar tidak akan dapat menghasilkan apa pun yang banyak digunakan. Akibatnya, pendekatan khas untuk prekondisi yang baik tidak berfungsi dalam kasus itu, dan itulah alasan yang mendasari mengapa tidak ada opsi yang benar-benar bagus dalam kasus Anda kecuali hanya melempar banyak prosesor pada masalah;

Wolfgang Bangerth
sumber
4

Hal-hal H-matrix dari Jack Poulson dan Lexing Ying adalah metode paling efisien yang saya tahu. Ini harus dirilis pada musim semi tetapi mereka telah memberikan presentasi tentang itu.

Matt Knepley
sumber
2
Saya mungkin harus memenuhi syarat pernyataan Anda dengan mengatakan bahwa kami telah memecahkan masalah frekuensi tinggi skala besar dengan pendekatan PML yang bergerak secara efisien pada ribuan core untuk frekuensi yang sangat tinggi, tetapi kami belum menguji strategi matriks-H pada skala itu. Alasannya adalah bahwa ia tidak memiliki banyak pembenaran teoretis dalam 3d, meskipun fakta bahwa itu akan lebih terukur dari perspektif komputasi paralel.
Jack Poulson