Combining the application of DD schemes and NURBS bases with isogeometric analysis (IGA) [1]}, [2]} for high-fidelity modeling of nonlinear PDEs [3]}, [4]}, [5]} have enjoyed recent success at scale. However, many of these implementations lack full support to handle multiple geometric patches in a distributed memory setting due to non-trivial requirements on continuity constraints at patch boundaries.
Directly imposing higher-order geometric continuity in IGA requires specialized parameterizations in order to preserve the approximation properties [6]}, which can be difficult to parallelize [7]} generally.
In a similar vein, using B-spline bases to compute the MFA in parallel, while maintaining higher-order continuity across subdomains has not been fully explored previously.
