-
Notifications
You must be signed in to change notification settings - Fork 191
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Update hypre_ParVectorPrintIJ
#1033
Conversation
src/parcsr_mv/par_vector.c
Outdated
for (j = 0; j < local_size; j++) | ||
{ | ||
hypre_fprintf(file, "%b %.14e\n", | ||
(HYPRE_BigInt) j + base_j, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Does this need to be BigInt? I thought local sizes were always int?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
These are global indices, so need to be BigInt
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@victorapm is base_j supposed to be BigInt then? Index j only goes from 0 to local_size, so it seems not a BigInt to me. base_j is passed as argument as HYPRE_Int.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Yes, let me update this, thanks for catching!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Done in 200bbf9
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@victorapm I left you a few comments for your review.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks Victor!!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks fine.
Should we extend the ReadIJ function to be able to read vectors that have been printed with this extended function?
Yes, I have a separate PR for that :) |
* Add multi-component vector support to hypre_ParVectorPrintIJ * Make hypre_ParVectorPrintIJ conforming with HYPRE_IJVectorPrint * Make hypre_ParVectorReadIJ conforming with HYPRE_IJVectorRead
Add multi-component vector support to
hypre_ParVectorPrintIJ
This feature is important for GEOS