Speeding up independent eigenvalue-calculations for many parametersets and/or for many different points

Question

I hope that someone could help me with my problem that is about solving eigenvalue problems for a matrix h_matrix(x, y, param1, param2) that depends on 2 variables