`#!/bin/tcsh`

set n = 128

set beta = `echo "(-0.6)*($n)*($n)" | bc -l`

set ppn = 16

set np = 31

set nodes = `echo "if ($np%$ppn==0) {$np/$ppn;} else {$np/$ppn+1}" | bc`

### One simple way to do simple floating point arithmetics in shell

Published September 15, 2014 linux Leave a CommentI have for long ignored blogging, but today I came across the useful linux command “paste” and I thought I should write it down. In the past I have used quite a few smart tricks to handle data or output files, but I failed to document these tricks. I really should do it this time.

So here is the story. I have an output file test.o1551739 which contains lines such as “iter = 125″ scattering everywhere. I need to get the sum of all the numbers in these lines. How can I achieve this without writing a program? The following is a little thought process.

First, I am so used to `grep`

that I can grab all such lines by doing:

`grep iter test.o1551739`

Then, use `sed`

to get rid of the suffix “iter = “:

`sed s/"iter = "//`

Now that I have a bunch of numbers, one in a line. How can I sum them up? Yah, here comes the usage of `paste`

. The command `paste`

allows me to merge the lines and put a delimiter in between. I chose to use the plus sign, because then it gives me an arithmetic expression. Try this:

`paste -sd+`

Once I have the arithmetic expression, I’ll just use `bc`

to calculate the result. So piping all these steps, here is the one line command:

`grep iter test.o1551739 | sed s/"iter = "// | paste -sd+ | bc`

Nice. I realize that there must be tons of ways to do the same job. But I really love my solution.

The `cases`

environment defined by the amsmath package sets each case to be inline math style. This will look ugly when the formulas in the cases are complicated, such as fractions or integrals. Consider the following sample code:

`\usepackage{amsmath}`

\[

f(x)=

\begin{cases}

\frac{1}{2(x-1)}, & x>1 \\

\frac{\Gamma(x)}{2^{x-1}}, & 0<x<1

\end{cases}

\]

The result looks like this:

The solution is to use the `dcases`

environment provided by the mathtools package instead. The prefix `d’ means `display’. It will set the cases in displayed math style (exactly the same as `\frac`

versus `\dfrac`

). So if we write the following code

`\usepackage{amsmath,mathtools}`

\[

f(x)=

\begin{dcases}

\frac{1}{2(x-1)}, & x>1 \\

\frac{\Gamma(x)}{2^{x-1}}, & 0<x<1

\end{dcases}

\]

then we will see the result which looks like this:

http://melusine.eu.org/syracuse/pstricks/

Many fantastic PSTricks examples.

A colleague of mine mentioned to me today the Bocher’s formula for computing the coefficients of the characteristic polynomial of a matrix. It seems that this formula does not appear too often in textbooks or literature. I’ll just write down the formula and the idea of a simple proof here.

Let the characteristic polynomial of a matrix be

Then the coefficients can be computed by

To prove the formula, note that the coefficient is the summation of all possible products of j eigenvalues, i.e.,

where denotes the j-combination of numbers from 1 to n, and the trace of is the sum of the th power of the eigenvalues, i.e.,

In addition, we have

The above indicates that the first part of cancels the second part of , whereas the second part of cancels the first part of The rest of proof becomes obvious now.

## Recent Comments