void getCellCorners(NDArray& cellCornerValues, const QString &hostName, quint16 port, const qint32& caseId, const quint32& gridIndex) { QString serverName = hostName; quint16 serverPort = port; QTcpSocket socket; socket.connectToHost(serverName, serverPort); if (!socket.waitForConnected(riOctavePlugin::connectTimeOutMilliSecs)) { error((("Connection: ") + socket.errorString()).toLatin1().data()); return; } // Create command and send it: QString command = QString("GetCellCorners %1 %2").arg(caseId).arg(gridIndex); QByteArray cmdBytes = command.toLatin1(); QDataStream socketStream(&socket); socketStream.setVersion(riOctavePlugin::qtDataStreamVersion); socketStream << (qint64)(cmdBytes.size()); socket.write(cmdBytes); // Get response. First wait for the header while (socket.bytesAvailable() < (int)(5 * sizeof(quint64))) { if (!socket.waitForReadyRead(riOctavePlugin::longTimeOutMilliSecs)) { error((("Waiting for header: ") + socket.errorString()).toLatin1().data()); return; } } quint64 cellCountI; quint64 cellCountJ; quint64 cellCountK; quint64 cellCount; quint64 byteCount; socketStream >> cellCount; socketStream >> cellCountI; socketStream >> cellCountJ; socketStream >> cellCountK; socketStream >> byteCount; if (!(byteCount && cellCount)) { error ("Could not find the requested data in ResInsight"); return; } dim_vector dv; dv.resize(5); dv(0) = cellCountI; dv(1) = cellCountJ; dv(2) = cellCountK; dv(3) = 8; dv(4) = 3; cellCornerValues.resize(dv); double* internalMatrixData = cellCornerValues.fortran_vec(); QStringList errorMessages; if (!RiaSocketDataTransfer::readBlockDataFromSocket(&socket, (char*)(internalMatrixData), byteCount, errorMessages)) { for (int i = 0; i < errorMessages.size(); i++) { error(errorMessages[i].toLatin1().data()); } OCTAVE_QUIT; } return; }
void getCellCorners(NDArray& cellCornerValues, const QString &hostName, quint16 port, const qint32& caseId, const quint32& gridIndex) { QString serverName = hostName; quint16 serverPort = port; QTcpSocket socket; socket.connectToHost(serverName, serverPort); if (!socket.waitForConnected(riOctavePlugin::connectTimeOutMilliSecs)) { error((("Connection: ") + socket.errorString()).toLatin1().data()); return; } // Create command and send it: QString command = QString("GetCellCorners %1 %2").arg(caseId).arg(gridIndex); QByteArray cmdBytes = command.toLatin1(); QDataStream socketStream(&socket); socketStream.setVersion(riOctavePlugin::qtDataStreamVersion); socketStream << (qint64)(cmdBytes.size()); socket.write(cmdBytes); // Get response. First wait for the header while (socket.bytesAvailable() < (int)(5 * sizeof(quint64))) { if (!socket.waitForReadyRead(riOctavePlugin::shortTimeOutMilliSecs)) { error((("Waiting for header: ") + socket.errorString()).toLatin1().data()); return; } } quint64 cellCountI; quint64 cellCountJ; quint64 cellCountK; quint64 cellCount; quint64 byteCount; socketStream >> cellCount; socketStream >> cellCountI; socketStream >> cellCountJ; socketStream >> cellCountK; socketStream >> byteCount; if (!(byteCount && cellCount)) { error ("Could not find the requested data in ResInsight"); return; } dim_vector dv; dv.resize(5); dv(0) = cellCountI; dv(1) = cellCountJ; dv(2) = cellCountK; dv(3) = 8; dv(4) = 3; cellCornerValues.resize(dv); while (socket.bytesAvailable() < (qint64)(byteCount)) { if (!socket.waitForReadyRead(riOctavePlugin::longTimeOutMilliSecs)) { error((("Waiting for data: ") + socket.errorString()).toLatin1().data()); return; } OCTAVE_QUIT; } double* internalMatrixData = cellCornerValues.fortran_vec(); #if 0 double val; for (octave_idx_type i = 0; i < valueCount; i++) { socketStream >> internalMatrixData[i]; } #else quint64 bytesRead = 0; bytesRead = socket.read((char*)(internalMatrixData), byteCount); if (byteCount != bytesRead) { error("Could not read binary double data properly from socket"); octave_stdout << "Cell count: " << cellCount << std::endl; } #endif return; }
void getActiveCellCenters(NDArray& cellCenterValues, const QString &hostName, quint16 port, const qint32& caseId, const QString& porosityModel) { QString serverName = hostName; quint16 serverPort = port; QTcpSocket socket; socket.connectToHost(serverName, serverPort); if (!socket.waitForConnected(riOctavePlugin::connectTimeOutMilliSecs)) { error((("Connection: ") + socket.errorString()).toLatin1().data()); return; } // Create command and send it: QString command = QString("GetActiveCellCenters %1 %2").arg(caseId).arg(porosityModel); QByteArray cmdBytes = command.toLatin1(); QDataStream socketStream(&socket); socketStream.setVersion(riOctavePlugin::qtDataStreamVersion); socketStream << (qint64)(cmdBytes.size()); socket.write(cmdBytes); // Get response. First wait for the header while (socket.bytesAvailable() < (int)(2 * sizeof(quint64))) { if (!socket.waitForReadyRead(riOctavePlugin::shortTimeOutMilliSecs)) { error((("Waiting for header: ") + socket.errorString()).toLatin1().data()); return; } } // Read timestep count and blocksize quint64 activeCellCount; quint64 byteCount; socketStream >> activeCellCount; socketStream >> byteCount; if (!(byteCount && activeCellCount)) { error ("Could not find the requested data in ResInsight"); return; } dim_vector dv; dv.resize(2); dv(0) = activeCellCount; dv(1) = 3; cellCenterValues.resize(dv); while (socket.bytesAvailable() < (qint64)(byteCount)) { if (!socket.waitForReadyRead(riOctavePlugin::longTimeOutMilliSecs)) { error((("Waiting for data: ") + socket.errorString()).toLatin1().data()); return; } OCTAVE_QUIT; } quint64 bytesRead = 0; double* internalMatrixData = cellCenterValues.fortran_vec(); bytesRead = socket.read((char*)(internalMatrixData), byteCount); if (byteCount != bytesRead) { error("Could not read binary double data properly from socket"); octave_stdout << "Active cell count: " << activeCellCount << std::endl; } return; }
void getGridProperty(NDArray& propertyFrames, const QString &serverName, quint16 serverPort, const int& caseId, int gridIdx, QString propertyName, const int32NDArray& requestedTimeSteps, QString porosityModel) { QTcpSocket socket; socket.connectToHost(serverName, serverPort); if (!socket.waitForConnected(riOctavePlugin::connectTimeOutMilliSecs)) { error((("Connection: ") + socket.errorString()).toLatin1().data()); return; } QDataStream socketStream(&socket); socketStream.setVersion(riOctavePlugin::qtDataStreamVersion); // Create command as a string with arguments , and send it: QString command; command += "GetGridProperty " + QString::number(caseId) + " " + QString::number(gridIdx) + " " + propertyName + " " + porosityModel; for (int i = 0; i < requestedTimeSteps.length(); ++i) { if (i == 0) command += " "; command += QString::number(static_cast<int>(requestedTimeSteps.elem(i)) - 1); // To make the index 0-based if (i != requestedTimeSteps.length() -1) command += " "; } QByteArray cmdBytes = command.toLatin1(); socketStream << (qint64)(cmdBytes.size()); socket.write(cmdBytes); // Get response. First wait for the header while (socket.bytesAvailable() < (int)(4*sizeof(quint64))) { if (!socket.waitForReadyRead(riOctavePlugin::longTimeOutMilliSecs)) { error((("Waiting for header: ") + socket.errorString()).toLatin1().data()); return; } } // Read sizes quint64 totalByteCount; quint64 cellCountI; quint64 cellCountJ; quint64 cellCountK; quint64 timestepCount; socketStream >> cellCountI; socketStream >> cellCountJ; socketStream >> cellCountK; socketStream >> timestepCount; totalByteCount = cellCountI*cellCountJ*cellCountK*timestepCount*sizeof(double); if (!(totalByteCount)) { error ("Could not find the requested data in ResInsight"); return; } dim_vector dv; dv.resize(4); dv(0) = cellCountI; dv(1) = cellCountJ; dv(2) = cellCountK; dv(3) = timestepCount; propertyFrames.resize(dv); // Wait for available data while (socket.bytesAvailable() < (int)totalByteCount) { if (!socket.waitForReadyRead(riOctavePlugin::longTimeOutMilliSecs)) { error(("Waiting for data : " + socket.errorString()).toLatin1().data()); return ; } OCTAVE_QUIT; } qint64 bytesRead = 0; double * internalMatrixData = propertyFrames.fortran_vec(); // Raw data transfer. Faster. bytesRead = socket.read((char*)(internalMatrixData ), totalByteCount); if ((int)totalByteCount != bytesRead) { error("Could not read binary double data properly from socket"); } QString tmp = QString("riGetGridProperty : Read %1").arg(propertyName); if (caseId < 0) { tmp += QString(" from current case,"); } else { tmp += QString(" from case with Id: %1,").arg(caseId); } tmp += QString(" grid index: %1, ").arg(gridIdx); octave_stdout << tmp.toStdString() << " I, J, K " << cellCountI << ", " << cellCountJ << ", " << cellCountK << ", Timesteps : " << timestepCount << std::endl; return; }