Some bugs about knot insert

[Brisk]

Dear all,

I found that when the inserted knots is too large, FEAP8.5 cannot obtain the correct control points. For example, when I insert 200 knot in a rectangular mesh, the coordinate of the last control point is incorrect. And the mesh obtained by FEAP cannot be recognized in analysis.

Attached is the test file, could you help me?

[Prof. R.L. Taylor]

We did not anticipate that users would use igafeap with knot vectors as large as 200.  There are still some fixed dimension arrays that are used to store the knot vectors during generation.  These are located in the include file "cnurb.h" and are named UU1, UU2 (store knot vector values) and PP1, PP2 which store the control point data (which in the Piegl algorithm include the weight and thus get dimensioned PP1(ndm+1,0:lastknot) and similar for PP2.  Thus if you have a 2-d geometry you run out of space at 200 knots. 

These need to be changed and then be sure to recompile the entire iga part of the code.  If you use SHOW COMMon options you need also to recompile the 'outary.f' file in the program directory (this is the only routine that references cnurb.h).

[Brisk]

Following your advice, I have fixed this problem. Really appreciate for your help!