Allow FieldMatrix and FieldVector to hold const SpatialFldPtr<FieldT>
In some cases, we may only have const SpatialFldPtr<FieldT> that we want to use to build the FieldMatrix and FieldVector objects, but we still want to perform operations like a linear solve. We currently don't support this.
@mahanse will need this very soon.