compute_temp_boundary.cc
View Options

			#include "compute_temp_boundary.hh"
			/* -------------------------------------------------------------------------- */
			ComputeTemperatureBoundary::ComputeTemperatureBoundary(){}

			/* -------------------------------------------------------------------------- */

			void ComputeTemperatureBoundary::compute(System& system) {
			for (auto & par : system) {

			auto& p = static_cast<MaterialPoint&>(par);

			//enforcing zero boundary conditions in the physical space at only the left and bottom borders
			//(because of periodicity)
			if (p.getPosition()[0] == -1.0 \|\| p.getPosition()[1] == -1.0){
			p.getTemperature() = 0.0;
			}
			}
			}

			/* -------------------------------------------------------------------------- */

SP4E_HW1/homework3/starting_point/compute_temp_boundary.cccb94a3305560master