-
Notifications
You must be signed in to change notification settings - Fork 16
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
[WIP] PETSc #165
base: master
Are you sure you want to change the base?
[WIP] PETSc #165
Conversation
we also need to add petsc to travis... |
8486dfc
to
6f1aa38
Compare
I got problem with the test:
I am going to work on other stuff and then I will fix it. |
@nicola-giuliani any idea?
|
|
se provo ad aggiungere compress:
|
compress dove lo hai aggiunto? |
a occhio mi pare che quando synchronize chiama il costruttore copia il vettore originale non sia stato compressed. Bisognerebbe capire l'ultima operazione fatta sul vettore e chiamare un compress prima di syncronize. Cosi spero che l'errore cambi...come a Pozioni nel primo semestre ;) |
Ho provato:
ma ho errori legati alle funzioni const e a solution che è const. se provo a togliere const è un macello. |
(lldb) run
Process 23235 launched: './tests/poisson_04.debug/poisson_04.debug' (x86_64)
Number of active cells: 256 (on 5 levels)
Number of degrees of freedom: 289(289)
computing consistent initial conditions with the option use_y_diff please be patient.
compute initial conditions: done.
0 ----> 0.01
0.01 ----> 0.02 5
0.02 ----> 0.03 6
################ restart #########
max_kelly > threshold
0.0100929 > 0.01
######################################
Number of active cells: 274 (on 6 levels)
Number of degrees of freedom: 331(331)
computing consistent initial conditions with the option use_y_dot please be patient.
compute initial conditions: done.
0.03 ----> 0.04
0.04 ----> 0.05 15
0.05 ----> 0.06 19
################ restart #########
max_kelly > threshold
0.0101164 > 0.01
######################################
Number of active cells: 328 (on 6 levels)
Number of degrees of freedom: 392(392)
computing consistent initial conditions with the option use_y_dot please be patient.
compute initial conditions: done.
################ restart #########
max_kelly > threshold
0.0107257 > 0.01
######################################
Number of active cells: 394 (on 6 levels)
Number of degrees of freedom: 462(462)
computing consistent initial conditions with the option use_y_dot please be patient.
compute initial conditions: done.
0.06 ----> 0.07
0.07 ----> 0.08 22
0.08 ----> 0.09 19
0.09 ----> 0.1
################ restart #########
max_kelly > threshold
0.0102897 > 0.01
######################################
Number of active cells: 442 (on 6 levels)
Number of degrees of freedom: 507(507)
computing consistent initial conditions with the option use_y_dot please be patient.
compute initial conditions: done.
0.1 ----> 0.11
iterations: 5
cells dofs u_L2 u_H1
442 507 7.964e-02 - 5.706e-01 -
*** The MPI_Attr_get() function was called after MPI_FINALIZE was invoked.
*** This is disallowed by the MPI standard.
*** Your MPI job will now abort.
[(null):23235] Local abort after MPI_FINALIZE completed successfully; not able to aggregate error messages, and not able to guarantee that all other processes were killed!
Process 23235 exited with status = 1 (0x00000001)
(lldb) up
error: invalid thread |
se il poisson_04.cc è quello che c'è in questa PR è normale che salti alla fine.. oppure l'errore non è legato al ciclo for? |
per curiosità.. puoi mettere |
|
Una cosa importante.... PETSc NON è multithreaded. Qui bisogna sempre forzare il numero di threads a uno. Se fai così funziona? |
Require mathLab/deal2lkit#292
I got some problem with reinit of LinearOperator when the LinearOperator is associated to a PETSc matrix.
I am trying to make a test similar to linear_operator_04 for PETSc. -> too messy. I will look for another approach.