Dear Steve & Lammps user

I want to calculate r_ij(f_ij*v_i) as a output qunatity. But it seems I am doing some thing wrong.

the code is as follows:

if (rsq < cutsq[itype][jtype]) {

force->pair->single(i,j,itype,jtype,rsq,factor_coul,factor_lj,0,one);

fforce = one.fforce;

//////////////////////

// multiplying f_ij*v_i

stress[i][0] -= v[i][0]*fforce;

stress[i][1] -= v[i][0]*fforce;

stress[i][2] -= v[i][0]*fforce;

stress[i][3] -= v[i][0]*fforce;

stress[i][4] -= v[i][0]*fforce;

stress[i][5] -= v[i][0]*fforce;

//calculating r_ij (f_ij*v_i)

for (i = 0; i < nlocal; i++) {

if (rsq < cutsq[itype][jtype]) {

stress[i][0] *= delx*delx;

stress[i][1] *= dely*dely;

stress[i][2] *= delz*delz;

stress[i][3] *= delx*dely;

stress[i][4] *= delx*delz;

stress[i][5] *= dely*delz;}

}

Is something wrong with the code? because I am not getting the result as I expected. any suggestions???

thanx and regards

Vijay