mfileswriter.hh

Go to the documentation of this file.

/*****************************************************************************
*   This program is part of the msbGrid software                            *
*                                                                           *
*   msbGrid stands for multi-structured block Grid generator                *
*                                                                           *
*   msbGrid is a free software: you can redistribute it and/or modify       *
*   it under the terms of the GNU General Public License as published by    *
*   the Free Software Foundation, either version 2 of the License, or       *
*   (at your option) any later version.                                     *
*                                                                           *
*   msbGrid is distributed in the hope that it will be useful,              *
*   but WITHOUT ANY WARRANTY; without even the implied warranty of          *
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the           *
*   GNU General Public License for more details.                            *
*                                                                           *
*   See the file COPYING for full copying permissions.                      *
*****************************************************************************/
#ifndef _M_FILES_WRITER_HH
#define _M_FILES_WRITER_HH

#include <fstream>
#include <iomanip>
#include <cstdlib>

namespace msbGrid
{
class MFilesWriter;
}

class msbGrid :: MFilesWriter
{
public:
  typedef msbGrid :: Point Point;
  typedef msbGrid :: InputData InputData;
  typedef msbGrid :: Color Color;

private:
  const InputData *inputData_;
  std :: ofstream outputFile_;

public:

  MFilesWriter(const InputData * inputData)
  {
    inputData_ = inputData;
  }

  ~MFilesWriter() {}

  void openFile(const std :: string &name)
  {
    std :: string fileName = inputData_->outputPath() + "m/" + name + ".m";
#if VERBOSE_LEVEL >= 2
    std :: cout << "fileName = " << fileName << "\n";
#endif
    outputFile_.open(fileName.c_str(), std :: ios :: out);
    if ( !outputFile_ ) {
      std :: cerr << "Can't open output file" << fileName << std :: endl;
      std :: exit(1);
    }

    outputFile_ << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n";
    outputFile_ << "%%  This file was generated by msbGrid software %%\n";
    outputFile_ << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n\n";
    outputFile_ << "clear all; close all; clc; hold on; axis square;\n\n";
    writeCommanLine_();
  }

  void closeFile()
  { outputFile_.close(); }

  void writeCommanLine_()
  {
    outputFile_ << "%% command line options : ";
    outputFile_ << " --outputPath=\"" << inputData_->outputPath() << "\""
                << " --outputDirName=\"" << inputData_->outputDirName() << "\""
                << " --xMin=\"" << inputData_->xMin() << "\""
                << " --xMax=\"" << inputData_->xMax() << "\""
                << " --bnCx=\"" << inputData_->blockNumberController() [ 0 ] << "\""
                << " --bnCy=\"" << inputData_->blockNumberController() [ 1 ] << "\""
                << " --dr=\"" << inputData_->dr() << "\""
                << " --rafC=\"" << inputData_->rafCoeff() << "\""
                << " --bCId=\"" << inputData_->blockCenterRandSeqId() << "\""
                << " --bAId=\"" << inputData_->blockAngleRandSeqId() << "\""
                << " --bCPC=\"" << inputData_->blockCenterPertCoeff() << "\""
                << " --d2bC=\"" << inputData_->dist2BoundCoeff() << "\""
                << " --intbC=\"" << inputData_->interBlockDistCoeff() << "\""
                << "\n\n";
  }

  void write(const std :: string &texte)
  { outputFile_ << texte; }

  void write(const Point &pt0, const bool &writeText)
  {
#if !( ( DEV_MOD == 1 ) && ( CHECK_MARKED_POINTS == 1 ) )
    if ( pt0.mark() ) {
      return;
    }

#endif

    outputFile_ << std :: setw(1)
                << std :: fixed
                << "iPt = " << pt0.Id() << "; "
                << "x(iPt) = " << pt0.pos() [ 0 ] << "; "
                << "y(iPt) = " << pt0.pos() [ 1 ] << "; "
                << "plot(x(iPt),y(iPt)," << pt0.marker().symbol
                << ",'MarkerSize'," << pt0.marker().size
                << ",'MarkerEdgeColor'," << colorSymbol(pt0.marker().edgeColor)
                << ",'MarkerFaceColor'," << colorSymbol(pt0.marker().faceColor)
                << "); ";

    if ( writeText ) {
      outputFile_ << "text(x(iPt), y(iPt), num2str(iPt));";
    }

    outputFile_ << std :: endl;
  }

