Make all vectors non-ghosted except current_solution. Removed many compress(insert) since they are not needed. Solution is saved in directory output. Removed const qualifier on output_results since the Timer does not work.
Initial condition for Pk seems to work. Renamed some local variables. Formatted screen output.