Hi,<br> <br>Using units and ublas for new code is great, but how about using both in legacy code? I’ve got some legacy numerical code using primitive types, i plan to use units for some reasons, i.e. to check performance using both primitive types and boost::units types, and to check my code at compile time which is its whole purpose after all. I spent a day and my code it’s already broken, specially because of the initializations and ublas, some self-explanatory code follows:<br>
 <br> <br>----- calculations.hpp -----------<br>//typedefs to get library-free code<br> <br>#include &lt;complex&gt;<br>#include &lt;iostream&gt;<br> <br>#if USING_BOOST_UNITS<br> <br> <br>#include &lt;boost/typeof/std/complex.hpp&gt; <br>
#include &lt;boost/units/systems/si/power.hpp&gt;<br>//...<br>using namespace boost::units;<br>using namespace boost::units::si;<br>using namespace std;<br> <br>typedef std::complex&lt;double&gt; complex_type;<br> <br>typedef quantity&lt;electric_potential,complex_type&gt; tPotencial;<br>
typedef quantity&lt;current,complex_type&gt; tIntensidad;<br>typedef quantity&lt;resistance,complex_type&gt; tResistancia;<br> <br>//...<br>#elif<br> <br>typedef double tPotencial;<br>typedef double tIntensidad;<br>typedef double tResistancia;<br>
 <br>//...<br>#endif<br>------------------------------------------<br> <br>-----  some_legacy_calculations.cpp ------<br> <br>#define USING_BOOST_UNITS <br>tPotencial tension = complex_type(12.5,0.0)*volts; //ok<br>tPotencial tension2 = complex_type(12.5,0.0); //auch!<br>
 <br>some_file &lt;&lt; tesion + tension2 &lt;&lt; std::endl;<br>some_file.close();<br> <br>//get Real part from &#39;tension&#39; and pretend it&#39;s a current<br>quantity&lt;current, double&gt; phony_tension = tension.value().real * amperes;//compiler doesn&#39;t complain, but how do I prevent this to compile? <br>
 <br>#undef USING_BOOST_UNITS<br>tPotencial tension = complex_type(12.5,0.0)*volts; //auch!<br>tPotencial tension2 = complex_type(12.5,0.0); //ok<br> <br>some_file.open(&quot;/*some_file*/&quot;); <br>some_file &gt;&gt; tesion + tension2 &gt;&gt; std::endl; //auch!<br>
 <br> <br>#define USING_BOOST_UNITS <br>//use ublas<br>ublas::matrix&lt;tAdmitancia, ublas::column_major&gt; admi(2, 2);      <br>ublas::matrix&lt;tPotencia, ublas::column_major&gt; ind(2, 1);   <br> <br>ind(0, 0) = complex_type (1, 0)*volts;<br>
ind(1, 0) = complex_type (-1, 0)*volts;<br> <br>admi(0, 0) = complex_type (0, -2)*siemens;<br>admi(0, 1) = complex_type (0, 2)*siemens;<br>admi(1, 0) = complex_type (0, 2)*siemens;<br>admi(1, 1) = complex_type (0, -2)*siemens;<br>
 <br> <br>int er = lapack::gesv(admi, ind);//auch!<br> <br>------------------------------------------<br> <br>For each of those problems i found some workaround but I&#39;d like to discuss whether it&#39;s ok or not.<br>Thanks for  your time.