# Concerning the charge densities I think you are doing it a bit backwards. Both the electric field and the charge density are scalar fields. I do not see the need for a class like ScalarPotential or ChargeDensity when for both of them you can use ScalarField. Moreover, I think it is just wrong to use ScalarPotential to represent a density of charge.
For instance something like this (if you want the superposition of an external field and a field created by some charge density): density = ScalarField(some_expression_defining_the_field_in_certain_coords) potential_internal = Laplacian(density) potential_external = ScalarFeild(another_expression) final_potential = potential_internal + potential_external .... stuff done with final_potential In CS terms, what I am saying is that you should not create new `types` and your module should not be strongly typed. It would be nice if some more CS-oriented people can comment on the above (Ronan possibly?) # Concerning degenerate distributions It would be great if we can switch from `ParametricRegion` to scalar field containing delta functions: for instance if the coordinates are r, phi in 2D a circle is both ParametricRegion(R*sin(phi), R*cos(phi)) and DiracDelta(r-R) -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy?hl=en-US. For more options, visit https://groups.google.com/groups/opt_out.
