fixed adp density routine to not crash ... harmonic case tested with CE ... shape not checked