Error happend when writing Velocites from CalculiX

Official adapters and tutorials
,
1

Hi, everyone

Recently, I tried to writing velocities from CalculiX to another solver. I carried 2D simulation, the Segmentation fault appear after the adapter writing coupling data, as follow:

图片
图片1159×276 67.4 KB

However, this error won’t appear in 3D simulation, and this situation is the same with this link :

So I want to know if this problem has been solved.

Thank you!

2

I find out that this may be due to the incomplete coding of writing quasi 2D3D data of velocities.
In PreciceInterface.c

case DISPLACEMENTDELTAS:
        if (isQuasi2D3D(interfaces[i]->quasi2D3D)) {
          getNodeDisplacementDeltas(interfaces[i]->nodeIDs, interfaces[i]->numNodes, interfaces[i]->dimCCX, sim->vold, sim->coupling_init_v, sim->mt, interfaces[i]->mappingQuasi2D3D->bufferVector3D);
          consistentVectorWrite(interfaces[i]->mappingQuasi2D3D, interfaces[i]->couplingMeshName, interfaces[i]->displacementDeltas);
        } else {
          getNodeDisplacementDeltas(interfaces[i]->nodeIDs, interfaces[i]->numNodes, interfaces[i]->dimCCX, sim->vold, sim->coupling_init_v, sim->mt, interfaces[i]->nodeVectorData);
          precicec_writeData(interfaces[i]->couplingMeshName, interfaces[i]->displacementDeltas, interfaces[i]->numNodes, interfaces[i]->preciceNodeIDs, interfaces[i]->nodeVectorData);
        }
        printf("Writing DISPLACEMENTDELTAS coupling data.\n");
        break;
case VELOCITIES:
        getNodeVelocities(interfaces[i]->nodeIDs, interfaces[i]->numNodes, interfaces[i]->dim, sim->veold, sim->mt, interfaces[i]->nodeVectorData);
        precicec_writeData(interfaces[i]->couplingMeshName, interfaces[i]->velocities, interfaces[i]->numNodes, interfaces[i]->preciceNodeIDs, interfaces[i]->nodeVectorData);
        printf("Writing VELOCITIES coupling data.\n");
        break;
3

I add some codes in PreciceInterface.c. I am not a expert of calculix adapter, so I’m not sure if I coding it right. From the simulation now, it seems ok. If anyone is familiar with this, please check it.

case VELOCITIES:
        if (isQuasi2D3D(interfaces[i]->quasi2D3D)) {
          getNodeVelocities(interfaces[i]->nodeIDs, interfaces[i]->numNodes, interfaces[i]->dimCCX, sim->veold, sim->mt, interfaces[i]->mappingQuasi2D3D->bufferVector3D);
          consistentVectorWrite(interfaces[i]->mappingQuasi2D3D, interfaces[i]->couplingMeshName, interfaces[i]->velocities);
        } else {
          getNodeVelocities(interfaces[i]->nodeIDs, interfaces[i]->numNodes, interfaces[i]->dim, sim->veold, sim->mt, interfaces[i]->nodeVectorData);
          precicec_writeData(interfaces[i]->couplingMeshName, interfaces[i]->velocities, interfaces[i]->numNodes, interfaces[i]->preciceNodeIDs, interfaces[i]->nodeVectorData);
        }
        printf("Writing VELOCITIES coupling data.\n");
        break;

Thank you~