import bpy

obj = bpy.context.active_object
me = obj.data

norms = [None] * len(me.loops)

me.calc_tangents()

for i in range(len(me.loops)):
    norms[i] = me.loops[i].normal

me.normals_split_custom_set(norms)