Handle DomainErrors

Also tweak domains for some tests
Add benchmarks for A.^b where b is scalar