/[escript]/trunk/finley/src/Mesh_rec8.c
ViewVC logotype

Diff of /trunk/finley/src/Mesh_rec8.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 149 by jgs, Thu Sep 1 03:31:39 2005 UTC revision 150 by jgs, Thu Sep 15 03:44:45 2005 UTC
# Line 1  Line 1 
1    /*
2     ******************************************************************************
3     *                                                                            *
4     *       COPYRIGHT  ACcESS 2003,2004,2005 -  All Rights Reserved              *
5     *                                                                            *
6     * This software is the property of ACcESS. No part of this code              *
7     * may be copied in any form or by any means without the expressed written    *
8     * consent of ACcESS.  Copying, use or modification of this software          *
9     * by any unauthorised person is illegal unless that person has a software    *
10     * license agreement with ACcESS.                                             *
11     *                                                                            *
12     ******************************************************************************
13    */
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /*   Finley: generates rectangular meshes */  /*   Finley: generates rectangular meshes */
# Line 7  Line 21 
21    
22  /**************************************************************/  /**************************************************************/
23    
 /*   Copyrights by ACcESS Australia 2003 */  
24  /*   Author: gross@access.edu.au */  /*   Author: gross@access.edu.au */
25  /*   Version: $Id$ */  /*   Version: $Id$ */
26    
27  /**************************************************************/  /**************************************************************/
28    
 #include "Common.h"  
 #include "Finley.h"  
 #include "Mesh.h"  
29  #include "RectangularMesh.h"  #include "RectangularMesh.h"
30    
31  /**************************************************************/  /**************************************************************/
# Line 49  Finley_Mesh* Finley_RectangularMesh_Rec8 Line 59  Finley_Mesh* Finley_RectangularMesh_Rec8
59        
60    sprintf(name,"Rectangular %d x %d mesh",N0,N1);    sprintf(name,"Rectangular %d x %d mesh",N0,N1);
61    out=Finley_Mesh_alloc(name,2,order);    out=Finley_Mesh_alloc(name,2,order);
62    if (Finley_ErrorCode!=NO_ERROR) return NULL;    if (! Finley_noError()) return NULL;
63    
64    out->Elements=Finley_ElementFile_alloc(Rec8,out->order);    out->Elements=Finley_ElementFile_alloc(Rec8,out->order);
65    if (useElementsOnFace) {    if (useElementsOnFace) {
# Line 60  Finley_Mesh* Finley_RectangularMesh_Rec8 Line 70  Finley_Mesh* Finley_RectangularMesh_Rec8
70       out->ContactElements=Finley_ElementFile_alloc(Line3_Contact,out->order);       out->ContactElements=Finley_ElementFile_alloc(Line3_Contact,out->order);
71    }    }
72    out->Points=Finley_ElementFile_alloc(Point1,out->order);    out->Points=Finley_ElementFile_alloc(Point1,out->order);
73    if (Finley_ErrorCode!=NO_ERROR) {    if (! Finley_noError()) {
74         Finley_Mesh_dealloc(out);         Finley_Mesh_dealloc(out);
75         return NULL;         return NULL;
76    }    }
# Line 70  Finley_Mesh* Finley_RectangularMesh_Rec8 Line 80  Finley_Mesh* Finley_RectangularMesh_Rec8
80    Finley_NodeFile_allocTable(out->Nodes,N0*N1);    Finley_NodeFile_allocTable(out->Nodes,N0*N1);
81    Finley_ElementFile_allocTable(out->Elements,NE0*NE1);    Finley_ElementFile_allocTable(out->Elements,NE0*NE1);
82    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
83    if (Finley_ErrorCode!=NO_ERROR) {    if (! Finley_noError()) {
84        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
85        return NULL;        return NULL;
86    }    }
# Line 275  Finley_Mesh* Finley_RectangularMesh_Rec8 Line 285  Finley_Mesh* Finley_RectangularMesh_Rec8
285    printf("timing: mesh generation: %.4e sec\n",Finley_timer()-time0);    printf("timing: mesh generation: %.4e sec\n",Finley_timer()-time0);
286    #endif    #endif
287    
288    if (Finley_ErrorCode!=NO_ERROR) {    if (! Finley_noError()) {
289        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
290        return NULL;        return NULL;
291    }    }
# Line 283  Finley_Mesh* Finley_RectangularMesh_Rec8 Line 293  Finley_Mesh* Finley_RectangularMesh_Rec8
293  }  }
294    
295  /*  /*
 * $Log$  
296  * Revision 1.3  2005/09/01 03:31:36  jgs  * Revision 1.3  2005/09/01 03:31:36  jgs
297  * Merge of development branch dev-02 back to main trunk on 2005-09-01  * Merge of development branch dev-02 back to main trunk on 2005-09-01
298  *  *
299    * Revision 1.2.2.2  2005/09/07 06:26:19  gross
300    * the solver from finley are put into the standalone package paso now
301    *
302  * Revision 1.2.2.1  2005/08/24 02:02:18  gross  * Revision 1.2.2.1  2005/08/24 02:02:18  gross
303  * timing output switched off. solver output can be swiched through getSolution(verbose=True) now.  * timing output switched off. solver output can be swiched through getSolution(verbose=True) now.
304  *  *

Legend:
Removed from v.149  
changed lines
  Added in v.150

  ViewVC Help
Powered by ViewVC 1.1.26