//************************************************************************
// --- Logiciel MEF++
// ---
// --- Copyright 1996-2000 GIREF, Université Laval
// --- TOUS DROITS RÉSERVÉS
// --- ALL RIGHTS RESERVED
// ---
// --- Ce logiciel est couvert par la Loi sur le droit d'auteur.
// --- L'utilisation, ou la modification de ce logiciel sous toutes ses
// --- formes, que ce soit code source ou code compilé, à des fins
// --- personnelles et non commerciales sont autorisées sans frais pour
// --- autant que la présente notice de copyright ainsi que cette
// --- permission apparaissent dans toutes les copies ainsi que dans la
// --- documentation.
// --- Le GIREF et l'Université Laval ne prétendent en aucune façon que
// --- ce logiciel convient à un emploi quelconque. Celui-ci est
// --- distribué sans aucune garantie implicite ou explicite.
//************************************************************************

//********************************************************************
//  Nom du fichier: conversionVU.cc
//
//********************************************************************
//  Description: Programme servant a faire la conversion du
//  format GIREF vers VU.
//
//********************************************************************

#include "ChaineCar.h"
#include "ChampGeoLin.h"
#include "ChampScalConstElem.h"
#include "ChampScalLin.h"
#include "ChampScalQuad.h"
#include "ChampTensO2SymLin.h"
#include "ChampVect2DLin.h"
#include "ChampVect2DQuad.h"
#include "ChampVect3DLin.h"
#include "ChampVect3DQuad.h"
#include "ExportVU.h"
#include "ExpressionAlgebrique.h"
#include "fonctionsChaineCar.h"
#include "InitialisationArgvArgc.h"
#include "InitialisationTraitementSignal.h"
#include "Maillage.h"
#include "PETScInitialisation.h"
#include "traducteurERMsg.h"
#include <cstring>
#include <ctime>
#include <fstream>


using namespace std;

ExportVU lExport; // pas de compression a priori car on va ajouter le temps dans le fichier

template <class PTTypeChamp>
void importeChamp(PTTypeChamp& pChamp,std::ifstream& pFichier,ChaineCar& pNomDuChamp, ChaineCar& pNomDuFichier){
  pChamp.asgnNom(pNomDuChamp);
  ERMsg lMsg = pChamp.importe(pFichier);
  if(lMsg) lExport.ajouteChamp(pChamp, pNomDuChamp);
//   if(lMsg) lMsg = lExport.exporteInfo(pNomDuFichier);
}

int main(int argc, char *argv[]) {

  ///

  //***********************************
  //* BLOC: Initialisations générales *
  //***********************************
  PETScInitialisation lPETScInit(&argc, &argv);
  sInitialisationArgvArgc.asgnArguments(argc, argv);
  initialiseTraitementSignal();

  const vector<string> lArgv = sInitialisationArgvArgc.reqArgv();

  // On déclare une variable qui contiendra potentiellement un message d'erreur
  ERMsg lMsg = ERMsg(ERMsg::OK);

  // On recupere le GroupeProcessus associe a PETSC_COMM_WORLD
  const PAGroupeProcessus lGroupeGlobal = PETScInitialisation::reqPetscCommWorld();

  // On déclare les variables qui seront extraites de la ligne de commande
  ChaineCar lNom;
  ChaineCar lTypeDeChamp1;
  ChaineCar lNomDuFichier1;
  ChaineCar lNomDuChamp1;
  std::ifstream lFichier1;
  ChaineCar lTypeDeChamp2;
  ChaineCar lNomDuFichier2;
  ChaineCar lNomDuChamp2;
  std::ifstream lFichier2;
  Entier lNombre(0);

  // Vérification des paramètres passés à la ligne de commande.
  const Entier lNbArg = lArgv.size();
  if (lNbArg < 3 || lNbArg > 4) {
    ChaineCar lNomProgramme = sInitialisationArgvArgc.reqNomProgramme();
    ChaineCar lMessageUsage = "USAGE: \n";
    lMessageUsage += enleveCheminAccesFichier(lNomProgramme);
    lMessageUsage += " Prefixe_fichier_mail nombre_fichier";
    lMsg.ajoute(ERMsg(ERMsg::ERREUR, lMessageUsage));
  }
  else {
    lNom           = lArgv[1];
    lNomDuFichier1 = lArgv[1];
    lNomDuChamp1   = "M";
    lNomDuFichier2 = lArgv[1];
    lNomDuChamp2   = "U";
    lNombre = atoi(argv[2]);
  }

  // Lecture du Maillage
  Maillage lMail;
  if(lMsg) lMsg = lMail.importe(lNom);

  for (int i = 1; i<lNombre+1;++i){
    ChaineCar lNom1(lNomDuFichier1);
    lNom1=ajouteNumeroALaFin(lNomDuFichier1,i,3);
    ChaineCar lNom2(lNom1);
    lNom1+= ".M.champscallin";
    lNom2 += ".U.champv3dlin";

    if(lMsg) {
      lFichier1.open(lNom1.c_str());
      if(lFichier1){
        importeString(lFichier1,lTypeDeChamp1);
        // on rewind
        lFichier1.seekg (0, std::ios::beg);
      } else {
        lMsg.ajoute(ERMsg(ERMsg::ERREUR,"ERREUR LECTURE FICHIER: "+ lNom1));
      }
    }
    if(lMsg) {
      lFichier2.open(lNom2.c_str());
      if(lFichier2){
        importeString(lFichier2,lTypeDeChamp2);
        // on rewind
        lFichier2.seekg (0, std::ios::beg);
      } else {
        lMsg.ajoute(ERMsg(ERMsg::ERREUR,"ERREUR LECTURE FICHIER: "+ lNom2));
      }
    }

    if(lMsg){
      // Creation d'un champ geometrique LINEAIRE
      ChampGeoLin lChampGeo(lMail);

      lExport.asgnChampGeometrique(lChampGeo);
      lExport.asgnDegreInterpolation(Exportation::MAILLAGE_LINEAIRE);

      ChampScalLin lChamp1(lChampGeo);
      lChamp1.asgnNom(lNomDuChamp1);
      lChamp1.importe(lFichier1);
      lExport.ajouteChamp(lChamp1,lNomDuChamp1);

      ChampVect3DLin lChamp2(lChampGeo);
      lChamp2.asgnNom(lNomDuChamp2);
      lChamp2.importe(lFichier2);
      lExport.ajouteChamp(lChamp2,lNomDuChamp2);

      if(lMsg) lMsg = lExport.exporteInfo(lNom1);

      lFichier1.close();
      lFichier2.close();
    } // if (lMsg)
  }
  if (!lMsg) {
    cerr << traducteurERMsg(lMsg);
  }

  return 0;
}