|
EpetraExt
Development
|
00001 subroutine izero ( n, x, incx ) 00002 00003 c ================================================================== 00004 c ================================================================== 00005 c ==== izero -- initialize integer vector to zero ==== 00006 c ================================================================== 00007 c ================================================================== 00008 00009 c purpose ... initializes integer vector to zero 00010 00011 c created ... mar. 8, 1985 00012 c last modified ... apr. 19, 1985 00013 00014 c ================================================================== 00015 00016 c -------------- 00017 c ... parameters 00018 c -------------- 00019 00020 integer n, incx 00021 00022 integer x (*) 00023 00024 c ------------------- 00025 c ... local variables 00026 c ------------------- 00027 00028 integer xaddr, i 00029 00030 c ================================================================== 00031 00032 if ( incx .eq. 1 ) then 00033 00034 c ---------------------------------- 00035 c ... unit increment (standard case) 00036 c ---------------------------------- 00037 00038 do 100 i = 1, n 00039 x(i) = 0 00040 100 continue 00041 00042 else 00043 00044 c ---------------------- 00045 c ... non-unit increment 00046 c ---------------------- 00047 00048 xaddr = 1 00049 if ( incx .lt. 0 ) then 00050 xaddr = (-n+1)*incx + 1 00051 endif 00052 00053 do 200 i = 1, n 00054 x (xaddr) = 0 00055 xaddr = xaddr + incx 00056 200 continue 00057 00058 endif 00059 00060 return 00061 00062 end
1.7.6.1