We present a grid-based non-parametric approach to obtain a triaxial three-dimensional luminosity density from its surface brightness distribution. Triaxial deprojection is highly degenerate and our approach illustrates the resulting difficulties. Fortunately, for massive elliptical galaxies, many deprojections for a particular line of sight can be discarded, because their projection along other lines of sight does not resemble elliptical galaxies. The near-elliptical isophotes of these objects imply near ellipsoidal intrinsic shapes. In fact, deprojection is unique for densities distributed on ellipsoidal shells. The constrained non-parametric deprojection method we present here relaxes this constraint and assumes that the contours of the luminosity density are boxy/discy ellipsoids with radially varying axis ratios. With this approach we are able to reconstruct the intrinsic triaxial densities of our test models, including one drawn from an $N$-body simulation. The method also allows to compare the relative likelihood of deprojections at different viewing angles. We show that the viewing orientations of individual galaxies with nearly ellipsoidal isophotal shapes can be constrained from photometric data alone.

We study the correlations between Supermassive Black Holes (BH) and their host galaxies, using a sample of 83 BH masses collected from the most recent and reliable spatially resolved estimates available from the literature. We confirm the mono- and bivariate correlations between SMBHs and the bulges of their host galaxies, confirming that the correlation with the effective velocity dispersion is not significantly improved by higher dimensionality. Instead, pseudobulges do not seem to correlate with their SMBHs, probably because their secular evolution is often unable to trigger accretion onto the central BH. We then present a novel approach aimed at finding the fundamental relation between SMBHs and their host galaxies. For the first time, we analytically combine BH masses with the Fundamental Plane (FP), showing that Mbh-sigma_e appears to be the fundamental relation rather than a putative "BH Fundamental Plane" of the kind Mbh-sigma_e-R_e. These results can be explained by a picture which sees the Mbh-sigma_e relation as a natural outcome of the change in AGN feedback from momentum- to energy-driven. The other scaling relations are then established through the FP.