From 0d267b1190281864f48a983c2c708aa9cfe8452e Mon Sep 17 00:00:00 2001
From: Ankang Gao <ankanggao@ustc.edu.cn>
Date: Mon, 1 Jul 2024 15:03:32 +0800
Subject: [PATCH] update
---
library/FieldUtils/ProcessModules/ProcessMRF.cpp | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/library/FieldUtils/ProcessModules/ProcessMRF.cpp b/library/FieldUtils/ProcessModules/ProcessMRF.cpp
index 2971562a62..49425dd4bd 100644
--- a/library/FieldUtils/ProcessModules/ProcessMRF.cpp
+++ b/library/FieldUtils/ProcessModules/ProcessMRF.cpp
@@ -58,6 +58,7 @@ ProcessMRF::~ProcessMRF()
void ProcessMRF::v_Process(po::variables_map &vm)
{
m_f->SetUpExp(vm);
+ ReadMRFData();
// Skip in case of empty partition
if (m_f->m_exp[0]->GetNumElmts() == 0)
@@ -67,7 +68,6 @@ void ProcessMRF::v_Process(po::variables_map &vm)
// Determine dimensions
m_spacedim = m_f->m_graph->GetMeshDimension() + m_f->m_numHomogeneousDir;
- ReadMRFData();
// transform coordinates
int nfields = m_f->m_variables.size();
@@ -192,6 +192,7 @@ void ProcessMRF::TransformVector(std::vector<Array<OneD, NekDouble>> &data)
for (int i = 0; i < dim; ++i)
{
tmp[i] = Array<OneD, NekDouble>(npoint);
+ Vmath::Vcopy(npoint, data[i], 1, tmp[i], 1);
}
Vmath::Svtsvtp(npoint, cz, tmp[0], 1, -sz, tmp[1], 1, data[0], 1);
Vmath::Svtsvtp(npoint, sz, tmp[0], 1, cz, tmp[1], 1, data[1], 1);
--
GitLab