summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/matrixOperations/diag/ddiagina.c
blob: c4093b2dc38f44a2e8f2d277e154a9e6e45c162a (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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
/*
 * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 * Copyright (C) 2008-2008 - INRIA - Allan SIMON
 *
 * 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
 *
 */

#include "diag.h"
#include<stdlib.h>  // Used for Absolute value of insert_post
void ddiagina(double *in, int _row,int _column,int insert_post,double *out)
{
    
     int i, j;
     if(_row == 1)
     {
       _column = _column;

     }
     else
     {
      _column = _row;

     }

      for(i=0;i < ((_column+abs(insert_post))*(_column+abs(insert_post)));i++)
      {
         out[i] = 0; 
  
      }
    
       if(insert_post < 0)
	{
      	for ( j = 0 ; j < _column ; j++ )
         {
                
                 
             
               out[((_column+abs(insert_post))*j)+abs(insert_post)+j] = in[j];   
          }

	}
    else
	{
      	for ( j = insert_post ; j < _column+insert_post ; j++ )
         {
                
                 
            out[((_column+insert_post)*j)-insert_post+j] = in[j-insert_post] ;  
                   
          }
    	}

}