/[escript]/trunk-mpi-branch/finley/src/CPPAdapter/MeshAdapterFactory.cpp
ViewVC logotype

Diff of /trunk-mpi-branch/finley/src/CPPAdapter/MeshAdapterFactory.cpp

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

revision 1222 by ksteube, Tue May 15 03:23:17 2007 UTC revision 1223 by gross, Fri Aug 3 02:40:39 2007 UTC
# Line 97  namespace finley { Line 97  namespace finley {
97              int periodic2,              int periodic2,
98              int integrationOrder,              int integrationOrder,
99                      int reducedIntegrationOrder,                      int reducedIntegrationOrder,
100              int useElementsOnFace)              int useElementsOnFace,
101                int useFullElementOrder)
102    {    {
103  //     cout << "n0=" << n0 << " n1=" << n1 << " n2=" << n2  //     cout << "n0=" << n0 << " n1=" << n1 << " n2=" << n2
104  //   << " order=" << order  //   << " order=" << order
# Line 118  namespace finley { Line 119  namespace finley {
119    
120      if (order==1) {      if (order==1) {
121        fMesh=Finley_RectangularMesh_Hex8(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,        fMesh=Finley_RectangularMesh_Hex8(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,
122                      useElementsOnFace) ;                      useElementsOnFace,useFullElementOrder) ;
123      }      }
124          else if (order==2) {          else if (order==2) {
125        fMesh=Finley_RectangularMesh_Hex20(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,        fMesh=Finley_RectangularMesh_Hex20(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,
126                       useElementsOnFace) ;                       useElementsOnFace,useFullElementOrder) ;
127      } else {      } else {
128        stringstream temp;        stringstream temp;
129        temp << "Illegal interpolation order: " << order;        temp << "Illegal interpolation order: " << order;
# Line 139  namespace finley { Line 140  namespace finley {
140              int periodic0,int periodic1,              int periodic0,int periodic1,
141              int integrationOrder,              int integrationOrder,
142                          int reducedIntegrationOrder,                          int reducedIntegrationOrder,
143              int useElementsOnFace)              int useElementsOnFace,
144                    int useFullElementOrder)
145    {    {
146      int numElements[]={n0,n1};      int numElements[]={n0,n1};
147      double length[]={l0,l1};      double length[]={l0,l1};
# Line 148  namespace finley { Line 150  namespace finley {
150      Finley_Mesh* fMesh=0;      Finley_Mesh* fMesh=0;
151      if (order==1) {      if (order==1) {
152        fMesh=Finley_RectangularMesh_Rec4(numElements, length,periodic,integrationOrder,reducedIntegrationOrder,        fMesh=Finley_RectangularMesh_Rec4(numElements, length,periodic,integrationOrder,reducedIntegrationOrder,
153                      useElementsOnFace);                      useElementsOnFace,useFullElementOrder);
154      }      }
155      else if (order==2) {      else if (order==2) {
156        fMesh=Finley_RectangularMesh_Rec8(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,        fMesh=Finley_RectangularMesh_Rec8(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,
157                      useElementsOnFace);                      useElementsOnFace,useFullElementOrder);
158      }      }
159      else {      else {
160        stringstream temp;        stringstream temp;
161        temp << "Illegal interpolation order: " << order;        temp << "Illegal interpolation order: " << order;
       setFinleyError(VALUE_ERROR,temp.str().c_str());  
     }  
     //  
     // Convert any finley errors into a C++ exception  
     checkFinleyError();  
     AbstractContinuousDomain* temp=new MeshAdapter(fMesh);  
     return temp;  
   }  
   AbstractContinuousDomain*  interval(int n0,int order,double l0,int periodic0,  
                int integrationOrder,  
                        int reducedIntegrationOrder,  
                int useElementsOnFace)  
   {  
     int numElements[]={n0};  
     double length[]={l0};  
     int periodic[]={periodic0};  
     Finley_Mesh* fMesh;  
     if (order==1) {  
       fMesh=Finley_RectangularMesh_Line2(numElements, length,periodic,integrationOrder,reducedIntegrationOrder,  
                      useElementsOnFace);  
     }  
     else if (order==2) {  
       fMesh=Finley_RectangularMesh_Line3(numElements,length,periodic,integrationOrder,reducedIntegrationOrder,  
                      useElementsOnFace);  
     }  
     else {  
       stringstream temp;  
       temp << "Illegal interpolation order: " << order;  
162        setFinleyError(VALUE_ERROR,temp.str().c_str());        setFinleyError(VALUE_ERROR,temp.str().c_str());
163      }      }
164      //      //

Legend:
Removed from v.1222  
changed lines
  Added in v.1223

  ViewVC Help
Powered by ViewVC 1.1.26