Comments (9)
IT would be great to have .prm that demonstrates the solver not working very well.
from aspect.
It's difficult -- the testcase I have uses plugins for materials and geometry that rely on external files. I could make the geometry simpler, but in any case I doubt that the testcase could be made simple enough for the testsuite. All of the runs I have done take half an hour upward.
For (my own) reference, the testcase that triggered all of this is the 'fishwick' case of Sarah Stamps' runs.
from aspect.
I am not surprised. Maybe we can construct a simple example (before starting to work on AMG) that has a strong viscosity contrast and shows the same issues. It would be good to have a test case that is easy and doesn't require running 30mins+.
from aspect.
I'll see what I can do when I start to address this.
from aspect.
I read through the paper by May and Moresi and it addresses a different question (at least that's my take-away message): what best to do to precondition the Schur complement (we just use a scaled mass matrix). This is an issue that pertains to the number of outer iterations, whereas doing something about how well AMG solves the top left block is something that primarily pertains to inner iterations.
Of course, we should address both issues -- I have testcases that demonstrate that each of the two can be a problem.
from aspect.
Our experience has been that the advantage of the LSC preconditioner in May & Moresi is specific to the unstable (Q1-P0) discretization used by Underworld. With a stable discretization like Q2-P1, the scaled mass matrix has performed better. As for solving the viscous block using AMG, you can do it with GAMG or ML, but you need to provide the rigid body modes and setup will be expensive. And with the Q2 discretization, performance tends to be very sensitive to thresholding because the unthresholded graph would coarsen much too fast. Our SC14 paper demonstrates that there is a very large benefit to geometric coarsening and to the use of matrix-free methods on the finest level(s).
from aspect.
@jedbrown : Thanks for the insight. We do provide the rigid body modes to ML. I'm not so much worried about the cost of setup -- we seem to have cases where ML just turns out to be a not very good preconditioner at all for the top left block (of course these are exactly the cases with strongly varying viscosity), and that I wanted to improve.
I am ready to believe that geometric multigrid is very useful, but we don't currently have the software infrastructure to make that happen in parallel :-(
from aspect.
I am ready to believe that geometric multigrid is very useful, but we don't currently have the software infrastructure to make that happen in parallel :-(
"soon" :-)
One way to cheat might be to apply some kind of averaging of the viscosity inside each cell.
from aspect.
The averaging has been implemented and vastly improves the situation. See the manual section committed as part of #306.
from aspect.
Related Issues (20)
- random failure in tests/world_builder_select_grains HOT 2
- Modify Stress limiter in VP? HOT 1
- Update to particle composition HOT 5
- An error and its solution "bash: /home/ubuntu/aspect/aspect: No such file or directory"
- Deformed cell for 3D mesh deformation in release mode HOT 4
- Test viscoplastic material model
- Warning: handling deprecated packages HOT 4
- deal.II-v9.5.2 missing bundled directory HOT 8
- Move FastScape restarting functionality into the ASPECT checkpointing framework
- Fixed compositional field (and temperature) boundary conditions on deforming boundaries HOT 1
- Continental extension cookbook fails with gmg solver HOT 4
- contrib/utilities/update_scripts/source_files out of date
- MPI error installation main HOT 3
- Build ASPECT in 'release mode only' causes external plugins to not compile HOT 2
- Document default boundary condition for Stokes solve HOT 1
- Fix spherical unit vectors documentation
- benchmarks/free_surface_tractions/viscoelastic out of date
- Stripping the restart.mesh.data file HOT 5
- add documentation about different docker images
- A lingering interuption issue running on our local sever at UC Davis
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 aspect.