node = hou.pwd()
geo = node.geometry()
outgeo = node.geometry()

for x in range(0, 10):
    xformverb = hou.sopNodeTypeCategory().nodeVerb("xform")
    xformverb.setParms({"r": hou.Vector3(0, 10, 0)})
    xformverb.execute(outgeo, [geo])
    geo = outgeo

import hou

objA = hou.node("/obj/nodeA")
objB = hou.node("/obj/nodeB")

# get world space transforms of objects
xformA = objA.worldTransform()
xformB = objB.worldTransform()

# convert to quaternions
qA = hou.Quaternion(xformA)
qB = hou.Quaternion(xformB)
# create an output quaternion to store the result
qOut = hou.Quaternion(xformA)

# you can compare two quaternions using the dot product, just like vectors.
# if dot(q1,q2) is close to 1, they're almost identical.
# if they're not identical enough, slowly increase the bias when interpolating
# from A to B until we're finally basically at B.

threshold = 0.99
increment = 0.01
bias = 0.01
while qB.dot(qOut) < threshold:
    qOut = qA.slerp(qB, bias)
    bias += increment

# set output rotations based on this quaternion.
rOut = qOut.extractEulerRotates()
objA.parm("rx").set(rOut[0])
objA.parm("ry").set(rOut[1])
objA.parm("rz").set(rOut[2])

import hou

def match_xforms(sources, dests):
    with hou.undos.group("Matching transforms"):
        for x in range(0, len(sources)):
            source_obj = hou.node(sources[x])
            dest_obj = hou.node(dests[x])
            source_xform = source_obj.worldTransform()
            dest_xform = dest_obj.worldTransform()
            
            # get rotation from source, translate and scale from dest.
            source_rot = source_xform.extractRotationMatrix3()
            orig_t = dest_xform.extractTranslates()
            orig_s = dest_xform.extractScales()
            
            # create an output transform that's just the rotation.
            # by avoiding euler angles we don't have to worry about 
            # rotation order.
            out_xform = hou.Matrix4(source_rot)
            dest_obj.setWorldTransform(out_xform)
            
            # set translate and scale directly.
            dest_obj.parm("tx").set(orig_t[0])
            dest_obj.parm("ty").set(orig_t[1])
            dest_obj.parm("tz").set(orig_t[2])
            dest_obj.parm("sx").set(orig_s[0])
            dest_obj.parm("sy").set(orig_s[1])
            dest_obj.parm("sz").set(orig_s[2])
        
sources = ["/obj/sourceA", "/obj/sourceB", "/obj/sourceC"]
dests = ["/obj/destA", "/obj/destB", "/obj/destC"]
match_xforms(sources, dests)