/[escript]/trunk/finley/src/finleyC/ElementFile_markNodes.c
ViewVC logotype

Annotation of /trunk/finley/src/finleyC/ElementFile_markNodes.c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 150 - (hide annotations)
Thu Sep 15 03:44:45 2005 UTC (13 years, 11 months ago) by jgs
Original Path: trunk/esys2/finley/src/finleyC/ElementFile_markNodes.c
File MIME type: text/plain
File size: 3499 byte(s)
Merge of development branch dev-02 back to main trunk on 2005-09-15

1 jgs 150 /*
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 jgs 82 /**************************************************************/
16    
17     /* Finley: ElementFile */
18    
19     /* mark the used nodes with offeset: */
20    
21     /**************************************************************/
22    
23 jgs 150 /* Copyrights by ACcESS Australia 2003,2004,2005 */
24     /* Author: gross@access.edu.au */
25     /* Version: $Id$ */
26 jgs 82
27     /**************************************************************/
28    
29     #include "ElementFile.h"
30    
31     /**************************************************************/
32    
33 jgs 123 void Finley_ElementFile_markNodes(index_t* mask,index_t offset,Finley_ElementFile* in,bool_t useLinear) {
34     dim_t i,NN,NN2,e;
35     index_t color,*lin_node;
36 jgs 82 if (in!=NULL) {
37 jgs 123 index_t id[in->ReferenceElement->Type->numNodes];
38 jgs 82 for (i=0;i<in->ReferenceElement->Type->numNodes;i++) id[i]=i;
39     if (useLinear) {
40     NN=in->LinearReferenceElement->Type->numNodes;
41     lin_node=in->ReferenceElement->Type->linearNodes;
42     } else {
43     NN=in->ReferenceElement->Type->numNodes;
44     lin_node=id;
45     }
46     NN2=in->ReferenceElement->Type->numNodes;
47 jgs 126 if ((in->maxColor-in->minColor+1)*NN<in->numElements) {
48     #pragma omp parallel private(color)
49     {
50     for (color=in->minColor;color<=in->maxColor;color++) {
51     #pragma omp for private(e,i) schedule(static)
52     for (e=0;e<in->numElements;e++) {
53     if (in->Color[e]==color) {
54     for (i=0;i<NN;i++) mask[in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset]=1;
55     }
56     }
57     }
58     #pragma omp barrier
59 jgs 82 }
60 jgs 126 } else {
61     #pragma omp parallel for private(e,i) schedule(static)
62     for (e=0;e<in->numElements;e++) {
63     for (i=0;i<NN;i++) mask[in->Nodes[INDEX2(lin_node[i],e,NN2)]-offset]=1;
64     }
65     }
66 jgs 82 }
67     }
68     /*
69     * $Log$
70 jgs 150 * Revision 1.4 2005/09/15 03:44:22 jgs
71     * Merge of development branch dev-02 back to main trunk on 2005-09-15
72     *
73     * Revision 1.3.2.1 2005/09/07 06:26:18 gross
74     * the solver from finley are put into the standalone package paso now
75     *
76 jgs 126 * Revision 1.3 2005/07/22 03:53:08 jgs
77     * Merge of development branch back to main trunk on 2005-07-22
78     *
79 jgs 123 * Revision 1.2 2005/07/08 04:07:50 jgs
80     * Merge of development branch back to main trunk on 2005-07-08
81 jgs 82 *
82 jgs 126 * Revision 1.1.1.1.2.2 2005/07/18 10:34:54 gross
83     * some informance improvements when reading meshes
84     *
85 jgs 123 * Revision 1.1.1.1.2.1 2005/06/29 02:34:49 gross
86     * some changes towards 64 integers in finley
87     *
88     * Revision 1.1.1.1 2004/10/26 06:53:57 jgs
89     * initial import of project esys2
90     *
91 jgs 82 * Revision 1.1.1.1 2004/06/24 04:00:40 johng
92     * Initial version of eys using boost-python.
93     *
94     *
95     */

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26