<div dir="ltr">Hello all,<div><br></div><div>I am writing to to let you know my current progress on the B-Spline geometry interface and propose some of my doubts.</div><div><br></div><div>I modified the B-Spline basis function calculation.  Instead of building the whole pyramid to get the basis functions we need on the top level, we can build an inverted triangular to merely compute the non-zero functions. Later, we only need to multiple these non-vanishing values with the corresponding control points, which is more efficient. In addition, I implemented a simple working NURBS geometry interface and added some test code in the main file. I tested on my computer, they both look fine.</div><div><br></div><div>Now I'm working on other interface functions. I have implemented the following functions:</div><div><div>GeometryType <span style="white-space:pre"> </span>type () const;</div><div>bool <span class="" style="white-space:pre">   </span>affine () const;</div><div>int <span class="" style="white-space:pre">       </span>corners () const;</div><div>GlobalCoordinate <span class="" style="white-space:pre"> </span>corner (int i) const;</div><div>GlobalCoordinate <span class="" style="white-space:pre">     </span>center () const;</div><div><br></div></div><div>All the codes are available via <a href="https://gitlab.dune-project.org/Xinyun.Li/dune-iga">https://gitlab.dune-project.org/Xinyun.Li/dune-iga</a>. I will be very grateful if you can review my code as some of my thoughts may possibly be wrong.The next step is  Jacobian related function. I have written some pseudocode and have got a clue on it. So hopefully I can finish it this week or early next week.</div><div><br></div><div>My doubts </div><div>1. How to implement the following functions?</div><div>ctype volume() const.</div><div>LocalCoordinate local(const GlobalCoordinate&) const;</div><div>I don't know how to define the volume of the B-Spline structure and the inverse mapping from the physical space to the element space also seems non-trivial.</div><div>2. Matrix operation</div><div>I will need the transpose and pseudo-inverse of a matrix. I think these functions should already been implemented in DDUNE. I had a quick look of the FieldMatrix class but not sure if it is enough.</div><div><br></div><div>Thanks for your time in advance.</div><div><br></div><div>Best regards,</div><div>Xinyun </div></div>