234 |
/* all elements have been read, now we have to identify the elements for finley */ |
/* all elements have been read, now we have to identify the elements for finley */ |
235 |
|
|
236 |
if (Finley_noError()) { |
if (Finley_noError()) { |
237 |
/* first we have to identify the elements to define Elementis and FaceElements */ |
/* first we have to identify the elements to define Elementis and FaceElements */ |
238 |
|
if (final_element_type == NoType) { |
239 |
|
if (numDim==1) { |
240 |
|
final_element_type=Line2; |
241 |
|
} else if (numDim==2) { |
242 |
|
final_element_type=Tri3; |
243 |
|
} else if (numDim==2) { |
244 |
|
final_element_type=Tet4; |
245 |
|
} |
246 |
|
} |
247 |
|
if (final_face_element_type == NoType) { |
248 |
|
if (numDim==1) { |
249 |
|
final_face_element_type=Point1; |
250 |
|
} else if (numDim==2) { |
251 |
|
final_face_element_type=Line2; |
252 |
|
} else if (numDim==2) { |
253 |
|
final_face_element_type=Tri3; |
254 |
|
} |
255 |
|
} |
256 |
mesh_p->Elements=Finley_ElementFile_alloc(final_element_type,mesh_p->order); |
mesh_p->Elements=Finley_ElementFile_alloc(final_element_type,mesh_p->order); |
257 |
mesh_p->FaceElements=Finley_ElementFile_alloc(final_face_element_type,mesh_p->order); |
mesh_p->FaceElements=Finley_ElementFile_alloc(final_face_element_type,mesh_p->order); |
258 |
mesh_p->ContactElements=Finley_ElementFile_alloc(Point1_Contact,mesh_p->order); |
mesh_p->ContactElements=Finley_ElementFile_alloc(Point1_Contact,mesh_p->order); |