| Line 175... |
Line 175... |
PetscValidScalarPointer(p,4);
|
PetscValidScalarPointer(p,4);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ip->innerproducts++;
|
ip->innerproducts++;
|
if (ip->matrix) {
|
if (ip->matrix) {
|
ierr = IPApplyMatrix(ip,x,ip->work);CHKERRQ(ierr);
|
ierr = IPApplyMatrix_Private(ip,x);CHKERRQ(ierr);
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDot(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecDot(ip->Bx,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecTDot(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecTDot(ip->Bx,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDot(x,y,p);CHKERRQ(ierr);
|
ierr = VecDot(x,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
| Line 225... |
Line 225... |
PetscValidScalarPointer(p,4);
|
PetscValidScalarPointer(p,4);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ip->innerproducts++;
|
ip->innerproducts++;
|
if (ip->matrix) {
|
if (ip->matrix) {
|
ierr = IPApplyMatrix(ip,x,ip->work);CHKERRQ(ierr);
|
ierr = IPApplyMatrix_Private(ip,x);CHKERRQ(ierr);
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDotBegin(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecDotBegin(ip->Bx,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecTDotBegin(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecTDotBegin(ip->Bx,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDotBegin(x,y,p);CHKERRQ(ierr);
|
ierr = VecDotBegin(x,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
| Line 277... |
Line 277... |
PetscValidScalarPointer(p,4);
|
PetscValidScalarPointer(p,4);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
if (ip->matrix) {
|
if (ip->matrix) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDotEnd(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecDotEnd(ip->Bx,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecTDotEnd(ip->work,y,p);CHKERRQ(ierr);
|
ierr = VecTDotEnd(ip->Bx,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecDotEnd(x,y,p);CHKERRQ(ierr);
|
ierr = VecDotEnd(x,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
| Line 334... |
Line 334... |
PetscValidScalarPointer(p,5);
|
PetscValidScalarPointer(p,5);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ip->innerproducts += n;
|
ip->innerproducts += n;
|
if (ip->matrix) {
|
if (ip->matrix) {
|
ierr = IPApplyMatrix(ip,x,ip->work);CHKERRQ(ierr);
|
ierr = IPApplyMatrix_Private(ip,x);CHKERRQ(ierr);
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDot(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDot(ip->Bx,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecMTDot(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMTDot(ip->Bx,n,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDot(x,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDot(x,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
| Line 386... |
Line 386... |
PetscValidScalarPointer(p,5);
|
PetscValidScalarPointer(p,5);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ip->innerproducts += n;
|
ip->innerproducts += n;
|
if (ip->matrix) {
|
if (ip->matrix) {
|
ierr = IPApplyMatrix(ip,x,ip->work);CHKERRQ(ierr);
|
ierr = IPApplyMatrix_Private(ip,x);CHKERRQ(ierr);
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDotBegin(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDotBegin(ip->Bx,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecMTDotBegin(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMTDotBegin(ip->Bx,n,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDotBegin(x,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDotBegin(x,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
| Line 440... |
Line 440... |
PetscValidScalarPointer(p,5);
|
PetscValidScalarPointer(p,5);
|
|
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
ierr = PetscLogEventBegin(IP_InnerProduct,ip,x,0,0);CHKERRQ(ierr);
|
if (ip->matrix) {
|
if (ip->matrix) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDotEnd(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDotEnd(ip->Bx,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|
ierr = VecMTDotEnd(ip->work,n,y,p);CHKERRQ(ierr);
|
ierr = VecMTDotEnd(ip->Bx,n,y,p);CHKERRQ(ierr);
|
}
|
}
|
} else {
|
} else {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
if (ip->bilinear_form == IPINNER_HERMITIAN) {
|
ierr = VecMDotEnd(x,n,y,p);CHKERRQ(ierr);
|
ierr = VecMDotEnd(x,n,y,p);CHKERRQ(ierr);
|
} else {
|
} else {
|