LCPQ
/
DEHam
mirror of https://github.com/LCPQ/DEHam
2
0
Fork 0
Code Issues Releases Wiki Activity

removed redundant printing

This commit is contained in:
vijay gopal chilkuri 2016-12-21 17:42:56 +01:00
parent 9a7a10c37e
commit 260daba4b2
1 changed files with 2 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/ex1.c
View File

@ -243,29 +243,16 @@ int main(int argc,char **argv)
ierr = VecGetValues(xr, nlocal, indxr, valxr);CHKERRQ(ierr);
for (ii=Istart; ii<Iend; ii+=1) {
// PetscPrintf(PETSC_COMM_WORLD," Element # = %d Value = %18f \n", ii, valxr[ii-Istart]);
iii = ii+1;
xmat = 0.0;
getdet_(&iii, ideter);
/*
if(mpiid==1 && i==0){
for(kk=0; kk<=7; kk++){
// PetscPrintf(PETSC_COMM_WORLD," %d ", ideter[kk]);
printf(" %d ", ideter[kk]);
}
// PetscPrintf(PETSC_COMM_WORLD,"id= %d\n",mpiid);
printf("id= %d\n",mpiid);
}
*/
if(1){
norm=norm+valxr[ii]*valxr[ii];
for(kko=0;kko<=3;kko++){
for(kok=kko;kok<=3;kok++){
if(kok == kko && ideter[kok] != 3){
xmat=xmat+(3.0/4.0)*(valxr[ii]*valxr[ii]);
// if(!i)PetscPrintf(PETSC_COMM_WORLD,"i = %d ii = %d iaa2 = %d xmat = %18f xymat = %18f\n", i, ii, iaa2, xmat, xymat);
}
else{
if(ideter[kko] == 1 && ideter[kok] == 1){
@ -299,7 +286,7 @@ int main(int argc,char **argv)
}
}
}
if(!i && mpiid == 1)printf("id = %d ii = %d xmat = %18f xymat = %18f valxr = %18f \n", mpiid, ii-Istart, xmat, xymat, valxr[ii-Istart]);
for(kko=4;kko<=7;kko++){
for(kok=kko;kok<=7;kok++){
if(kok == kko && ideter[kok] != 3){
@ -337,7 +324,7 @@ int main(int argc,char **argv)
}
}
}
if(!i && mpiid == 1)printf("id = %d ii = %d xmat = %18f xymat = %18f\n", mpiid, ii-Istart, xmat, xymat);
for(kko=0;kko<=3;kko++){
for(kok=4;kok<=7;kok++){
if(kok == kko && ideter[kok] != 3){
@ -359,7 +346,6 @@ int main(int argc,char **argv)
ideter2[kok]=1;
adr_(ideter2, &iaa2);
iaa2 = iaa2 - 1;
if(!(i) && mpiid == 1)printf("iaa2 = %d\n",iaa2);
xmat=xmat+valxr[ii]*valxr[iaa2];
}
if(ideter[kko] == 2 && ideter[kok] == 1){
@ -371,7 +357,6 @@ int main(int argc,char **argv)
ideter2[kok]=2;
adr_(ideter2, &iaa2);
iaa2 = iaa2 - 1;
if(!(i) && mpiid == 1)printf("2 iaa2 = %d\n",iaa2);
xmat=xmat+valxr[ii]*valxr[iaa2];
}
}
@ -379,7 +364,6 @@ int main(int argc,char **argv)
}
//---------------------------------------
if(!i && mpiid == 1)printf("id = %d ii = %d xmat = %18f xymat = %18f\n", mpiid, ii-Istart, xmat, xymat);
xymat=xymat+xmat;
}
}