[Zhen ZHANG]

Hi all,

I got a problem when setting values at the boundary to the volume data container.

I created a VectorType member sens_inlet in the CEulerSolver class and resized it to nPoint.

I then registered the Tt of all vertices at the inlet boundary. After extracting the sensitivity, I set them to the sens_inlet using the following codes:

for(unsigned short iMarker=0; iMarker < config->GetnMarker_All(); iMarker++) {
  for (unsigned long iVertex = 0; iVertex < geometry->GetnVertex(iMarker); iVertex++) {
    unsigned long iPoint = geometry->nodes->GetGlobalIndex(geometry->vertex[iMarker][iVertex]->GetNode());
    nodes->SetSens_Inlet_T(iPoint,SU2_TYPE::GetDerivative(Inlet_Ttotal[iMarker][iVertex]))}}

The orders of registration and extraction are the same.

However, in the result file, the non-zero values of the sens_inlet do not appear at the inlet boundary (see the attached figure). If running with MPI, it returns an error "Assertion `i>=0 && i<m_rows' failed".

Could anybody give me some clues why this problem happen? Thanks in advance!

[Zhen ZHANG]

Well, I changed the iPoint calculation to

unsigned long iPoint = geometry->vertex[iMarker][iVertex]->GetNode();

and it works.
What is the difference between these two calculations?

Thanks,
Zhen