IRELIEF-SIGMA(1) User Commands IRELIEF-SIGMA(1)NAME
irelief-sigma - Command line interface to irelief-sigma in mlpy (version 2.2.0)
SYNOPSIS
irelief-sigma [options]
OPTIONS -h, --help
show this help message and exit
-d FILE, --data=FILE
data - required
-t THETA, --theta=THETA
theta (default 0.001)
-s, --standardize
standardize data
-n, --normalize
normalize data
-m MIN, --min=MIN
min value (default -5)
-M MAX, --max=MAX
max value (default 5)
-p STEPS, --steps=STEPS
steps (default 11)
-e SCALE, --scale=SCALE
scale (lin or log, default log)
mlpy July 2010 IRELIEF-SIGMA(1)
Check Out this Related Man Page
CLAR1V(l) ) CLAR1V(l)
NAME
CLAR1V - compute the (scaled) r-th column of the inverse of the sumbmatrix in rows B1 through BN of the tridiagonal matrix L D L^T - sigma
I
SYNOPSIS
SUBROUTINE CLAR1V( N, B1, BN, SIGMA, D, L, LD, LLD, GERSCH, Z, ZTZ, MINGMA, R, ISUPPZ, WORK )
INTEGER B1, BN, N, R
REAL MINGMA, SIGMA, ZTZ
INTEGER ISUPPZ( * )
REAL D( * ), GERSCH( * ), L( * ), LD( * ), LLD( * ), WORK( * )
COMPLEX Z( * )
PURPOSE
CLAR1V computes the (scaled) r-th column of the inverse of the sumbmatrix in rows B1 through BN of the tridiagonal matrix L D L^T - sigma
I. The following steps accomplish this computation : (a) Stationary qd transform, L D L^T - sigma I = L(+) D(+) L(+)^T, (b) Progressive qd
transform, L D L^T - sigma I = U(-) D(-) U(-)^T, (c) Computation of the diagonal elements of the inverse of
L D L^T - sigma I by combining the above transforms, and choosing
r as the index where the diagonal of the inverse is (one of the)
largest in magnitude.
(d) Computation of the (scaled) r-th column of the inverse using the
twisted factorization obtained by combining the top part of the
the stationary and the bottom part of the progressive transform.
ARGUMENTS
N (input) INTEGER
The order of the matrix L D L^T.
B1 (input) INTEGER
First index of the submatrix of L D L^T.
BN (input) INTEGER
Last index of the submatrix of L D L^T.
SIGMA (input) REAL
The shift. Initially, when R = 0, SIGMA should be a good approximation to an eigenvalue of L D L^T.
L (input) REAL array, dimension (N-1)
The (n-1) subdiagonal elements of the unit bidiagonal matrix L, in elements 1 to N-1.
D (input) REAL array, dimension (N)
The n diagonal elements of the diagonal matrix D.
LD (input) REAL array, dimension (N-1)
The n-1 elements L(i)*D(i).
LLD (input) REAL array, dimension (N-1)
The n-1 elements L(i)*L(i)*D(i).
GERSCH (input) REAL array, dimension (2*N)
The n Gerschgorin intervals. These are used to restrict the initial search for R, when R is input as 0.
Z (output) COMPLEX array, dimension (N)
The (scaled) r-th column of the inverse. Z(R) is returned to be 1.
ZTZ (output) REAL
The square of the norm of Z.
MINGMA (output) REAL
The reciprocal of the largest (in magnitude) diagonal element of the inverse of L D L^T - sigma I.
R (input/output) INTEGER
Initially, R should be input to be 0 and is then output as the index where the diagonal element of the inverse is largest in mag-
nitude. In later iterations, this same value of R should be input.
ISUPPZ (output) INTEGER array, dimension (2)
The support of the vector in Z, i.e., the vector Z is nonzero only in elements ISUPPZ(1) through ISUPPZ( 2 ).
WORK (workspace) REAL array, dimension (4*N)
FURTHER DETAILS
Based on contributions by
Inderjit Dhillon, IBM Almaden, USA
Osni Marques, LBNL/NERSC, USA
Ken Stanley, Computer Science Division, University of
California at Berkeley, USA
LAPACK version 3.0 15 June 2000 CLAR1V(l)
I have a huge matrix file containing some 1.5 million rows and 6000 columns. The matrix looks something like this:
1 2 3
4 5 6
7 8 9
3 4 5
I want to add all the numbers in the columns of this matrix and display the result to my stdout. This means that the numbers in the first column are:
... (2 Replies)
Hi,
Just trying to get to grips with sed and awk for some reporting for work and I need some assistance:
I have a file that lists policy names on the first line and then on the second line whether the policy is active or not.
Policy Name: Policy1
Active: yes
Policy... (8 Replies)
version info :
vi availabe with RHEL 5.4
I have a text file with 10,000 lines. I want to copy lines from 5000th line to 7000th and redirect to a file. Any idea how I can do this?
Note:
The above scenario is just an example. In my actual requirement, the file has 14 million lines and I want... (9 Replies)
Hi everyone,
I know the following questions are noobish questions but I am asking them because I am confused about the basics of history behind UNIX and LINUX.
Ok onto business, my questions are-:
Was/Is UNIX ever an open source operating system ?
If UNIX was... (21 Replies)
Hello,
I couldn't find an actual introduction thread, so I decided to just put this here.
I go by d0wngrade online. I have been programming in multiple languages for about 15+ years. I started with standard web design languages like HTML and CSS, but I then advanced from design to development... (2 Replies)
Hi guys...
The first active code line in AudioScope.sh is set -u .
This causes a complete exit if a variable is used/found but has not been allocated at the start of the program.
However, apart from writing code to do the task, is there a switch to to check which variables have been... (17 Replies)
Hi.
In thread https://www.unix.com/shell-programming-and-scripting/267833-grouping-counting.html rovf and I had a mini-discussion on grep and awk.
Here is a demo script that compares the awk and grep approaches for this single problem:
#!/usr/bin/env bash
# @(#) s2 Demonstrate group... (1 Reply)
Hello,
I have to fish out some specific columns from a file based on the header value. I have the list of columns I need in a different file. I thought I could read in the list of headers I need,
# file with header names of required columns in required order
headers_file=$2
# read contents... (11 Replies)
For those interested in installing dash shell on OSX Lion to help test POSIX compliancy of shell scripts, it is quite easy. I did it like this:
If you don't have gcc on your system:
0. Download and install the Command Line Tools for Xcode package from Sign In - Apple *
1. Download the dash... (2 Replies)
Hello and thanks in advance for any help anyone can offer me
I'm trying to learn the find command and thought I was understanding it... Apparently I was wrong. I was doing compound searches and I started getting weird results with the -size test. I was trying to do a search on a 1G file owned by... (14 Replies)
I have data of an excel files as given below,
file1
org1_1 1 1 2.5 100
org1_2 1 2 5.5 98
org1_3 1 3 7.2 88
file2
org2_1 1 1 2.5 100
org2_2 1 2 5.5 56
org2_3 1 3 7.2 70
I have multiple excel files as above shown.
I have to copy column 1, column 4 and paste into a new excel file as... (26 Replies)
Dear All,
Taking a break from Vue.js coding for the site, SEO and YT videos; and hopefully addressing some well deserved criticism from some here that I have been too focused on the visual aspects of the forums versus the substance and the community....
While the "current generation... (9 Replies)
Hi all...
Well guys and gals, I jumped in at the deep end and found things that PERL cannot do by default.
Many tricky terminal escape codes are not catered for so I had to create workarounds.
One thing I searched for was this:
Passing perl variable to shell command
AND, @Neo this was... (15 Replies)