blob: 3a3cc0e54edccd74a1611cfab1c9991d01542d39 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
// Copyright (C) 2016 - IIT Bombay - FOSSEE
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
// Author: Siddhesh Wani
// Organization: FOSSEE, IIT Bombay
// Email: toolbox@scilab.in
//Function to extract lower triagular entries from given matrix
#include "triu.h"
void dtriua (double *in, int row, int column, double diag, double *out)
{
int row_counter, col_counter = 0;
for(row_counter=0; row_counter < row; row_counter++)
{
for(col_counter=0; col_counter < column; col_counter++)
{
if((double)(row_counter + diag) <= (double)col_counter)
{
out[col_counter*row+row_counter] = in[col_counter*row+row_counter];
}
else
{
out[col_counter*row+row_counter] = 0;
}
}
}
}
|
