E INVALID BC FOUND in genfast
Hi all, I am getting the following error when I try to run nek : *0 3912 BC FAIL* *EXIT: E INVALID BC FOUND in genfast 3912* * * I grep-ed to the source code to get an idea of what the error actually means, but I couldn't comprehend whats going on in the code (fast3d.f) . Could anybody tell me what this error actually means ? I checked the Boundary condition of the particular element (3912) and everything seems to be OK. Thanks for any help. Regards Shriram Jagannathan
Shriram, It's indicating that one of the bcs for element 3912 is not valid. What are those? Paul On Fri, 25 Jun 2010, [email protected] wrote:
Hi all,
I am getting the following error when I try to run nek :
*0 3912 BC FAIL* *EXIT: E INVALID BC FOUND in genfast 3912*
* * I grep-ed to the source code to get an idea of what the error actually means, but I couldn't comprehend whats going on in the code (fast3d.f) . Could anybody tell me what this error actually means ? I checked the Boundary condition of the particular element (3912) and everything seems to be OK. Thanks for any help.
Regards Shriram Jagannathan
Hi Paul, I am also getting something similar, however it is referring to many elements with a "BC fail" and then gives the E INVALID BC FOUND. I have attached the case and output log for referrence. I checked the elements that it is referring to and the BC's are good, so I'm not sure why it tells me this. I looked in fast3d.f and also couldn't figure out what "ierr and ierrmx" refer to. Could explain what this error means or what it is checking for? Thanks for any help! - Michael ----- Original Message ----- From: [email protected] To: [email protected] Sent: Friday, June 25, 2010 1:39:43 PM GMT -06:00 US/Canada Central Subject: Re: [Nek5000-users] E INVALID BC FOUND in genfast Shriram, It's indicating that one of the bcs for element 3912 is not valid. What are those? Paul On Fri, 25 Jun 2010, [email protected] wrote:
Hi all,
I am getting the following error when I try to run nek :
*0 3912 BC FAIL* *EXIT: E INVALID BC FOUND in genfast 3912*
* * I grep-ed to the source code to get an idea of what the error actually means, but I couldn't comprehend whats going on in the code (fast3d.f) . Could anybody tell me what this error actually means ? I checked the Boundary condition of the particular element (3912) and everything seems to be OK. Thanks for any help.
Regards Shriram Jagannathan
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Michael, Shriram, . ierr is set to the element that has a problem . ierrmx is the maximum of ierr, taken over all processors. If more than one element has a problem, you won't know it because it will only flag the element having the largest number -- thus indicating that there is _a_ problem, but not the extent of the problem. (Flagging errors is tricky when you have millions of elements.) My general approach at this point would be something like the following. Suppose the troublesome element is 3192. Add the following code to the gen_fast routine: mid = gllnid(3192) ie = gllel (3192) if (mid.eq.nid) write(6,6) mid,ie,(cbc(k,ie,1),k=1,6) 6 format(2i8,6('.',a3,'.'),' BOUNDARY') when you then run the code and grep for BOUNDARY in the logfile you should see something like: 20 5 .v ..E ..E ..E ..E ..E . and not (say) 20 5 .v .. E .. E .. E .. E .. E . It's imperative that the BCs are 3 characters each and that the spaces are in the correct place. Only by printing the bcs can you see what bcs are really being passed to the code. Paul On Wed, 30 Jun 2010, [email protected] wrote:
Hi Paul,
I am also getting something similar, however it is referring to many elements with
a "BC fail" and then gives the E INVALID BC FOUND. I have attached the case
and output log for referrence.
I checked the elements that it is referring to and the BC's are good, so I'm not sure
why it tells me this. I looked in fast3d.f and also couldn't figure out what "ierr and ierrmx"
refer to. Could explain what this error means or what it is checking for? Thanks
for any help!
- Michael
----- Original Message ----- From: [email protected] To: [email protected] Sent: Friday, June 25, 2010 1:39:43 PM GMT -06:00 US/Canada Central Subject: Re: [Nek5000-users] E INVALID BC FOUND in genfast
Shriram,
It's indicating that one of the bcs for element 3912 is not valid.
What are those?
Paul
On Fri, 25 Jun 2010, [email protected] wrote:
Hi all,
I am getting the following error when I try to run nek :
*0 3912 BC FAIL* *EXIT: E INVALID BC FOUND in genfast 3912*
* * I grep-ed to the source code to get an idea of what the error actually means, but I couldn't comprehend whats going on in the code (fast3d.f) . Could anybody tell me what this error actually means ? I checked the Boundary condition of the particular element (3912) and everything seems to be OK. Thanks for any help.
Regards Shriram Jagannathan
_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Ah, thanks very much for the insight! I was thinking that it was either a formatting issue or perhaps a tolerance issue, but after looking at my rea I already see where it may have gone wrong. I'll add the below, and try again. Thanks again Paul! - Michael ----- Original Message ----- From: [email protected] To: [email protected] Sent: Wednesday, June 30, 2010 7:29:36 PM GMT -06:00 US/Canada Central Subject: Re: [Nek5000-users] E INVALID BC FOUND in genfast Michael, Shriram, . ierr is set to the element that has a problem . ierrmx is the maximum of ierr, taken over all processors. If more than one element has a problem, you won't know it because it will only flag the element having the largest number -- thus indicating that there is _a_ problem, but not the extent of the problem. (Flagging errors is tricky when you have millions of elements.) My general approach at this point would be something like the following. Suppose the troublesome element is 3192. Add the following code to the gen_fast routine: mid = gllnid(3192) ie = gllel (3192) if (mid.eq.nid) write(6,6) mid,ie,(cbc(k,ie,1),k=1,6) 6 format(2i8,6('.',a3,'.'),' BOUNDARY') when you then run the code and grep for BOUNDARY in the logfile you should see something like: 20 5 .v ..E ..E ..E ..E ..E . and not (say) 20 5 .v .. E .. E .. E .. E .. E . It's imperative that the BCs are 3 characters each and that the spaces are in the correct place. Only by printing the bcs can you see what bcs are really being passed to the code. Paul On Wed, 30 Jun 2010, [email protected] wrote: > > > Hi Paul, > > > > I am also getting something similar, however it is referring to many elements with > > a "BC fail" and then gives the E INVALID BC FOUND. I have attached the case > > and output log for referrence. > > > > I checked the elements that it is referring to and the BC's are good, so I'm not sure > > why it tells me this. I looked in fast3d.f and also couldn't figure out what "ierr and ierrmx" > > refer to. Could explain what this error means or what it is checking for? Thanks > > for any help! > > > > - Michael > > > ----- Original Message ----- > From: [email protected] > To: [email protected] > Sent: Friday, June 25, 2010 1:39:43 PM GMT -06:00 US/Canada Central > Subject: Re: [Nek5000-users] E INVALID BC FOUND in genfast > > > Shriram, > > It's indicating that one of the bcs for element 3912 is not valid. > > What are those? > > Paul > > > On Fri, 25 Jun 2010, [email protected] wrote: > >> Hi all, >> >> I am getting the following error when I try to run nek : >> >> *0 3912 BC FAIL* >> *EXIT: E INVALID BC FOUND in genfast 3912* >> >> * >> * >> I grep-ed to the source code to get an idea of what the error actually >> means, but I couldn't comprehend whats going on in the code (fast3d.f) . >> Could anybody tell me what this error actually means ? I checked the >> Boundary condition of the particular element (3912) and everything seems to >> be OK. >> Thanks for any help. >> >> Regards >> Shriram Jagannathan >> > _______________________________________________ > Nek5000-users mailing list > [email protected] > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
Paul, Thanks for the reply. Earlier that day you had mentioned about the format in connect2.f and I looked at that and implemented in the rea file and that took care of the error. As you had pointed out earlier, the b.c was not a 3-character string and that is why I had the problem. Thanks again. Regards Shriram
participants (1)
-
nek5000-users@lists.mcs.anl.gov