1453 |
const MeshAdapter& newDomain=dynamic_cast<const MeshAdapter&>(*(new_x.getFunctionSpace().getDomain())); |
const MeshAdapter& newDomain=dynamic_cast<const MeshAdapter&>(*(new_x.getFunctionSpace().getDomain())); |
1454 |
if (newDomain!=*this) |
if (newDomain!=*this) |
1455 |
throw FinleyAdapterException("Error - Illegal domain of new point locations"); |
throw FinleyAdapterException("Error - Illegal domain of new point locations"); |
1456 |
tmp = new_x.getDataC(); |
if ( new_x.getFunctionSpace() == continuousFunction(asAbstractContinuousDomain()) ) { |
1457 |
Finley_Mesh_setCoordinates(mesh,&tmp); |
tmp = new_x.getDataC(); |
1458 |
|
Finley_Mesh_setCoordinates(mesh,&tmp); |
1459 |
|
} else { |
1460 |
|
escript::Data new_x_inter=escript::Data( new_x, continuousFunction(asAbstractContinuousDomain()) ); |
1461 |
|
tmp = new_x_inter.getDataC(); |
1462 |
|
Finley_Mesh_setCoordinates(mesh,&tmp); |
1463 |
|
} |
1464 |
checkFinleyError(); |
checkFinleyError(); |
1465 |
} |
} |
1466 |
|
|