lhy11009

Okay, so that would be helpful indeed. However, I am sure that I am using the "block AMG" here all through.Here are the files. New ones: [log.txt](https://github.com/geodynamics/aspect/files/9620478/log.txt) [original.prm.txt](https://github.com/geodynamics/aspect/files/9620509/original.prm.txt) [statistics.txt](https://github.com/geodynamics/aspect/files/9620510/statistics.txt) Here...

@naliboff Sorry that I forgot to git add

Okay, thanks for your review, @naliboff . In this case, it makes sense to me to add test cases using the different compositions (the pyrolite, basalt, and harzburgite compositions). I...

@naliboff No pressure. I just like to hear from you about the tests to implement. I am free to meet, so if you feel more explanation from me is helpful,...

@bobmyhill, I agree with @bobmyhill , this is really neat solution.

@bobmyhill I just looked at the composite visco plastic material model and It seems to me the issue it has is different. It computes the composite viscosity directly by harmonic...

@MFraters you may also be interested in what we are doing here.

Hi Rene Thanks for your answer. > Repeated coarsening and refinement of a region, combined with addition/deletion of particles can cause this smearing This is also what I think the...

I am glad to be part of the conversation and I'll post what I find out from my 2-D & 3-D subduction models here.

I added `set Load balancing strategy = none` and `set Allow cells without particles = true` Here is the new result of the crustal composition at t =0.8 Ma ![eba_cdpt_SA80...