Various minor tweaks.
Read in spontaneous polarization.
Add skeleton assembly of system matrix and system rhs.