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;