it gives all the points of a Coords a different rotation.
per every point , you can specifify a differnt rotation axis, a differnt angle and initial point around which to rotate.

I put the rotation matrix inside the function and did all what the normal rotate was doing by calling different methods because

arraytools.rotationMatrix cretae only a 3x3 matrix (now is 3x3x ncoords)
Coords.rotate make a check if rotationMatrix is 3x3 , error otherwise
affine uses the dot product which does not allow to perform multiple rotation I want (every rotmatrix slice X evry coord)

I think by modyfieng the mentioned funtions we can just extend rotate to perform multiple rotations

I created this function where
angle,axis,around
can have values of the coords array. I think rotate affine and other funtions in coords can be all modified to be extended to this approach so that the function can use already existing rotation fuunctions. I am afraid to modify the core so i post the code to get suggestions abbout how it can be embedded in coords.

def rotmany(coords,angle,axis,around=None):
    out=coords
    if around is not None:
        around = asarray(around)
        out = out.translate(-around)

    angle=asarray(angle,float)
    axis=asarray(axis).reshape(-1,3)

    a = angle # in degree
    c = cosd(a)
    s = sind(a)
    t = 1.-c

    X,Y,Z = axis.T
    rot = [ [ t*X*X + c  , t*X*Y + s*Z, t*X*Z - s*Y ],
              [ t*Y*X - s*Z, t*Y*Y + c  , t*Y*Z + s*X ],
              [ t*Z*X + s*Y, t*Z*Y - s*X, t*Z*Z + c   ] ]
    rot=asarray(rot)

    out=out.T.reshape(3,1,-1)
    out=(rot*out).sum(axis=0).T

    return Coords(out)