Skip to content
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

Add second version of TensorProductMatrixSymmetricSum::apply_inverse() #14178

Merged
merged 1 commit into from Aug 23, 2022

Conversation

peterrum
Copy link
Member

@peterrum peterrum commented Aug 7, 2022

No description provided.

Copy link
Member

@drwells drwells left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Otherwise looks good

/**
* Same as above but the user can provide a user-owned temporal array,
* resolving the two issues described above.
*/
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It looks like the temporary array needs to be of length 2 * n: could you document that here?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The array is set to the right size internally.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Got it. Could we document that instead? We typically require that users set output arrays to the right sizes so we should clarify what to do for work arrays.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have added a comment!

Copy link
Member

@kronbichler kronbichler left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good, but there are a few typos in the text.

include/deal.II/lac/tensor_product_matrix.h Outdated Show resolved Hide resolved
include/deal.II/lac/tensor_product_matrix.h Outdated Show resolved Hide resolved
@peterrum
Copy link
Member Author

@kronbichler Updated!

@kronbichler
Copy link
Member

We run into timeouts for some tests affected by this:

The following tests FAILED:
	4171 - lac/tensor_product_matrix_04.debug (Timeout)
	4200 - lac/tensor_product_matrix_vectorized_02.debug (Timeout)
	4218 - lac/tensor_product_matrix_02.debug (Timeout)
	4249 - lac/tensor_product_matrix_vectorized_04.debug (Timeout)
	4259 - lac/tensor_product_matrix_05.debug (Timeout)
	4261 - lac/tensor_product_matrix_vectorized_03.debug (Timeout)
	4290 - lac/tensor_product_matrix_06.debug (Timeout)
	4322 - lac/tensor_product_matrix_vectorized_01.debug (Timeout)
	4334 - lac/tensor_product_matrix_03.debug (Timeout)
	4338 - lac/tensor_product_matrix_01.debug (Timeout)

Can you take a look?

@peterrum
Copy link
Member Author

@kronbichler Should be fixed now. I din't move the mutex but duplicated it. So that the code dead locked.

Copy link
Member

@kronbichler kronbichler left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Good.

@kronbichler kronbichler merged commit 2fa8d12 into dealii:master Aug 23, 2022
mkghadban pushed a commit to OpenFCST/dealii that referenced this pull request Sep 8, 2022
Add second version of TensorProductMatrixSymmetricSum::apply_inverse()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants