Type 3 involves 2 vectors of positions, one in signal space, one in k space. Their respective derivatives should be related to type 1 and type 2. This may simply require evaluating type3 nufft a few times, but may require some understanding of the underlying procedure (there is likely a grid involved as an intermediate step and the size of this grid may be automatically determined)