/* Do that, you can't compile it with VC++ when <complex> included
#if ( defined(__BORLANDC__) || _MSC_VER ) && !defined( __GNUG__ ) 
inline float abs (float v) { return (float)fabs( v); } 
inline double abs (double v) { return fabs( v); } 
inline long double abs (long double v) { return fabsl( v); }
#endif
*/

#include <complex>
using namespace std;
#include "matrix.h"
 
typedef complex<double> Complex;
typedef math::matrix<Complex>  CMat;
 
int main()
{
    CMat m(2,2);
    m(0,0) = Complex(0,1);
    m(0,1) = m(1,0) = 0.0;
    m(1,1) = Complex(0,1);
 
    cout << m.Det() << endl;
 
    return 0;
}
// prints (-1,0)