zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge branch 'ensight-timeFormat' into 'develop'

Browse Source 
support user specification of ensight time format/precision (#2999)

See merge request Development/openfoam!634
This commit is contained in:
Andrew Heather
2023-10-19 10:39:07 +00:00
parent 8ec35d0b26 ba2e774aac
commit 2b34412591
28 changed files with 546 additions and 290 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/OpenFOAM/db/Time/Time.C
View File

@ -73,7 +73,10 @@ Foam::Time::writeControlNames
});
Foam::Time::fmtflags Foam::Time::format_(Foam::Time::fmtflags::general);
Foam::IOstreamOption::floatFormat Foam::Time::format_
(
IOstreamOption::floatFormat::general
);
int Foam::Time::precision_(6);

14
src/OpenFOAM/db/Time/Time.H
View File

@ -104,15 +104,6 @@ public:
saUnknown //!< Dummy no-op. Do not change current value.
};
//- Supported time directory name formats
enum fmtflags
{
general = 0, //!< default float notation
fixed = ios_base::fixed, //!< fixed-point notation
scientific = ios_base::scientific //!< scientific notation
};
//- Names for writeControls
static const Enum<writeControls> writeControlNames;
@ -179,9 +170,8 @@ protected:
//- Signal handler for write and clean exit upon signal
sigStopAtWriteNow sigStopAtWriteNow_;
//- Time directory name format
static fmtflags format_;
//- Format for time directory names (general | fixed | scientific)
static IOstreamOption::floatFormat format_;
//- Time directory name precision
static int precision_;

25
src/OpenFOAM/db/Time/TimeIO.C
View File

@ -329,29 +329,8 @@ void Foam::Time::readDict()
}
}
if (controlDict_.found("timeFormat"))
{
const word formatName(controlDict_.get<word>("timeFormat"));
if (formatName == "general")
{
format_ = fmtflags::general;
}
else if (formatName == "fixed")
{
format_ = fmtflags::fixed;
}
else if (formatName == "scientific")
{
format_ = fmtflags::scientific;
}
else
{
WarningInFunction
<< "Unsupported time format " << formatName
<< endl;
}
}
format_ =
IOstreamOption::floatFormatEnum("timeFormat", controlDict_, format_);
controlDict_.readIfPresent("timePrecision", precision_);

