Skip to content

Commit 442d638

Browse files
CHG: better error message to avoid people using periodic BC with non linear elements
1 parent 53f2f40 commit 442d638

File tree

1 file changed

+3
-1
lines changed

1 file changed

+3
-1
lines changed

src/c/modules/ModelProcessorx/CreateNodes.cpp

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -604,7 +604,9 @@ void CreateNodes(Nodes* nodes, IoModel* iomodel,int analysis,int finite_element,
604604
int *vertex_pairing = NULL;
605605
int numvertex_pairing;
606606
iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.stressbalance.vertex_pairing");
607-
_assert_(numvertex_pairing==0 || finite_element==P1Enum);
607+
if(numvertex_pairing>0){
608+
if(finite_element!=P1Enum) _error_("Periodic boundary conditions only implemented for linear elements");
609+
}
608610
for(int i=0;i<numvertex_pairing;i++){
609611
int nid1 = vertex_pairing[2*i+0]-1;
610612
int nid2 = vertex_pairing[2*i+1]-1;

0 commit comments

Comments
 (0)