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

fixed preconditioning #13

Merged
merged 1 commit into from
Oct 18, 2016
Merged

Conversation

barakovic
Copy link
Contributor

Fixed order of repetitions for each kernels.

Checked that the norm of the matrix A generated with 'doNormalizeKernels' = False

norm_A1 = np.linalg.norm(A,axis=0)`

is approximately equal to the norm calculate by COMMIT generate with 'doNormalizeKernels' = True

norm1 = np.repeat(mit.KERNELS['wmr_norm'],nF)
norm2 = np.repeat(mit.KERNELS['wmh_norm'],nE)
norm3 = np.repeat(mit.KERNELS['iso_norm'],nV)
norm_fib = np.kron(np.ones(mit.KERNELS['wmr'].shape[0]), mit.DICTIONARY['TRK']['norm'])

norm_A2 =  np.hstack( (norm1*norm_fib,norm2,norm3) )

norm_A1 ≈ norm_A2

@daducci daducci merged commit 4f8844d into daducci:master Oct 18, 2016
@daducci daducci mentioned this pull request Oct 26, 2016
@daducci daducci added the bug label Oct 26, 2016
This was referenced Oct 26, 2016
@barakovic barakovic deleted the fixed_preconditioning branch March 14, 2017 12:40
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants