CFD/hypre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed struct_matvec for case x=y

Browse Source
This commit is contained in:
Ulrike Yang 2015-08-17 17:14:51 -07:00
parent c28f3c9145
commit 5e19d62b87
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
struct_mv/struct_matvec.c
View File

@ -125,6 +125,8 @@ hypre_StructMatvecCompute( void *matvec_vdata,
HYPRE_Complex temp;
HYPRE_Int compute_i, i;
hypre_StructVector *x_tmp = NULL;
/*-----------------------------------------------------------------------
* Initialize some things
*-----------------------------------------------------------------------*/
@ -168,6 +170,11 @@ hypre_StructMatvecCompute( void *matvec_vdata,
return hypre_error_flag;
}
if (x == y)
{
x_tmp = hypre_StructVectorClone(y);
x = x_tmp;
}
/*-----------------------------------------------------------------------
* Do (alpha != 0.0) computation
*-----------------------------------------------------------------------*/
@ -275,6 +282,12 @@ hypre_StructMatvecCompute( void *matvec_vdata,
}
}
if (x_tmp)
{
hypre_StructVectorDestroy(x_tmp);
x = y;
}
return hypre_error_flag;
}