<div dir="ltr"><div><div><div>Thanks for your suggestion. I think I would have to stick to <br>Eigen interface at this moment. However, is it possible to <br>add such provision as it is done for Eigen in CppAD?<br></div></div></div><br>From web search, I find some links where numpy routines <br>could be extended.<br><br><a href="https://docs.scipy.org/doc/numpy/reference/arrays.dtypes.html">https://docs.scipy.org/doc/numpy/reference/arrays.dtypes.html</a><br><a href="http://www.scipy-lectures.org/advanced/advanced_numpy/#universal-functions">http://www.scipy-lectures.org/advanced/advanced_numpy/#universal-functions</a><br></div>