Comments (2)
Using this approximation, I've observed convergence of ILUT preconditioned MCSA with 21 entries per row in the composite SPN operator, only 3 times as many as the original SPN operator, a big improvement over 2000 entries per row. This makes sense because if ILUT parameters were used to get a pretty good factorization, then those domain entries remaining after doing the approximation will therefore be the largest and have the greatest contribution to the Neumann series, effectively capturing most of the important pieces of the problem.
from chimera.
I've tested the reduced domain approximation with ParaSails preconditioning for SP1, 1-group problem 3. Here are the results with full weight recovery:
We note first that we can't reduce the fill level as low as we could with the ILUT preconditioning. ParaSails doesn't do as good of a job preconditioning in the first place so throwing away enough domain entries ultimately destroys convergence. The good news here is reducing the domain doesn't particularly reduce iterative performance when the method converges with an ideal fill level of 200 entries per row.
Now the weight recovery is set to 0:
This helps iterative performance a tremendous amount and now we can converge with a significantly reduced domain (at the cost of iterative performance). This is different than the initial ILUT observations where here zero weight recovery aids robustness. My thinking here is that the reduced domain from a more poorly conditioned problem (with ParaSails in this case), should not be pushed as hard in the Monte Carlo and hence reducing the size of the contributions to the estimator by reducing weights gives a similar effect to using a Neumann relaxation parameter of less than 1.
from chimera.
Related Issues (20)
- Prelim Exam Slides for Monte Carlo Linear Methods
- Prelim Exam Slides for Nonlinear Methods
- Prelim Exam Slides for Research Plan
- Add MCLS to Denovo Stratimikos solvers. HOT 3
- Design a series of computational experiments for domain overlap studies.
- Design a series of computational experiments for domain replication studies.
- Build MCLS/Denovo on Titan
- Spectral Analysis of SPn Equations HOT 10
- Drekar Access HOT 2
- Implement Point/Block Jacobi Preconditioning HOT 10
- SPN Convergence Analysis for Reactor Problems HOT 36
- Expected Value Estimator Anlaysis HOT 6
- ParaSails Preconditioner Analysis for the SPN Equations HOT 5
- Proof-of-Principle Calculation for the FANM Method HOT 2
- Nonlinear Benchmark 1: Thermal Convection Cavity
- Nonlinear Benchmark 2: Lid-Driven Cavity
- Nonlinear Benchmark 3: Backward-Facing Step
- Asymmetric Benchmark 1: Problem 3
- Asymmetric Benchmark 2: Problem 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from chimera.