unutbu did great. I suggest one equivalent without where (but with numpy anyway)
import numpy as np s1=[2,1,2,5,4,6] s2=[1,2,4,5,7,8] s3=[0.1,0.4,0.5,0.6,0.1,0.1] res = [xv if c else yv for (c,xv,yv) in zip([si1<si2 for si1,si2 in zip(s1,s2)], list(np.arcsin(s3)), list(np.arccos(s3)))]
If you print zip() you get this list
>>> [(False, 0.1001674211615598, 1.4706289056333368), (True, 0.41151684606748806, 1.1592794807274085), (True, 0.52359877559829893, 1.0471975511965979), (False, 0.64350110879328437, 0.9272952180016123), (True, 0.1001674211615598, 1.4706289056333368), (True, 0.1001674211615598, 1.4706289056333368)]
Take the first element (False, 0.1001674211615598, 1.4706289056333368) : 2<1 really false. So you will have 1.4706289056333368 as the first value in res .
Result
>>> res [1.4706289056333368, 0.41151684606748806, 0.52359877559829893, 0.9272952180016123, 0.1001674211615598, 0.1001674211615598]