83
src/OpenFOAM/db/options/IOstreamOption.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2018-2020 OpenCFD Ltd.
Copyright (C) 2018-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -35,6 +35,17 @@ License
const Foam::IOstreamOption::versionNumber Foam::IOstreamOption::currentVersion;
const Foam::Enum
<
Foam::IOstreamOption::floatFormat
>
Foam::IOstreamOption::floatFormatNames
({
{ floatFormat::general, "general" },
{ floatFormat::fixed, "fixed" },
{ floatFormat::scientific, "scientific" },
});
const Foam::Enum
<
Foam::IOstreamOption::streamFormat
@ -48,27 +59,72 @@ Foam::IOstreamOption::formatNames
// * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * * //
Foam::IOstreamOption::streamFormat
Foam::IOstreamOption::formatEnum
Foam::IOstreamOption::floatFormat
Foam::IOstreamOption::floatFormatEnum
(
const word& formatName,
const streamFormat deflt
const word& fmtName,
const floatFormat deflt
)
{
// Handle bad input graciously. A no-op for an empty string
if (!formatName.empty())
if (!fmtName.empty())
{
if (formatNames.contains(formatName))
const auto iter = floatFormatNames.cfind(fmtName);
if (iter.good())
{
return formatNames[formatName];
return iter.val();
}
// Fall-through to warning
WarningInFunction
<< "Unknown format specifier '" << formatName
<< "', using '" << formatNames[deflt] << "'\n";
<< "Unknown float format specifier '" << fmtName
<< "' using '" << floatFormatNames[deflt]
<< "' from " << floatFormatNames << nl;
}
return deflt;
}
Foam::IOstreamOption::floatFormat
Foam::IOstreamOption::floatFormatEnum
(
const word& key,
const dictionary& dict,
const floatFormat deflt
)
{
return floatFormatNames.getOrDefault(key, dict, deflt, true); // warnOnly
}
Foam::IOstreamOption::streamFormat
Foam::IOstreamOption::formatEnum
(
const word& fmtName,
const streamFormat deflt
)
{
// Handle bad input graciously. A no-op for an empty string
if (!fmtName.empty())
{
const auto iter = formatNames.cfind(fmtName);
if (iter.good())
{
return iter.val();
}
// Fall-through to warning
WarningInFunction
<< "Unknown stream format specifier '" << fmtName
<< "' using '" << formatNames[deflt]
<< "' from " << formatNames << nl;
}
return deflt;
@ -83,7 +139,7 @@ Foam::IOstreamOption::formatEnum
const streamFormat deflt
)
{
return formatNames.getOrDefault(key, dict, deflt, true); // failsafe=true
return formatNames.getOrDefault(key, dict, deflt, true); // warnOnly
}
@ -114,8 +170,7 @@ Foam::IOstreamOption::compressionEnum
WarningInFunction
<< "Unknown compression specifier '" << compName
<< "', using compression "
<< (deflt ? "on" : "off" ) << nl;
<< "' using compression " << (deflt ? "on" : "off") << nl;
}
return deflt;
@ -132,7 +187,7 @@ Foam::IOstreamOption::compressionEnum
{
return
(
Switch(key, dict, Switch(bool(deflt)), true) // failsafe=true
Switch(key, dict, Switch(bool(deflt)), true) // warnOnly
? compressionType::COMPRESSED
: compressionType::UNCOMPRESSED
);

135
src/OpenFOAM/db/options/IOstreamOption.H
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2015 OpenFOAM Foundation
Copyright (C) 2018-2022 OpenCFD Ltd.
Copyright (C) 2018-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -48,6 +48,7 @@ SourceFiles
#define Foam_IOstreamOption_H
#include "word.H"
#include <ios>
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -97,6 +98,17 @@ public:
ATOMIC //!< atomic = true
};
//- Float formats (eg, time directory name formats)
enum class floatFormat : unsigned
{
//! default float notation
general = unsigned(0),
//! fixed-point notation
fixed = unsigned(std::ios_base::fixed),
//! scientific notation
scientific = unsigned(std::ios_base::scientific)
};
//- Representation of a major/minor version number
class versionNumber
@ -134,7 +146,7 @@ public:
explicit versionNumber(const token& tok);
//- Failsafe construct from dictionary lookup.
versionNumber(const word& keyword, const dictionary& dict);
versionNumber(const word& key, const dictionary& dict);
// Member Functions
@ -168,13 +180,16 @@ public:
// Positive when 'this' is greater than other.
int compare(const versionNumber& other) const noexcept
{
return number_ - other.number_;
return (number_ - other.number_);
}
};
// Public Static Data
//- Names for float formats (general, fixed, scientific)
static const Enum<floatFormat> floatFormatNames;
//- Stream format names (ascii, binary)
static const Enum<streamFormat> formatNames;
@ -182,6 +197,75 @@ public:
static const versionNumber currentVersion;
// Static Helpers
//- Lookup floatFormat enum corresponding to the string
//- (general | fixed | scientific).
//
// If the string is not recognized, emit warning and return default.
// Silent if the string itself is empty.
//
// \note Can be used as constructor substitute for the enumeration
static floatFormat floatFormatEnum
(
const word& fmtName,
const floatFormat deflt = floatFormat::general
);
//- getOrDefault floatFormat from dictionary,
//- warn only on bad enumeration.
static floatFormat floatFormatEnum
(
const word& key, //!< Lookup key. Uses LITERAL (not REGEX)
const dictionary& dict, //!< dictionary
const floatFormat deflt = floatFormat::general
);
//- Lookup streamFormat enum corresponding to the string
//- (ascii | binary).
//
// If the string is not recognized, emit warning and return default.
// Silent if the string itself is empty.
//
// \note Can be used as constructor substitute for the enumeration
static streamFormat formatEnum
(
const word& fmtName,
const streamFormat deflt = streamFormat::ASCII
);
//- getOrDefault streamFormat from dictionary,
//- warn only on bad enumeration.
static streamFormat formatEnum
(
const word& key, //!< Lookup key. Uses LITERAL (not REGEX)
const dictionary& dict, //!< dictionary
const streamFormat deflt = streamFormat::ASCII
);
//- The compression enum corresponding to the string.
// Expects switch values (true/false, on/off, ...)
//
// If the string is not recognized, emit warning and return default.
// Silent if the string itself is empty.
//
// \note Can be used as constructor substitute for the enumeration
static compressionType compressionEnum
(
const word& compName,
const compressionType deflt = compressionType::UNCOMPRESSED
);
//- getOrDefault compressionType from dictionary,
//- warn only on bad enumeration.
static compressionType compressionEnum
(
const word& key, //!< Lookup key. Uses LITERAL (not REGEX)
const dictionary& dict, //!< dictionary
const compressionType deflt = compressionType::UNCOMPRESSED
);
private:
// Private Data
@ -252,51 +336,6 @@ public:
{}
// Static Member Functions
//- The stream format enum corresponding to the string
//- (ascii | binary).
//
// If the string is not recognized, emit warning and return default.
// Silent if the string itself is empty.
//
// \note Can be used as constructor substitute for the enumeration
static streamFormat formatEnum
(
const word& formatName,
const streamFormat deflt = streamFormat::ASCII
);
//- Failsafe construct streamFormat from optional dictionary lookup
static streamFormat formatEnum
(
const word& key, //!< Lookup key. Uses LITERAL (not REGEX)
const dictionary& dict, //!< dictionary
const streamFormat deflt = streamFormat::ASCII
);
//- The compression enum corresponding to the string.
// Expects switch values (true/false, on/off, ...)
//
// If the string is not recognized, emit warning and return default.
// Silent if the string itself is empty.
//
// \note Can be used as constructor substitute for the enumeration
static compressionType compressionEnum
(
const word& compName,
const compressionType deflt = compressionType::UNCOMPRESSED
);
//- Failsafe construct compressionType from optional dictionary lookup
static compressionType compressionEnum
(
const word& key, //!< Lookup key. Uses LITERAL (not REGEX)
const dictionary& dict, //!< dictionary
const compressionType deflt = compressionType::UNCOMPRESSED
);
// Member Functions
//- Get the current stream format

