Optimizing wave energy parks with over 1000 interacting point-absorbers using an approximate analytical method