Revision 1261

trunk/src/linad99/dmat42.cpp (revision 1261)
28 28

  
29 29
static const int  maxiter = 40;
30 30

  
31
/**
32
 * Description not yet available.
33
 * \param
34
 */
31
/*
35 32
static double pythag(double a, double b)
36 33
{
37 34
  double fa=fabs(a);
......
41 38
  else
42 39
    return fb*sqrt(1.0+square(fa/fb));
43 40
}
41
*/
44 42

  
45 43
/*
46 44
class sing_val_decomp
trunk/src/linad99/dmat3.cpp (revision 1261)
82 82
*/
83 83
dmatrix inv(const dmatrix& m1,const double& _ln_det, const int& _sgn)
84 84
{
85
  double d;
85
  double d = 0.0;
86 86
  double& ln_det=(double&)(_ln_det);
87 87
  ln_det=0.0;
88 88
  int& sgn=(int&)(_sgn);
......
430 430
*/
431 431
double det(const dmatrix& m1)
432 432
{
433
  double d;
433
  double d = 0.0;
434 434
  dmatrix a(m1.rowmin(),m1.rowmax(),m1.rowmin(),m1.rowmax());
435 435

  
436 436
  if (m1.rowmin()!=m1.colmin()||m1.rowmax()!=m1.colmax())
......
469 469
*/
470 470
double ln_det(const dmatrix& m1, const int& _sgn)
471 471
{
472
  double d;
472
  double d = 0.0;
473 473
  int& sgn=(int&)_sgn;
474 474
  dmatrix a(m1.rowmin(),m1.rowmax(),m1.rowmin(),m1.rowmax());
475 475

  

Also available in: Unified diff