Joining curves removes non active object's materials.
I fixed it by porting code for mesh join operator (ED_mesh_join_objects_exec and join_mesh_single in meshtools.c) to curve.

I used the blend file uploaded in T88138