Stack frames: 10 0: print_trace(std::ostream&) 1: FEBase::compute_single_point_map(std::vector<double, std::allocator<double> > const&, Elem const*, unsigned int) 2: FEBase::compute_map(std::vector<double, std::allocator<double> > const&, Elem const*) 3: FE<3u, (libMeshEnums::FEFamily)0>::reinit(Elem const*, std::vector<Point, std::allocator<Point> > const*) 4: compute_residual(NumericVector<double> const&, NumericVector<double>&) 5: __libmesh_petsc_snes_residual 6: SNESComputeFunction(_p_SNES*, _p_Vec*, _p_Vec*) 7: SNESSolve_LS(_p_SNES*) 8: SNESSolve(_p_SNES*, _p_Vec*, _p_Vec*) 9: PetscNonlinearSolver<double>::solve(SparseMatrix<double>&, NumericVector<double>&, NumericVector<double>&, double, unsigned int)