Unstructured block grids
Function1<double> Sin(sin); // definition of functions Function1<double> Sinh(sinh); X_coordinate X(grid); Y_coordinate Y(grid); Z_coordinate Z(grid); Local_stiffness_matrix<double> Laplace(grid); Local_stiffness_matrix<double> Helm(grid); Laplace.Calculate(grad(v_())*grad(w_())); Helm.Calculate(v_()*w_()); u_exakt = Sin(X)*Sinh(Y)*Z; // exact solution u = 0.0;