  void write(const msbGrid :: Edge &ed0, const bool &writeText)
  {
    if ( ed0.mark() ) {
      return;
    }

    outputFile_ << std :: setw(1)
                << std :: fixed
                << "iEd = " << ed0.Id() << "; "
                << "u = [x(" << ed0.vPt() [ 0 ]->Id() << ") x(" << ed0.vPt() [ 1 ]->Id() << ")]; "
                << "v = [y(" << ed0.vPt() [ 0 ]->Id() << ") y(" << ed0.vPt() [ 1 ]->Id() << ")]; "
                << "plot(u,v," << colorSymbol( ed0.color() ) << "); ";

    if ( writeText ) {
      outputFile_ << "xbc(iEd) = " << ed0.pos() [ 0 ] << "; "
                  << "ybc(iEd) = " << ed0.pos() [ 1 ] << "; "
                  << "text(xbc(iEd), ybc(iEd), num2str(iEd));";
    }

    outputFile_ << std :: endl;
  }

  void write(const msbGrid :: Triangle &tr0, const bool &writeText)
  {
    if ( tr0.mark() ) {
      return;
    }

    const int &Id0 = tr0.vPt() [ 0 ]->Id();
    const int &Id1 = tr0.vPt() [ 1 ]->Id();
    const int &Id2 = tr0.vPt() [ 2 ]->Id();

    outputFile_ << std :: setw(1)
                << std :: fixed
                << "iTr =" << tr0.Id() << "; "
                << "u = [x(" << Id0 << ") x(" << Id1 << ") x(" << Id2 << ") x(" << Id0 << ")]; "
                << "v = [y(" << Id0 << ") y(" << Id1 << ") y(" << Id2 << ") y(" << Id0 << ")]; "
                << "plot(u, v, 'k'); ";

    if ( tr0.color() != Color :: Black ) {
      outputFile_ << "fill(u,v," << colorSymbol( tr0.color() ) << "); ";
    }

    if ( writeText ) {
      outputFile_ << "xbc(iTr) = ( " << tr0.pos() [ 0 ] << "); "
                  << "ybc(iTr) = ( " << tr0.pos() [ 1 ] << "); "
                  << "text(xbc(iTr), ybc(iTr), num2str(iTr));";
    }

    outputFile_ << std :: endl;
  }

  void write(const msbGrid :: Quadrangle &qd0, const bool &writeText)
  {
    if ( qd0.mark() ) {
      return;
    }

    const int &Id0 = qd0.vPt() [ 0 ]->Id();
    const int &Id1 = qd0.vPt() [ 1 ]->Id();
    const int &Id2 = qd0.vPt() [ 2 ]->Id();
    const int &Id3 = qd0.vPt() [ 3 ]->Id();

#if 0
    if ( qd0.Id() == 44 ) {
      std :: cout << "Id0 = " << Id0 << " , "
                  << "Id1 = " << Id1 << " , "
                  << "Id2 = " << Id2 << " , "
                  << "Id3 = " << Id3 << "\n";
    }

#endif

    outputFile_ << std :: setw(1)
                << std :: fixed
                << "iQd =" << qd0.Id() << "; "
                << "u = [x(" << Id0 << ") x(" << Id1 << ") x(" << Id2 << ") x(" << Id3 << ") x(" << Id0 << ")]; "
                << "v = [y(" << Id0 << ") y(" << Id1 << ") y(" << Id2 << ") y(" << Id3 << ") y(" << Id0 << ")]; "
                << "plot(u, v, 'k'); ";

    if ( qd0.color() != Color :: Black ) {
      outputFile_ << "fill(u,v," << colorSymbol( qd0.color() ) << "); ";
    }

    if ( writeText ) {
      outputFile_ << "xbc(iQd) = ( " << qd0.pos() [ 0 ] << "); "
                  << "ybc(iQd) = ( " << qd0.pos() [ 1 ] << "); "
                  << "text(xbc(iQd), ybc(iQd), num2str(iQd));";
    }

    outputFile_ << std :: endl;
  }

  template< typename EntityT >
  void write(const msbGrid :: Block< EntityT > &bl0, const bool &writeText = false)
  {
    if ( bl0.size() == 0 ) {
      return;
    }

    for ( int iEnt = 0; iEnt < bl0.size(); ++iEnt ) {
      write(bl0.vEnt() [ iEnt ], writeText);
    }
  }
};

#endif  /* #ifndef _M_FILES_WRITER_HH */