Discrete equivalence of adjoint neumann–dirichlet div-grad and grad-div equations in curvilinear 3d domains