A parallel technique for an adaptive solution of 3D boundary value problems is described. It incorporates a parallel mesh generation and a parallel iterative solution of the corresponding discrete problem. Both generation and solution are problem independent and may be considered as black-boxes.