Commit 21e7ad98 by Francois Gygi

fixed comments


git-svn-id: http://qboxcode.org/svn/qb/trunk@954 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 5751e9b0
......@@ -82,13 +82,12 @@ void CGCellStepper::compute_new_cell(double e0, const valarray<double>& sigma,
}
// convert descent direction dcell to cell space from the stress tensor
// dcell = - omega * sigma * A
// dcell = sigma * A
valarray<double> dcell(9); // descent direction in cell space
assert(sigma.size()==6);
// next line: local copy of sigma to circumvent compiler error
valarray<double> sigma_loc(sigma);
cell.smatmult3x3(&sigma_loc[0],cell.amat(),&dcell[0]);
// dcell *= -cell.volume();
for ( int i = 0; i < dcell.size(); i++ )
f0[3*nat_+i] = dcell[i];
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment