Numerical methods for black box software

A number of requirements are formulated to the numerical algorithms for black box software intended for mathematical modeling in continuum mechanics. An analysis of applied properties of the classical multigrid methods and robust multigrid technique in the framework of “robustness-efficiency-parallelism” problem is performed. It is shown that a close-to-optimal complexity with the least number of problem-dependent components and high parallel efficiency can be achieved with the robust multigrid technique on globally structured grids. Application of unstructured grids requires the accurate definition of two problem-dependent components (intergrid operators) that strongly affect on the complexity of an algorithm.