// protected override Band_Data Get_ChemPot_From_External(Band_Data density)
// {
// Save_Data(density, dens_filename);
//
// return Get_Data_From_External(flexpde_script, initcalc_result_filename);
// }
/// <summary>
/// Calculates the Laplacian for the given band structure of the input potential
/// ie. returns d(eps * d(input))
/// NOTE: the input should be a potential so make sure you divide all band energies by q_e
/// </summary>
public Band_Data Calculate_Phi_Laplacian(Band_Data input)
{
DoubleTriDiagMatrix lap_mat = Generate_Laplacian(exp.Layers);
return(new Band_Data(MatrixFunctions.Product(lap_mat, input.vec)));
}
