Obtaining the integral of a healpy map -

i have healpy map deflection (gradient of lensing potential), obtained during particular lensing simulation of cmb. want obtain map of lensing potential using healpy if possible. notice there healpy function alm2map_der1() give me healpy map , first derivative given map's alms. assuming first derivative gradient of map - please correct me if wrong. want know if can use healpy backwards process of this. want remove gradient , want lensing potential.

so far, attempt has been use relation between deflection , lensing potential power spectra; cls of deflection = l(l+1) * cls of lensing potential, , rearranging to: cls of lensing potential = cls of deflection / l(l+1), using synfast convert map. not seem getting correct map.

is there better way trying do? maybe not using healpy?

i can't first part, know converting cls destroys orientation information. 'synfast' provides map power spectrum input random orientation. if run 'synfast' on list of cls dipole power, random dipole orientation every time run:

import healpy healpy.mollview(healpy.synfast([0,1],32,lmax=1)) 

i suggest working alms if want map in end.