compute_temperature.cc
79ad179388bdmaster

compute_temperature.cc
View Options

			#include "compute_temperature.hh"
			#include "fft.hh"
			#include "material_point.hh"
			#include <cmath>

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

			void ComputeTemperature::compute(System& system) {
			// create temperature field matrix --> skech of operations

			// casting the matrices
			// preparing the temp and heatrate vectors
			std::vector<Real> temperatures;
			std::vector<Real> heat_rates;
			// getting the vectors
			for (auto& npar : system){
			auto& part = static_cast<MaterialPoint&>(npar);
			temperatures.push_back(part.getTemperature());
			heat_rates.push_back(part.getHeatRate());
			}
			// creating the matrices out of the vectors
			// UInt nb_particles = system.getNbParticles();
			// UInt size = std::sqrt(nb_particles);


			auto temperatures_matrix = makeMatrix(temperatures);
			auto hv_matrix = makeMatrix(heat_rates);

			// step 1: forward fft
			auto temperaturesfft_matrix = FFT::transform(temperatures_matrix);
			auto hvfft_matrix = FFT::transform(hv_matrix);
			/* auto qx == FFT::computeFrequencies()...
			auto qy == FFT::transform()...
			// step 2: dteta^/dt = ...
			qx2 = ...
			qy2 = ...
			auto kTe1 = mat_mult_sm(k, temperaturesfft_matrix);
			auto dTe_dt_fft = matrix_sub_mm(hvfft_matrix ,mat_mult_mm(kTe1 ,matrix_add_mm(qx2, qy2)));

			// step 3: backward fft
			auto dTe_dt = FFT::itransform(dTe_dt_fft)

			// step 4: temperatures_matrix+1
			auto tempertures_matrix_1 = matrix_sum_mm(temperatures_matrix, matrix_mult_sm(dt, dTe_dt))

			*/
			}

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