mass_weight_matrix=mass_weight_frame(mass_matrix,gaussian_matrix)

function mass_weight_matrix=mass_weight_frame(mass_matrix,gaussian_matrix)

N=length(mass_matrix);

reshape_gaussian_matrix=reshape(gaussian_matrix,[N,3*length(gaussian_matrix(1,1,:))]);

reshape_mass_weight_gaussian_matrix=diag(sqrt(mass_matrix))*reshape_gaussian_matrix;

mass_weight_matrix=reshape(reshape_mass_weight_gaussian_matrix,[N,3,length(gaussian_matrix(1,1,:))]);