geometrie/src/fem_solution.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       fem_solution.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include <math.h>
#include <iostream>
#include <stdlib.h>
#include <string.h>

#include "fem_solution.h"
#include "fem_maillage.h"
#include "ot_mathematique.h"


#include <stdio.h> // pour la fonction remove() qui remplace DeleteFile


FEM_SOLUTION::FEM_SOLUTION(FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite):MG_IDENTIFICATEUR(),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite)
{
    nom_fichier=new char[strlen(chemin)+2];
    legende=new std::string[nb];
    for (int i=0;i<nb;i++)
    {
        char mess[15];
        sprintf(mess,"Champs_%d",i);
        legende[i]=mess ;
    }
    strcpy(nom_fichier,chemin);
    if (code!=0)
    {
        in=fopen(nom_fichier,"wb");
        for (int i=0;i<code;i++)
            for (int j=0;j<nb;j++)
            {
                double val=0.0;
                fwrite(&val,sizeof(double),1,in);
            }
        fclose(in);
    }
    in=fopen(nom_fichier,"r+b");
}

FEM_SOLUTION::FEM_SOLUTION(unsigned long num,FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite):MG_IDENTIFICATEUR(num),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite)
{
    nom_fichier=new char[strlen(chemin)+2];
    legende=new std::string[nb];
    for (int i=0;i<nb;i++)
    {
        char mess[15];
        sprintf(mess,"Champs_%d",i);
        legende[i]=mess ;
    }
    strcpy(nom_fichier,chemin);
    if (code!=0)
    {
        in=fopen(nom_fichier,"wb");
        for (int i=0;i<code;i++)
            for (int j=0;j<nb;j++)
            {
                double val=0.0;
                fwrite(&val,sizeof(double),1,in);
            }
        fclose(in);
    }
    in=fopen(nom_fichier,"r+b");
}

FEM_SOLUTION::~FEM_SOLUTION()
{
    delete [] nom_fichier;
    delete [] legende;
    if (in!=NULL) fclose(in);
}

void FEM_SOLUTION::efface(void)
{
    fclose(in);
    remove(nom_fichier);
// DeleteFile(nom_fichier);
    in=NULL;
}

int FEM_SOLUTION::get_type_solution(void)
{
    return typeentite;
}

double FEM_SOLUTION::lire(int i,int j)
{
    long pos=(nb_champs*i+j)*sizeof(double);
    fseek(in,pos,SEEK_SET);
    double val;
    size_t res=fread(&val,sizeof(double),1,in);
    return (val);
}

void FEM_SOLUTION::ecrire(int i,int j,double val)
{
    long pos=(nb_champs*i+j)*sizeof(double);
    fseek(in,pos,SEEK_SET);
    fwrite(&val,sizeof(double),1,in);
    fflush(in);
}


void FEM_SOLUTION::change_legende(int num,std::string val)
{
    legende[num]=val;
}

std::string FEM_SOLUTION::get_legende(int num)
{
    return legende[num];
}

int FEM_SOLUTION::get_nb_champ(void)
{
    return nb_champs;
}

FEM_MAILLAGE* FEM_SOLUTION::get_maillage(void)
{
    return femmai;
}

double FEM_SOLUTION::get_legende_min(void)
{
    return solmin;
}

double FEM_SOLUTION::get_legende_max(void)
{
    return solmax;
}

void FEM_SOLUTION::active_solution(int num)
{
    int nb;
    solmax=-1e308;
    solmin=1e308;
    int i=0;
    if (typeentite==ENTITE_NOEUD)
    {
        LISTE_FEM_NOEUD::iterator it;
        for (FEM_NOEUD* no=femmai->get_premier_noeud(it);no!=NULL;no=femmai->get_suivant_noeud(it))
        {
            double val=lire(i,num);
            no->change_solution(val);
            if (val<solmin) solmin=val;
            if (val>solmax) solmax=val;
            i++;
        }
    }
    if (typeentite==ENTITE_ELEMENT1)
    {
        LISTE_FEM_ELEMENT1::iterator it;
        for (FEM_ELEMENT1* seg=femmai->get_premier_element1(it);seg!=NULL;seg=femmai->get_suivant_element1(it))
        {
            double val=lire(i,num);
            seg->change_solution(val);
            if (val<solmin) solmin=val;
            if (val>solmax) solmax=val;
            i++;
        }
    }
    if (typeentite==ENTITE_ELEMENT2)
    {
        LISTE_FEM_ELEMENT2::iterator it;
        for (FEM_ELEMENT2* tri=femmai->get_premier_element2(it);tri!=NULL;tri=femmai->get_suivant_element2(it))
        {
            double val=lire(i,num);
            tri->change_solution(val);
            if (val<solmin) solmin=val;
            if (val>solmax) solmax=val;
            i++;
        }
    }
    if (typeentite==ENTITE_ELEMENT3)
    {
        LISTE_FEM_ELEMENT3::iterator it;
        for (FEM_ELEMENT3* tet=femmai->get_premier_element3(it);tet!=NULL;tet=femmai->get_suivant_element3(it))
        {
            double val=lire(i,num);
            tet->change_solution(val);
            if (val<solmin) solmin=val;
            if (val>solmax) solmax=val;
            i++;
        }
    }
    if (OPERATEUR::egal(fabs(solmin-solmax),0.,1e-10)==1)
        solmax=solmin+1e-10;
}

