3D OGRSpatialReference::Transform is failed with PROJ 9.3 (regression!, it works with PROJ 9.2)

[aharondavid]

GDAL Version: 3.1.2
PROJ Version: 9.3
OS: Windows 11

Code to reproduced:

std::string SourceWKT ="GEOGCS["WGS84 Coordinate System",DATUM["WGS_1984",SPHEROID["WGS 1984",6378137,298.257223563],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]";

std::string DestWKT ="GEOCCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Geocentric X",OTHER],AXIS["Geocentric Y",OTHER],AXIS["Geocentric Z",NORTH],AUTHORITY["EPSG","4978"]]";

OGRSpatialReference sourceSR;
sourceSR.importFromWkt(SourceWKT.c_str());
sourceSR.PromoteTo3D(NULL);

OGRSpatialReference destSR;
destSR.importFromWkt(DestWKT.c_str());
destSR.PromoteTo3D(NULL);



OGRCoordinateTransformation*  ogrTransformation = OGRCreateCoordinateTransformation(sourceSR,destSR);

double X	=-0.31036012488063702;
double Y=	39.611038380093397;
double Z	=28.253883366008200;

int ret = ogrTransformation->Transform(X,Y,Z);
//Result:
//ret equal to 0
//X = inf
//Y = inf
//Z = inf

[rouault]

Transfering as a PROJ issue, and I don't confirm it is a regression. Using projinfo, I can reproduce the issue at least from PROJ 8.0.0 (didn't try older ones)