58
src/fileFormats/ensight/file/ensightCase.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -56,6 +56,44 @@ Foam::word Foam::ensightCase::padded(const int nwidth, const label value)
}
void Foam::ensightCase::setTimeFormat
(
OSstream& os,
IOstreamOption::floatFormat timeFmt,
const int timePrec
)
{
os.setf(std::ios_base::left);
os.setf
(
std::ios_base::fmtflags(timeFmt),
std::ios_base::floatfield
);
if (timePrec > 0)
{
os.precision(timePrec);
}
}
void Foam::ensightCase::setTimeFormat
(
OSstream& os,
const ensightCase::options& opts
)
{
os.setf(std::ios_base::left);
os.setf
(
std::ios_base::fmtflags(opts.timeFormat()),
std::ios_base::floatfield
);
os.precision(opts.timePrecision());
}
void Foam::ensightCase::printTimeset
(
OSstream& os,
@ -184,7 +222,7 @@ Foam::fileName Foam::ensightCase::dataDir() const
void Foam::ensightCase::initialize()
{
if (Pstream::master())
if (UPstream::master())
{
// EnSight and EnSight/data directories must exist
@ -211,11 +249,7 @@ void Foam::ensightCase::initialize()
// The case file is always ASCII
os_.reset(new OFstream(ensightDir_/caseName_, IOstreamOption::ASCII));
// Format options
os_->setf(ios_base::left);
os_->setf(ios_base::scientific, ios_base::floatfield);
os_->precision(5);
ensightCase::setTimeFormat(*os_, *options_); // Format options
writeHeader();
}
@ -468,7 +502,7 @@ Foam::ensightCase::createDataFile
const word& name
) const
{
if (Pstream::master())
if (UPstream::master())
{
// The data/ITER subdirectory must exist
// Note that data/ITER is indeed a valid ensight::FileName
@ -490,7 +524,7 @@ Foam::ensightCase::createCloudFile
const word& name
) const
{
if (Pstream::master())
if (UPstream::master())
{
// Write
// eg -> "data/********/lagrangian/<cloudName>/positions"
@ -584,7 +618,7 @@ void Foam::ensightCase::setTime(const scalar value, const label index)
timeIndex_ = index;
timeValue_ = value;
if (Pstream::master())
if (UPstream::master())
{
// The data/ITER subdirectory must exist
// Note that data/ITER is indeed a valid ensight::FileName
@ -810,7 +844,7 @@ Foam::ensightCase::newGeometry
{
autoPtr<Foam::ensightGeoFile> output;
if (Pstream::master())
if (UPstream::master())
{
// Set the path of the ensight file
fileName path;
@ -844,7 +878,7 @@ Foam::ensightCase::newCloud
{
autoPtr<Foam::ensightFile> output;
if (Pstream::master())
if (UPstream::master())
{
output = createCloudFile(cloudName, "positions");

74
src/fileFormats/ensight/file/ensightCase.H
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -57,6 +57,7 @@ namespace Foam
// Forward Declarations
class bitSet;
class dictionary;
class ensightCase;
class instant;
class OSstream;
@ -333,6 +334,21 @@ public:
// Output Helpers
//- Set output time format for ensight case file
static void setTimeFormat
(
OSstream& os,
IOstreamOption::floatFormat timeFmt,
const int timePrec
);
//- Set output time format for ensight case file
static void setTimeFormat
(
OSstream& os,
const ensightCase::options& opts
);
//- Print time-set for ensight case file with a single time
static void printTimeset
(
@ -388,7 +404,7 @@ class ensightCase::options
{
// Private Data
//- Ascii/Binary file output
//- The output file format (ascii/binary)
IOstreamOption::streamFormat format_;
//- Remove existing directory and sub-directories on creation
@ -400,8 +416,14 @@ class ensightCase::options
//- Write clouds into their own directory
bool separateCloud_;
//- Time format for case file (default: scientific)
IOstreamOption::floatFormat timeFormat_;
//- Time precision for case file (default: 5)
int timePrecision_;
//- Width of mask for subdirectories
label width_;
int width_;
//- The '*' mask appropriate for subdirectories
word mask_;
@ -415,16 +437,39 @@ public:
// Constructors
//- Construct with the specified format (default is binary)
options(IOstreamOption::streamFormat fmt = IOstreamOption::BINARY);
explicit options
(
IOstreamOption::streamFormat fmt = IOstreamOption::BINARY
);
//- If present, construct with the format specified in the dictionary
//- or use default (binary)
options
(
//! The lookup name for the format, typically 'format'
//! or 'writeFormat' etc.
const word& formatKeyword,
const dictionary& dict,
IOstreamOption::streamFormat fmt = IOstreamOption::BINARY
);
// Member Functions
// Access
//- Ascii/Binary file output
//- The output file format (ascii/binary)
IOstreamOption::streamFormat format() const noexcept { return format_; }
//- Time format for case file (general/fixed/scientific)
IOstreamOption::floatFormat timeFormat() const noexcept
{
return timeFormat_;
}
//- Time precision for case file
int timePrecision() const noexcept { return timePrecision_; }
//- The '*' mask appropriate for sub-directories
const word& mask() const noexcept { return mask_; }
@ -432,7 +477,7 @@ public:
word padded(const label i) const;
//- Return current width of mask and padded.
label width() const noexcept { return width_; }
int width() const noexcept { return width_; }
//- Remove existing directory and sub-directories on creation
bool overwrite() const noexcept { return overwrite_; }
@ -445,7 +490,22 @@ public:
//- Set width of mask and padded.
// Default width is 8 digits, max width is 31 digits.
void width(const label i);
void width(const int i);
//- Set the time format for case file
void timeFormat(IOstreamOption::floatFormat fmt) noexcept
{
timeFormat_ = fmt;
}
//- Set the time precision for case file
void timePrecision(int prec) noexcept { timePrecision_ = prec; }
//- Set the time format for case file
void timeFormat(const word& key, const dictionary& dict);
//- Set the time precision for case file
void timePrecision(const word& key, const dictionary& dict);
//- Remove existing directory and sub-directories on creation
void overwrite(bool on) noexcept { overwrite_ = on; }

49
src/fileFormats/ensight/file/ensightCaseOptions.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -26,15 +26,24 @@ License
\*---------------------------------------------------------------------------*/
#include "ensightCase.H"
#include "dictionary.H"
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::ensightCase::options::options(IOstreamOption::streamFormat fmt)
:
format_(fmt),
format_
(
// Can only be ASCII or BINARY
(fmt == IOstreamOption::streamFormat::ASCII)
? IOstreamOption::streamFormat::ASCII
: IOstreamOption::streamFormat::BINARY
),
overwrite_(false),
nodeValues_(false),
separateCloud_(false),
timeFormat_(IOstreamOption::floatFormat::scientific),
timePrecision_(5),
width_(0),
mask_(),
printf_()
@ -43,6 +52,17 @@ Foam::ensightCase::options::options(IOstreamOption::streamFormat fmt)
}
Foam::ensightCase::options::options
(
const word& formatKeyword,
const dictionary& dict,
IOstreamOption::streamFormat fmt
)
:
options(IOstreamOption::formatEnum(formatKeyword, dict, fmt))
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
Foam::word Foam::ensightCase::options::padded(const label i) const
@ -58,7 +78,7 @@ Foam::word Foam::ensightCase::options::padded(const label i) const
}
void Foam::ensightCase::options::width(const label n)
void Foam::ensightCase::options::width(const int n)
{
// Enforce min/max sanity limits
if (n < 1 || n > 31)
@ -74,4 +94,27 @@ void Foam::ensightCase::options::width(const label n)
}
void Foam::ensightCase::options::timeFormat
(
const word& key,
const dictionary& dict
)
{
timeFormat_ = IOstreamOption::floatFormatEnum(key, dict, timeFormat_);
}
void Foam::ensightCase::options::timePrecision
(
const word& key,
const dictionary& dict
)
{
if (!key.empty())
{
dict.readIfPresent(key, timePrecision_, keyType::LITERAL);
}
}
// ************************************************************************* //

24
src/fileFormats/ensight/file/ensightCaseTemplates.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2020 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -40,19 +40,16 @@ Foam::ensightCase::newData
{
autoPtr<ensightFile> output;
if (Pstream::master())
if (UPstream::master())
{
const ensight::VarName varName(name);
output = createDataFile(varName);
// Description
output().write
output().writeString
(
string
(
padded(timeIndex_) / varName
+ " <" + pTraits<Type>::typeName + ">"
)
padded(timeIndex_) / varName
+ " <" + pTraits<Type>::typeName + ">"
);
output().newline();
@ -91,19 +88,16 @@ Foam::ensightCase::newCloudData
{
autoPtr<ensightFile> output;
if (Pstream::master())
if (UPstream::master())
{
const ensight::VarName varName(name);
output = createCloudFile(cloudName, varName);
// Description
output().write
output().writeString
(
string
(
padded(timeIndex_) / cloudName / varName
+ " <" + pTraits<Type>::typeName + ">"
)
padded(timeIndex_) / cloudName / varName
+ " <" + pTraits<Type>::typeName + ">"
);
output().newline();

46
src/fileFormats/ensight/file/ensightFile.C
View File

@ -147,7 +147,7 @@ float Foam::ensightFile::undefValue(float value) noexcept
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
Foam::Ostream& Foam::ensightFile::writeString(const char* str, size_t len)
void Foam::ensightFile::writeString(const char* str, size_t len)
{
// Output 79 chars (ASCII) or 80 chars (BINARY)
char buf[80];
@ -165,41 +165,47 @@ Foam::Ostream& Foam::ensightFile::writeString(const char* str, size_t len)
else
{
buf[79] = 0; // Max 79 in ASCII
// TBD: Extra safety - trap newline in ASCII?
// char* p = ::strchr(buf, '\n');
// if (p) *p = 0;
stdStream() << buf;
syncState();
}
return *this;
}
Foam::Ostream& Foam::ensightFile::writeString(const char* str)
void Foam::ensightFile::writeString(const char* str)
{
return writeString(str, strlen(str));
writeString(str, strlen(str));
}
Foam::Ostream& Foam::ensightFile::writeString(const std::string& str)
void Foam::ensightFile::writeString(const std::string& str)
{
return writeString(str.data(), str.size());
writeString(str.data(), str.size());
}
Foam::Ostream& Foam::ensightFile::write(const char* str)
{
return writeString(str, strlen(str));
writeString(str, strlen(str));
return *this;
}
Foam::Ostream& Foam::ensightFile::write(const word& str)
{
return writeString(str.data(), str.size());
writeString(str.data(), str.size());
return *this;
}
Foam::Ostream& Foam::ensightFile::write(const std::string& str)
{
return writeString(str.data(), str.size());
writeString(str.data(), str.size());
return *this;
}
@ -304,10 +310,9 @@ void Foam::ensightFile::newline()
}
Foam::Ostream& Foam::ensightFile::writeUndef()
void Foam::ensightFile::writeUndef()
{
write(undefValue_);
return *this;
}
@ -330,14 +335,27 @@ Foam::Ostream& Foam::ensightFile::writeKeyword(const keyType& key)
}
Foam::Ostream& Foam::ensightFile::writeBinaryHeader()
void Foam::ensightFile::writeBinaryHeader()
{
if (format() == IOstreamOption::BINARY)
{
writeString("C Binary");
// Is binary: newline() is a no-op
}
}
return *this;
void Foam::ensightFile::beginTimeStep()
{
writeString("BEGIN TIME STEP");
newline();
}
void Foam::ensightFile::endTimeStep()
{
writeString("END TIME STEP");
newline();
}

16
src/fileFormats/ensight/file/ensightFile.H
View File

@ -146,19 +146,25 @@ public:
// Output
//- Write "C Binary" string for binary files (eg, geometry/measured)
Ostream& writeBinaryHeader();
void writeBinaryHeader();
//- Write "BEGIN TIME STEP" string and newline
void beginTimeStep();
//- Write "END TIME STEP" string and newline
void endTimeStep();
//- Write character/string content as "%79s" or as binary (max 80 chars)
Ostream& writeString(const char* str, size_t len);
void writeString(const char* str, size_t len);
//- Write C-string as "%79s" or as binary (max 80 chars)
Ostream& writeString(const char* str);
void writeString(const char* str);
//- Write string as "%79s" or as binary (max 80 chars)
Ostream& writeString(const std::string& str);
void writeString(const std::string& str);
//- Write undef value
Ostream& writeUndef();
void writeUndef();
//- Write element keyword with trailing newline,

19
src/fileFormats/ensight/file/ensightGeoFile.C
View File

@ -34,23 +34,32 @@ License
void Foam::ensightGeoFile::init()
{
writeBinaryHeader();
beginGeometry();
}
void Foam::ensightGeoFile::beginGeometry()
{
// Description line 1
write("Ensight Geometry File");
writeString("Ensight Geometry File");
newline();
// Description line 2
write(string("Written by OpenFOAM " + std::to_string(foamVersion::api)));
writeString("Written by OpenFOAM " + std::to_string(foamVersion::api));
newline();
write("node id assign");
writeString("node id assign");
newline();
write("element id assign");
writeString("element id assign");
newline();
}
void Foam::ensightGeoFile::endGeometry()
{}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::ensightGeoFile::ensightGeoFile
@ -96,7 +105,7 @@ Foam::Ostream& Foam::ensightGeoFile::writeKeyword(const keyType& key)
void Foam::ensightGeoFile::beginPart
(
const label index,
const string& description
const std::string& description
)
{
beginPart(index);

11
src/fileFormats/ensight/file/ensightGeoFile.H
View File

@ -52,16 +52,21 @@ class ensightGeoFile
{
// Private Member Functions
//- Initialize outputs the header information
//- Initialize outputs the header information and beginGeometry
void init();
//- Start of geometry information
void beginGeometry();
//- End of geometry information
void endGeometry();
//- No copy construct
ensightGeoFile(const ensightGeoFile&) = delete;
//- No copy assignment
void operator=(const ensightGeoFile&) = delete;
public:
// Static Functions
@ -111,7 +116,7 @@ public:
using ensightFile::beginPart;
//- Begin a "part" (0-based index), with a description.
void beginPart(const label index, const string& description);
void beginPart(const label index, const std::string& description);
//- Begin a "coordinates" block
void beginCoordinates(const label npoints);

49
src/fileFormats/ensight/file/ensightWriterCaching.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -27,6 +27,7 @@ License
#include "ensightWriterCaching.H"
#include "ListOps.H"
#include "OTstream.H"
#include "Fstream.H"
// * * * * * * * * * * * * * * * Local Functions * * * * * * * * * * * * * * //
@ -129,15 +130,9 @@ Foam::label Foam::ensightOutput::writerCaching::readPreviousTimes
const scalar timeValue
)
{
// In 1906 and earlier, the fieldsDict contained "meshes" and "times"
// entries, each with their own time values.
// This makes it more difficult to define the exact correspondence
// between geometry intervals and times.
//
// Now track the used geometry intervals as a bitSet.
// Track the used geometry intervals as a bitSet
// Only called from master
// Note: only called from master
label timeIndex = 0;
cache_.clear();
@ -158,21 +153,6 @@ Foam::label Foam::ensightOutput::writerCaching::readPreviousTimes
// Convert indices to bitSet entries
geoms_.set(geomIndices);
}
else if (cache_.readIfPresent("meshes", meshTimes))
{
WarningInFunction
<< nl
<< "Setting geometry timeset information from time values"
<< " (cache from an older OpenFOAM version)." << nl
<< "This may not be fully reliable." << nl
<< nl;
for (const scalar meshTime : meshTimes)
{
const label geomIndex = findTimeIndex(times_, meshTime);
geoms_.set(geomIndex);
}
}
// Make length consistent with time information.
// We read/write the indices instead of simply dumping the bitSet.
@ -263,17 +243,19 @@ bool Foam::ensightOutput::writerCaching::update
geoms_.set(timeIndex);
}
// Update time/geometry information in dictionary
cache_.set("times", times_);
cache_.set("geometry", geoms_.sortedToc());
// Debugging, or if needed for older versions:
//// cache_.set
//// (
//// "meshes",
//// IndirectList<scalar>(times_, geoms_.sortedToc())
//// );
// Note: to avoid inadvertent loss of precision,
// generate output tokens for the list of times directly
{
// Same as: cache_.set("times", times_);
OTstream os;
os << times_;
tokenList toks(std::move(os.tokens()));
cache_.set(new primitiveEntry("times", std::move(toks)));
}
cache_.set("geometry", geoms_.sortedToc());
// Add field information to dictionary
dictionary& dict = fieldDict(fieldName);
@ -293,6 +275,7 @@ bool Foam::ensightOutput::writerCaching::update
if (stateChanged)
{
OFstream os(dictFile);
os.precision(16); // increased precision to avoid rounding
os << "// State file for writer output" << nl << nl;
cache_.write(os, false);

10
src/functionObjects/utilities/ensightWrite/ensightWrite.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -97,10 +97,7 @@ Foam::functionObjects::ensightWrite::ensightWrite
:
fvMeshFunctionObject(name, runTime, dict),
writeOpts_(),
caseOpts_
(
IOstreamOption::formatEnum("format", dict, runTime.writeFormat())
),
caseOpts_("format", dict, IOstreamOption::BINARY),
outputDir_(),
consecutive_(false),
meshState_(polyMesh::TOPO_CHANGE),
@ -180,6 +177,9 @@ bool Foam::functionObjects::ensightWrite::read(const dictionary& dict)
caseOpts_.width(dict.getOrDefault<label>("width", 8));
caseOpts_.overwrite(dict.getOrDefault("overwrite", false));
caseOpts_.timeFormat("timeFormat", dict);
caseOpts_.timePrecision("timePrecision", dict);
// Output directory

9
src/functionObjects/utilities/ensightWrite/ensightWrite.H
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2016-2022 OpenCFD Ltd.
Copyright (C) 2016-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -42,6 +42,9 @@ Description
writeInterval 1;
format binary;
timeFormat scientific;
timePrecision 5;
overwrite true;
width 12;
@ -95,11 +98,13 @@ Description
\heading Ensight Output Options
\table
Property | Description | Required | Default
format | ascii or binary format | no | same as simulation
format | ascii or binary format | no | binary
width | Mask width for \c data/XXXX | no | 8
directory | The output directory name | no | postProcessing/NAME
overwrite | Remove existing directory | no | false
consecutive | Consecutive output numbering | no | false
timeFormat | Time format (ensight case) | no | scientific
timePrecision | Time precision (ensight case) | no | 5
\endtable
\heading Output Selection

14
src/meshTools/coordSet/writers/ensight/ensightCoordSetWriter.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation
Copyright (C) 2021-2022 OpenCFD Ltd.
Copyright (C) 2021-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -176,7 +176,7 @@ void Foam::coordSetWriters::ensightWriter::writeGeometry
Foam::coordSetWriters::ensightWriter::ensightWriter()
:
coordSetWriter(),
writeFormat_(IOstreamOption::ASCII),
caseOpts_(IOstreamOption::BINARY),
collateTimes_(true),
caching_("fieldsDict") // Historic name
{}
@ -185,13 +185,13 @@ Foam::coordSetWriters::ensightWriter::ensightWriter()
Foam::coordSetWriters::ensightWriter::ensightWriter(const dictionary& options)
:
coordSetWriter(options),
writeFormat_
(
IOstreamOption::formatEnum("format", options, IOstreamOption::ASCII)
),
caseOpts_("format", options, IOstreamOption::BINARY),
collateTimes_(options.getOrDefault("collateTimes", true)),
caching_("fieldsDict") // Historic name
{}
{
caseOpts_.timeFormat("timeFormat", options);
caseOpts_.timePrecision("timePrecision", options);
}
Foam::coordSetWriters::ensightWriter::ensightWriter

12
src/meshTools/coordSet/writers/ensight/ensightCoordSetWriter.H
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation
Copyright (C) 2021-2022 OpenCFD Ltd.
Copyright (C) 2021-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -44,8 +44,11 @@ Description
Format options:
\table
Property | Description | Required | Default
format | ascii/binary | no | ascii
format | ascii/binary | no | binary
collateTimes | use common geometry for times | no | true
timeFormat | Time format (ensight case) | no | scientific
timePrecision | Time precision (ensight case) | no | 5
\endtable
\endtable
SourceFiles
@ -57,6 +60,7 @@ SourceFiles
#define Foam_coordSetWriters_ensightWriter_H
#include "coordSetWriter.H"
#include "ensightCase.H"
#include "ensightWriterCaching.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -91,8 +95,8 @@ class ensightWriter
// Private Data
//- Output format option (default: ASCII)
IOstreamOption::streamFormat writeFormat_;
//- Ensight case options
ensightCase::options caseOpts_;
//- Collate times (default: true)
bool collateTimes_;

12
src/meshTools/coordSet/writers/ensight/ensightCoordSetWriterCollated.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2022 OpenCFD Ltd.
Copyright (C) 2022-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -148,7 +148,7 @@ Foam::fileName Foam::coordSetWriters::ensightWriter::writeCollated
(
geomFile.path(),
geomFile.name(),
writeFormat_
caseOpts_.format()
);
writeGeometry(osGeom, elemOutput);
@ -160,7 +160,7 @@ Foam::fileName Foam::coordSetWriters::ensightWriter::writeCollated
(
dataDir,
varName,
writeFormat_
caseOpts_.format()
);
if (verbose_)
@ -177,11 +177,7 @@ Foam::fileName Foam::coordSetWriters::ensightWriter::writeCollated
if (stateChanged)
{
OFstream osCase(outputFile, IOstreamOption::ASCII);
// Format options
osCase.setf(ios_base::left);
osCase.setf(ios_base::scientific, ios_base::floatfield);
osCase.precision(5);
ensightCase::setTimeFormat(osCase, caseOpts_); // time-format
if (verbose_)
{

12
src/meshTools/coordSet/writers/ensight/ensightCoordSetWriterUncollated.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2022 OpenCFD Ltd.
Copyright (C) 2022-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -101,13 +101,13 @@ Foam::fileName Foam::coordSetWriters::ensightWriter::writeUncollated
(
baseDir,
baseName + ".00000000.mesh",
writeFormat_
caseOpts_.format()
);
ensightFile osField
(
baseDir,
baseName + ".00000000." + varName,
writeFormat_
caseOpts_.format()
);
writeGeometry(osGeom, elemOutput);
@ -119,11 +119,7 @@ Foam::fileName Foam::coordSetWriters::ensightWriter::writeUncollated
// Update case file
{
OFstream osCase(outputFile, IOstreamOption::ASCII);
// Format options
osCase.setf(ios_base::left);
osCase.setf(ios_base::scientific, ios_base::floatfield);
osCase.precision(5);
ensightCase::setTimeFormat(osCase, caseOpts_); // time-format
osCase
<< "FORMAT" << nl

62
src/randomProcesses/noise/noiseModels/noiseModel/noiseModel.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -185,29 +185,34 @@ Foam::scalar Foam::noiseModel::checkUniformTimeStep
const scalarList& times
) const
{
scalar deltaT = -1.0;
scalar deltaT = 1;
if (times.size() > 1)
{
// Uniform time step
deltaT = (times.last() - times.first())/scalar(times.size() - 1);
deltaT = (times.back() - times.front())/scalar(times.size() - 1);
for (label i = 1; i < times.size(); ++i)
bool nonUniform = false;
for (label i = 1; i < times.size() && !nonUniform; ++i)
{
scalar dT = times[i] - times[i-1];
const scalar dT = times[i] - times[i-1];
if (mag(dT/deltaT - 1) > 1e-8)
{
FatalErrorInFunction
<< "Unable to process data with a variable time step"
<< exit(FatalError);
}
nonUniform = (mag(dT/deltaT - 1) > 1e-8);
}
if (nonUniform)
{
WarningInFunction
<< "Detected non-uniform time step:"
<< " resulting FFT may be incorrect"
<< " (or the deltaT is just very small): " << deltaT
<< endl;
}
}
else
{
FatalErrorInFunction
<< "Unable to create FFT with a single value"
<< "Unable to create FFT with 0 or 1 values"
<< exit(FatalError);
}
@ -617,6 +622,7 @@ Foam::noiseModel::noiseModel
nSamples_(65536),
fLower_(25),
fUpper_(10000),
sampleFreq_(0),
startTime_(0),
windowModelPtr_(),
SPLweighting_(weightingType::none),
@ -654,10 +660,16 @@ bool Foam::noiseModel::read(const dictionary& dict)
return false;
}
// Re-assign defaults (to avoid stickiness)
fLower_ = 25;
fUpper_ = 10000;
sampleFreq_ = 0;
dict.readIfPresent("rhoRef", rhoRef_);
dict.readIfPresent("N", nSamples_);
dict.readIfPresent("fl", fLower_);
dict.readIfPresent("fu", fUpper_);
dict.readIfPresentCompat("minFreq", {{"fl", 2312}}, fLower_);
dict.readIfPresentCompat("maxFreq", {{"fu", 2312}}, fUpper_);
dict.readIfPresent("sampleFreq", sampleFreq_);
dict.readIfPresent("startTime", startTime_);
dict.readIfPresent("minPressure", minPressure_);
dict.readIfPresent("maxPressure", maxPressure_);
@ -666,29 +678,39 @@ bool Foam::noiseModel::read(const dictionary& dict)
if (fLower_ < 0)
{
FatalIOErrorInFunction(dict)
<< "fl: lower frequency bound must be greater than zero"
<< "Lower frequency bound must be greater than zero"
<< exit(FatalIOError);
}
if (fUpper_ < 0)
{
FatalIOErrorInFunction(dict)
<< "fu: upper frequency bound must be greater than zero"
<< "Upper frequency bound must be greater than zero"
<< exit(FatalIOError);
}
if (fUpper_ < fLower_)
{
FatalIOErrorInFunction(dict)
<< "fu: upper frequency bound must be greater than lower "
<< "frequency bound (fl)"
<< "Upper frequency bound (" << fUpper_
<< ") must be greater than lower frequency bound ("
<< fLower_ << ")"
<< exit(FatalIOError);
}
Info<< " Frequency bounds:" << nl
<< " lower: " << fLower_ << nl
<< " upper: " << fUpper_ << endl;
<< " upper: " << fUpper_ << nl
<< " sample: ";
if (sampleFreq_ > 0)
{
Info<< sampleFreq_ << nl;
}
else
{
Info<< "auto" << nl;
}
weightingTypeNames_.readIfPresent("SPLweighting", dict, SPLweighting_);

17
src/randomProcesses/noise/noiseModels/noiseModel/noiseModel.H
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -34,8 +34,8 @@ Description
\verbatim
rhoRef 1;
N 4096;
fl 25;
fu 10000;
minFreq 25;
maxFreq 10000;
startTime 0;
outputPrefix "test1";
@ -58,8 +58,9 @@ Description
Property | Description | Required | Default value
rhoRef | Reference density | no | 1
N | Number of samples in sampling window | no | 65536 (2^16)
fl | Lower frequency bounds | no | 25
fu | Upper frequency bounds | no | 10000
minFreq | Lower frequency bounds (fl) | no | 25
maxFreq | Upper frequency bounds (fu) | no | 10000
sampleFreq | Sample frequency | no | 0
startTime | Start time | no | 0
outputPrefix | Prefix applied to output files| no | ''
SPLweighting | Weighting: dBA, dBB, dBC, DBD | no | none
@ -71,6 +72,9 @@ Description
writeOctaves | Write octaves data | no | yes
\endtable
If provided, the sampleFreq is used to define the deltaT between samples.
Otherwise the deltaT is determined from the time samples themselves.
SourceFiles
noiseModel.C
@ -158,6 +162,9 @@ protected:
//- Upper frequency limit, default = 10kHz
scalar fUpper_;
//- Prescribed sample frequency
scalar sampleFreq_;
//- Start time, default = 0s
scalar startTime_;

7
src/randomProcesses/noise/noiseModels/pointNoise/pointNoise.C
View File

@ -101,7 +101,12 @@ void pointNoise::processData
Info<< "Creating noise FFT" << endl;
const scalar deltaT = checkUniformTimeStep(t);
const scalar deltaT =
(
sampleFreq_ > 0
? (1.0/sampleFreq_)
: checkUniformTimeStep(t)
);
// Apply conversions
p *= rhoRef_;

18
src/randomProcesses/noise/noiseModels/surfaceNoise/surfaceNoise.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -96,15 +96,21 @@ void surfaceNoise::initialise(const fileName& fName)
// Note: all processors should call the windowing validate function
label nRequiredTimes = windowModelPtr_->validate(nAvailableTimes);
if (Pstream::master())
if (UPstream::master())
{
// Restrict times
times_.setSize(nRequiredTimes);
forAll(times_, timeI)
times_.resize_nocopy(nRequiredTimes);
forAll(times_, timei)
{
times_[timeI] = allTimes[timeI + startTimeIndex_].value();
times_[timei] = allTimes[timei + startTimeIndex_].value();
}
deltaT_ = checkUniformTimeStep(times_);
deltaT_ =
(
sampleFreq_ > 0
? (1.0/sampleFreq_)
: checkUniformTimeStep(times_)
);
// Read the surface geometry
// Note: hard-coded to read mesh from first time index

14
src/surfMesh/writers/ensight/ensightSurfaceWriter.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2014 OpenFOAM Foundation
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -55,7 +55,7 @@ namespace surfaceWriters
Foam::surfaceWriters::ensightWriter::ensightWriter()
:
surfaceWriter(),
writeFormat_(IOstreamOption::ASCII),
caseOpts_(IOstreamOption::BINARY),
collateTimes_(true),
caching_("fieldsDict") // Historic name
{}
@ -67,13 +67,13 @@ Foam::surfaceWriters::ensightWriter::ensightWriter
)
:
surfaceWriter(options),
writeFormat_
(
IOstreamOption::formatEnum("format", options, IOstreamOption::ASCII)
),
caseOpts_("format", options, IOstreamOption::BINARY),
collateTimes_(options.getOrDefault("collateTimes", true)),
caching_("fieldsDict") // Historic name
{}
{
caseOpts_.timeFormat("timeFormat", options);
caseOpts_.timePrecision("timePrecision", options);
}
Foam::surfaceWriters::ensightWriter::ensightWriter

11
src/surfMesh/writers/ensight/ensightSurfaceWriter.H
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011 OpenFOAM Foundation
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -44,12 +44,14 @@ Description
Format options for ensight:
\table
Property | Description | Required | Default
format | ascii/binary | no | ascii
format | ascii/binary | no | binary
collateTimes | Use common geometry for times | no | true
scale | Output geometry scaling | no | 1
transform | Output coordinate transform | no |
fieldLevel | Subtract field level before scaling | no | empty dict
fieldScale | Output field scaling | no | empty dict
timeFormat | Time format (ensight case) | no | scientific
timePrecision | Time precision (ensight case) | no | 5
\endtable
The collated format maintains an internal list of the known times
@ -65,6 +67,7 @@ SourceFiles
#define Foam_surfaceWriters_ensightWriter_H
#include "surfaceWriter.H"
#include "ensightCase.H"
#include "ensightWriterCaching.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -84,8 +87,8 @@ class ensightWriter
{
// Private Data
//- Output format option (default: ASCII)
IOstreamOption::streamFormat writeFormat_;
//- Ensight case options
ensightCase::options caseOpts_;
//- Collate times (default: true)
bool collateTimes_;

12
src/surfMesh/writers/ensight/ensightSurfaceWriterCollated.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2014 OpenFOAM Foundation
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -163,7 +163,7 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeCollated
(
geomFile.path(),
geomFile.name(),
writeFormat_
caseOpts_.format()
);
part.write(osGeom); // serial
}
@ -173,7 +173,7 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeCollated
(
dataDir,
varName,
writeFormat_
caseOpts_.format()
);
if (verbose_)
@ -191,11 +191,7 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeCollated
if (stateChanged)
{
OFstream osCase(outputFile, IOstreamOption::ASCII);
// Format options
osCase.setf(ios_base::left);
osCase.setf(ios_base::scientific, ios_base::floatfield);
osCase.precision(5);
ensightCase::setTimeFormat(osCase, caseOpts_); // time-format
if (verbose_)
{

18
src/surfMesh/writers/ensight/ensightSurfaceWriterUncollated.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2014 OpenFOAM Foundation
Copyright (C) 2015-2022 OpenCFD Ltd.
Copyright (C) 2015-2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -73,7 +73,7 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeUncollated()
(
outputDir,
baseName + ".00000000.mesh",
writeFormat_
caseOpts_.format()
);
ensightOutputSurface part
@ -85,7 +85,9 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeUncollated()
part.write(osGeom); // serial
// Update case file
OFstream osCase(outputFile);
OFstream osCase(outputFile, IOstreamOption::ASCII);
ensightCase::setTimeFormat(osCase, caseOpts_); // time-format
osCase
<< "FORMAT" << nl
<< "type: ensight gold" << nl
@ -176,13 +178,13 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeUncollated
(
baseDir,
baseName + ".00000000.mesh",
writeFormat_
caseOpts_.format()
);
ensightFile osField
(
baseDir,
baseName + ".00000000." + varName,
writeFormat_
caseOpts_.format()
);
// Ensight Geometry
@ -203,11 +205,7 @@ Foam::fileName Foam::surfaceWriters::ensightWriter::writeUncollated
// Update case file
{
OFstream osCase(outputFile, IOstreamOption::ASCII);
// Format options
osCase.setf(ios_base::left);
osCase.setf(ios_base::scientific, ios_base::floatfield);
osCase.precision(5);
ensightCase::setTimeFormat(osCase, caseOpts_); // time-format
osCase
<< "FORMAT" << nl