ITK  5.3.0
Insight Toolkit
Examples/IO/XML/itkParticleSwarmOptimizerSAXWriter.cxx
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#include "itkParticleSwarmOptimizerSAXWriter.h"
#include "itksys/SystemTools.hxx"
namespace itk
{
int
ParticleSwarmOptimizerSAXWriter::CanWriteFile(const char * name)
{
std::ofstream ofs(name);
int yes = ofs.is_open();
if (yes)
{
ofs.close();
}
return yes;
}
int
ParticleSwarmOptimizerSAXWriter::WriteFile()
{
try
{
if (!this->CanWriteFile(this->m_Filename.c_str()))
{
ExceptionObject e(__FILE__, __LINE__);
std::string message = "Cannot write to ";
message += this->m_Filename;
message += "!\n";
e.SetDescription(message.c_str());
throw e;
}
if (this->m_InputObject == nullptr)
{
itkExceptionMacro("Object to be written is null!\n");
}
std::ofstream ofs(this->m_Filename.c_str());
// start the tag 'optimizer'
ofs << "<optimizer type=\"ParticleSwarmOptimizer\"";
ofs << " NumberOfParticles=\""
<< this->m_InputObject->GetNumberOfParticles() << "\"";
ofs << " MaximumNumberOfIterations=\""
<< this->m_InputObject->GetMaximalNumberOfIterations() << "\"";
ofs << " InertiaCoefficient=\""
<< this->m_InputObject->GetInertiaCoefficient() << "\"";
ofs << " GlobalCoefficient=\""
<< this->m_InputObject->GetGlobalCoefficient() << "\"";
ofs << " PersonalCoefficient=\""
<< this->m_InputObject->GetPersonalCoefficient() << "\"";
ofs << " FunctionConvergenceTolerance=\""
<< this->m_InputObject->GetFunctionConvergenceTolerance() << "\"";
ofs << " ConvergedPercentageToStop=\""
<< this->m_InputObject->GetPercentageParticlesConverged() << "\"";
ofs << ">";
ofs << "\n";
// write the lower bound
ofs << " <bound id=\"lower\"";
ofs << " value=\"";
for (const auto & bound : this->m_InputObject->GetParameterBounds())
{
ofs << " " << bound.first;
}
ofs << "\"";
ofs << "/>";
ofs << "\n";
// write the upper bound
ofs << " <bound id=\"upper\"";
ofs << " value=\"";
for (const auto & bound : this->m_InputObject->GetParameterBounds())
{
ofs << " " << bound.second;
}
ofs << "\"";
ofs << "/>";
ofs << "\n";
// write the ParametersConvergenceTolerance
ofs << " <ParametersConvergenceTolerance>";
Array<double> ptols =
this->m_InputObject->GetParametersConvergenceTolerance();
// Note: The data-cast to unsigned int is required
// because itk::Array only supports 'unsigned int' number of
// elements.
for (unsigned int i = 0; i < ptols.GetSize(); i++)
{
ofs << " " << ptols[i];
}
ofs << "</ParametersConvergenceTolerance>";
ofs << "\n";
// close the tag 'optimizer'
ofs << "</optimizer>";
ofs << "\n";
ofs.close();
return 1;
}
catch (...)
{
return 0;
}
}
} // namespace itk
itk
The "itk" namespace contains all Insight Segmentation and Registration Toolkit (ITK) classes....
Definition: itkAnnulusOperator.h:24
itk::Math::e
static constexpr double e
Definition: itkMath.h:54