Make the pullback for split_network more general

[mtfishman]

Right now the pullback for split_network assumes that the network was split by priming one index of a pair of link indices, and also that none of the tensors of the original network before splitting were primed. It would be good to generalize this assumption.

This could be done by analyzing which pairs of indices get split in the forward pass, and how they were split, and use that information in the reverse pass.

ITensorNetworkAD.jl/src/Optimizations/itensor_network.jl

Lines 5 to 15 in 00d2811

	function ChainRulesCore.rrule(
	::typeof(split_network),
	tn::Matrix{ITensor};
	projector_center=default_projector_center(tn),
	)
	function pullback(dtn_split::Matrix{ITensor})
	dtn = map(t -> replaceprime(t, 1 => 0), dtn_split)
	return (NoTangent(), dtn, NoTangent())
	end
	return split_network(tn; projector_center=projector_center), pullback
	end