Commit e89b40e6 authored by Ron Rahaman's avatar Ron Rahaman

Transfer to Nek working

parent af2c7951
......@@ -52,6 +52,14 @@ namespace Nek5000
double flux_moose, temp_nek;
} test_passing_;
extern struct {
double coeff_tij[];
} expansion_tcoef_;
extern struct {
double coeff_fij[];
} expansion_fcoef_;
// subroutine nek_init(intracomm)
void FORTRAN_CALL(nek_init)(const int&); // Ron likes this better now
void FORTRAN_CALL(nek_init_step)();
......
......@@ -9,6 +9,8 @@
#include "libmesh/meshfree_interpolation.h"
#include "libmesh/system.h"
#include "NekInterface.h"
// Define the input parameters
template<>
InputParameters validParams<MultiAppPolynomialToNek>()
......@@ -56,8 +58,11 @@ MultiAppPolynomialToNek::execute()
_console << _source_variable_names[i] << '\n';
auto & solution_values = source_variables[i]->sln();
for (auto j = beginIndex(solution_values); j < solution_values.size(); ++j)
for (auto j = beginIndex(solution_values); j < solution_values.size(); ++j)
{
_console << solution_values[j] << ' ';
Nek5000::expansion_fcoef_.coeff_fij[i*100+j] = solution_values[j];
}
_console << '\n';
}
}
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment