OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fault Reactivation: parse element sets from INP file.

Browse Source
This commit is contained in:
Kristian Bendiksen 2023-11-06 14:34:53 +01:00 committed by jonjenssen
parent a61e62707b
commit 065e37a437
2 changed files with 56 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
ApplicationLibCode/GeoMech/OdbReader/RifInpReader.cpp
View File

@ -84,8 +84,9 @@ bool RifInpReader::readFemParts( RigFemPartCollection* femParts )
std::map<int, std::string> parts;
std::map<int, std::vector<std::pair<int, cvf::Vec3d>>> nodes;
std::map<int, std::vector<std::pair<int, std::vector<int>>>> elements;
std::map<int, std::vector<std::pair<std::string, std::vector<size_t>>>> elementSets;
read( m_stream, parts, nodes, elements );
read( m_stream, parts, nodes, elements, elementSets );
RiaLogging::debug( QString( "Read FEM parts: %1 nodes: %2 elements: %3" ).arg( parts.size() ).arg( nodes.size() ).arg( elements.size() ) );
@ -147,6 +148,13 @@ bool RifInpReader::readFemParts( RigFemPartCollection* femParts )
femPart->appendElement( elmType, elmId, indexBasedConnectivities.data() );
}
// read element sets
auto elementSetsForPart = elementSets[partId];
for ( auto [setName, elementSet] : elementSetsForPart )
{
femPart->addElementSet( setName, elementSet );
}
femPart->setElementPartId( femParts->partCount() );
femParts->addFemPart( femPart );
@ -162,7 +170,8 @@ bool RifInpReader::readFemParts( RigFemPartCollection* femParts )
void RifInpReader::read( std::istream& stream,
std::map<int, std::string>& parts,
std::map<int, std::vector<std::pair<int, cvf::Vec3d>>>& nodes,
std::map<int, std::vector<std::pair<int, std::vector<int>>>>& elements )
std::map<int, std::vector<std::pair<int, std::vector<int>>>>& elements,
std::map<int, std::vector<std::pair<std::string, std::vector<size_t>>>>& elementSets )
{
std::string line;
@ -202,6 +211,13 @@ void RifInpReader::read( std::istream&
skipComments( stream );
elements[partId] = readElements( stream );
}
else if ( uppercasedLine.starts_with( "*ELSET," ) )
{
skipComments( stream );
std::string setName = parseLabel( line, "elset" );
auto elementSet = readElementSet( stream );
elementSets[partId].push_back( { setName, elementSet } );
}
continue;
}
@ -284,6 +300,32 @@ std::vector<std::pair<int, std::vector<int>>> RifInpReader::readElements( std::i
return partElements;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<size_t> RifInpReader::readElementSet( std::istream& stream )
{
std::vector<size_t> elementSet;
// Read until we find a new section (which should start with a '*').
while ( stream.peek() != '*' && stream.peek() != EOF )
{
// Read entire line of comma-separated values
std::string line;
std::getline( stream, line );
// Process the comma-separated values
auto parts = RiaStdStringTools::splitString( line, ',' );
for ( auto part : parts )
{
int elementId = RiaStdStringTools::toInt( part ) - 1;
elementSet.push_back( elementId );
}
}
return elementSet;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------

4
ApplicationLibCode/GeoMech/OdbReader/RifInpReader.h
View File

@ -77,11 +77,13 @@ private:
static std::string parseLabel( const std::string& line, const std::string& labelName );
static std::vector<std::pair<int, cvf::Vec3d>> readNodes( std::istream& stream );
static std::vector<std::pair<int, std::vector<int>>> readElements( std::istream& stream );
static std::vector<size_t> readElementSet( std::istream& stream );
static void read( std::istream& stream,
std::map<int, std::string>& parts,
std::map<int, std::vector<std::pair<int, cvf::Vec3d>>>& nodes,
std::map<int, std::vector<std::pair<int, std::vector<int>>>>& elements );
std::map<int, std::vector<std::pair<int, std::vector<int>>>>& elements,
std::map<int, std::vector<std::pair<std::string, std::vector<size_t>>>>& elementSets );
private:
std::map<int, std::vector<std::string>> m_partElementSetNames;