Improve if conditionals and set a default case to evaluate viscosity in calc_visc

ggciag / mandyoc

MANDYOC is a finite element code written on top of the PETSc library to simulate thermo-chemical convection of the Earth's mantle

https://ggciag.github.io/mandyoc/

BSD 3-Clause "New" or "Revised" License

27 stars 5 forks source link

Open rafaelmds opened 2 years ago

rafaelmds commented 2 years ago

Currently, in calc_visc.cpp there is a series of ifs conditionals to match a rheology and then to evaluate the viscosity.

This need to be improved. A first improvement could be using a switch-case pattern with a default case.