std::string FEM_SOLUTION::get_nom(void)
{
    return nomsolution;
}

void FEM_SOLUTION::enregistrer(std::ostream& o)
{
    o << "%" << get_id() << "=FEM_SOLUTION(" << typeentite << "," << nomsolution << ",$" << femmai->get_id() << "," << nb_champs << "," << nom_fichier << ",(";
    for (int i=0;i<nb_champs;i++)
    {
        o << legende[i];
        if (i!=nb_champs-1) o<<",";
        else o<<")";
    }


    o << ");" << std::endl;
}


std::ostream& operator << (std::ostream& o,FEM_SOLUTION& sol)
{
    sol.enregistrer(o);
    return o;
}


Revision:	309
Committed:	Tue Jan 24 15:25:16 2012 UTC (13 years, 3 months ago) by francois
File size:	6181 byte(s)
Log Message:	Preparation du maillage FEM pour inclure de nouveaux élements non necessairement triangle/tetra
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// fem_solution.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include <math.h>
27			#include <iostream>
28			#include <stdlib.h>
29			#include <string.h>
30
31			#include "fem_solution.h"
32			#include "fem_maillage.h"
33			#include "ot_mathematique.h"
34
35
36			#include <stdio.h> // pour la fonction remove() qui remplace DeleteFile
37
38
39
40
41
42
43			FEM_SOLUTION::FEM_SOLUTION(FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite):MG_IDENTIFICATEUR(),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite)
44			{
45			nom_fichier=new char[strlen(chemin)+2];
46			legende=new std::string[nb];
47			for (int i=0;i<nb;i++)
48			{
49			char mess[15];
50			sprintf(mess,"Champs_%d",i);
51			legende[i]=mess ;
52			}
53			strcpy(nom_fichier,chemin);
54			if (code!=0)
55			{
56			in=fopen(nom_fichier,"wb");
57			for (int i=0;i<code;i++)
58			for (int j=0;j<nb;j++)
59			{
60			double val=0.0;
61			fwrite(&val,sizeof(double),1,in);
62			}
63			fclose(in);
64			}
65			in=fopen(nom_fichier,"r+b");
66			}
67
68			FEM_SOLUTION::FEM_SOLUTION(unsigned long num,FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite):MG_IDENTIFICATEUR(num),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite)
69			{
70			nom_fichier=new char[strlen(chemin)+2];
71			legende=new std::string[nb];
72			for (int i=0;i<nb;i++)
73			{
74			char mess[15];
75			sprintf(mess,"Champs_%d",i);
76			legende[i]=mess ;
77			}
78			strcpy(nom_fichier,chemin);
79			if (code!=0)
80			{
81			in=fopen(nom_fichier,"wb");
82			for (int i=0;i<code;i++)
83			for (int j=0;j<nb;j++)
84			{
85			double val=0.0;
86			fwrite(&val,sizeof(double),1,in);
87			}
88			fclose(in);
89			}
90			in=fopen(nom_fichier,"r+b");
91			}
92
93			FEM_SOLUTION::~FEM_SOLUTION()
94			{
95			delete [] nom_fichier;
96			delete [] legende;
97			if (in!=NULL) fclose(in);
98			}
99
100			void FEM_SOLUTION::efface(void)
101			{
102			fclose(in);
103			remove(nom_fichier);
104			// DeleteFile(nom_fichier);
105			in=NULL;
106			}
107
108			int FEM_SOLUTION::get_type_solution(void)
109			{
110			return typeentite;
111			}
112
113			double FEM_SOLUTION::lire(int i,int j)
114			{
115			long pos=(nb_champsi+j)sizeof(double);
116			fseek(in,pos,SEEK_SET);
117			double val;
118			size_t res=fread(&val,sizeof(double),1,in);
119			return (val);
120			}
121
122			void FEM_SOLUTION::ecrire(int i,int j,double val)
123			{
124			long pos=(nb_champsi+j)sizeof(double);
125			fseek(in,pos,SEEK_SET);
126			fwrite(&val,sizeof(double),1,in);
127			fflush(in);
128			}
129
130
131			void FEM_SOLUTION::change_legende(int num,std::string val)
132			{
133			legende[num]=val;
134			}
135
136			std::string FEM_SOLUTION::get_legende(int num)
137			{
138			return legende[num];
139			}
140
141			int FEM_SOLUTION::get_nb_champ(void)
142			{
143			return nb_champs;
144			}
145
146			FEM_MAILLAGE* FEM_SOLUTION::get_maillage(void)
147			{
148			return femmai;
149			}
150
151			double FEM_SOLUTION::get_legende_min(void)
152			{
153			return solmin;
154			}
155
156			double FEM_SOLUTION::get_legende_max(void)
157			{
158			return solmax;
159			}
160
161			void FEM_SOLUTION::active_solution(int num)
162			{
163			int nb;
164			solmax=-1e308;
165			solmin=1e308;
166			int i=0;
167			if (typeentite==ENTITE_NOEUD)
168			{
169			LISTE_FEM_NOEUD::iterator it;
170			for (FEM_NOEUD* no=femmai->get_premier_noeud(it);no!=NULL;no=femmai->get_suivant_noeud(it))
171			{
172			double val=lire(i,num);
173			no->change_solution(val);
174			if (val<solmin) solmin=val;
175			if (val>solmax) solmax=val;
176			i++;
177			}
178			}
179	francois	309	if (typeentite==ENTITE_ELEMENT1)
180	francois	283	{
181	francois	309	LISTE_FEM_ELEMENT1::iterator it;
182			for (FEM_ELEMENT1* seg=femmai->get_premier_element1(it);seg!=NULL;seg=femmai->get_suivant_element1(it))
183	francois	283	{
184			double val=lire(i,num);
185			seg->change_solution(val);
186			if (val<solmin) solmin=val;
187			if (val>solmax) solmax=val;
188			i++;
189			}
190			}
191	francois	309	if (typeentite==ENTITE_ELEMENT2)
192	francois	283	{
193	francois	309	LISTE_FEM_ELEMENT2::iterator it;
194			for (FEM_ELEMENT2* tri=femmai->get_premier_element2(it);tri!=NULL;tri=femmai->get_suivant_element2(it))
195	francois	283	{
196			double val=lire(i,num);
197			tri->change_solution(val);
198			if (val<solmin) solmin=val;
199			if (val>solmax) solmax=val;
200			i++;
201			}
202			}
203	francois	309	if (typeentite==ENTITE_ELEMENT3)
204	francois	283	{
205	francois	309	LISTE_FEM_ELEMENT3::iterator it;
206			for (FEM_ELEMENT3* tet=femmai->get_premier_element3(it);tet!=NULL;tet=femmai->get_suivant_element3(it))
207	francois	283	{
208			double val=lire(i,num);
209			tet->change_solution(val);
210			if (val<solmin) solmin=val;
211			if (val>solmax) solmax=val;
212			i++;
213			}
214			}
215			if (OPERATEUR::egal(fabs(solmin-solmax),0.,1e-10)==1)
216			solmax=solmin+1e-10;
217			}
218
219			std::string FEM_SOLUTION::get_nom(void)
220			{
221			return nomsolution;
222			}
223
224			void FEM_SOLUTION::enregistrer(std::ostream& o)
225			{
226			o << "%" << get_id() << "=FEM_SOLUTION(" << typeentite << "," << nomsolution << ",$" << femmai->get_id() << "," << nb_champs << "," << nom_fichier << ",(";
227			for (int i=0;i<nb_champs;i++)
228			{
229			o << legende[i];
230			if (i!=nb_champs-1) o<<",";
231			else o<<")";
232			}
233
234
235			o << ");" << std::endl;
236			}
237
238
239
240
241			std::ostream& operator << (std::ostream& o,FEM_SOLUTION& sol)
242			{
243			sol.enregistrer(o);
244			return o;
245			}
246
247