bugfix in Pmap for multiple s or zeta positions.

bugfix in Pmap (solve th+lamdba-th2=0),th2 was not correctly set if multiple s or zeta positions were evaluated. Extended the domain tests to check all combinations of 1d arrays of s,th,ze with multiple / single element