220 likes | 293 Views
Fast Method of Fundamental Solution. Xinrong Jiang, PhD candidate Wen Chen, Prof. C.S. Chen, Prof. NTU, Dec. 8, Taipei. Hohai University. Outline. Motivation Methodology Numerical Example Conclusion. Motivation. Radial Basis Functions (RBFs) : Domain Type Kansa’s Method
E N D
Fast Method of Fundamental Solution Xinrong Jiang, PhD candidate Wen Chen, Prof. C.S. Chen, Prof. NTU, Dec. 8, Taipei Hohai University
Outline • Motivation • Methodology • Numerical Example • Conclusion
Motivation • Radial Basis Functions (RBFs) : • Domain Type • Kansa’s Method • Local Method of Particular Solution (LMPS) • Boundary Type • Method of Fundamental Solution (MFS) • Regularized Meshless Method (RMM) • Boundary Knot Method (BKM) • Boundary Particle Method (BPM) • Singular Boundary Method (SBM)
Motivation • Dense Matrix • Large-scale problem • Infinite domain … • Speed up ^_^ iterative method • Mainframe T_T • expensive, computer volume large, electricity • PC T_T • Memory, CPU flops
Motivation • Fast algorithm • Save memory • Fast computing • Efficiency Accuracy
Methodology • Fast Multipole Method (FMM) • 1987, 1997 new version, Rokhlin and Greengard • Nlog(N)->N • Adaptive • Fast Fourier Transform (FFT) • Precorrected FFT, J White • N*log(N) • Uniform • Hierarchical Matrix,Adaptive Cross-Approximation, etc
Methodology collocation, source Points: N Matrix: N*N
Methodology • Iterative Method • ill-conditioned
Methodology • Krylov Subspace method::GMRES • Generalized minimal residual method • FMM-BEM • Fail in FMM-MFS • Iterator: open issue
MFS easy-to-program, exponential convergence, highly accuracy, geometric flexibility and so on infinite domain problems, large deformation problems, dynamic crack propagation etc FMM-MFS
GPU • Further acceleration
GPU::CUDA • CUDA: Compute Unified Device Architecture
Parallel • Tree Structure • Further Study
Conclusion • FMM-MFS • implement successfully • high precision and speed • high wave number requires wide band • high frequency • large domain with low frequency • ill-conditioned requires suitable iterative method
End Thanks for your attention! Comments? 姜欣榮 Xinrong Jiang hhujiangxr@163.com 08/12/ 2011