#2485: include complex and real scalar evolution code from Canuda in ET
Reporter: Roland Haas
Status: open
Milestone:
Version:
Type: task
Priority: major
Component:
Comment (by Giuseppe Ficarra):
ScalarEvolve:
> In current version, in Scalar\_calc\_Tmunu, "compute\_fluxes==1" and "use\_jacobian=false" is not available....Should we add the comments in param file ?
Do you mean “use\_jacobian = true“? If you have a look at line 254 of Scalar\_calc\_Tmunu.F90, you can see that there is an if statement that uses the combination "compute\_fluxes==1" and "use\_jacobian=false".
> in schedule, flux grid function is allocated, and Scalar\_zero\_densities is called when compute\_fluxes is true. But, in Scalar\_zero\_densities, flux grid functions are initialized only when use\_jacobian is false. So, if compute\_fluxes is true and use\_jacobian is true, flux grid function is not initialized.
fixed in [https://bitbucket.org/canuda/scalar/commits/cb03ec3b4d1bd2e32787ffbf2d32f4d…
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2485/include-complex-a…
#2629: MoL_PseudoEvolution vs ANALYSIS
Reporter: Gabriele Bozzola
Status: new
Milestone:
Version:
Type: bug
Priority: major
Component:
Comment (by Gabriele Bozzola):
Upon further reflection, and reading more carefully the previous message, my impression is that there’s no easy way to implement a `compute_every` function. I originally implemented the function following `WeylScal4`, which simply takes the modulo of the current iteration and the `compute_every` parameter. But, if I understand what is going on, this is not the correct thing to do, as I can see in my test runs \(unless `compute_every` is 1\).
So, should `compute_every` be removed from the `canuda` codes and `WeylScal4` ? \(Or at least, forced to be 1 for backwards compatibility\)
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2629/mol_pseudoevoluti…
#2629: MoL_PseudoEvolution vs ANALYSIS
Reporter: Gabriele Bozzola
Status: new
Milestone:
Version:
Type: bug
Priority: major
Component:
Comment (by Gabriele Bozzola):
It runs, but I am not sure it produces the correct results, so the PRs shouldn’t be merged at the moment.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2629/mol_pseudoevoluti…
#2629: MoL_PseudoEvolution vs ANALYSIS
Reporter: Gabriele Bozzola
Status: new
Milestone:
Version:
Type: bug
Priority: major
Component:
Comment (by Roland Haas):
In the call I said I knew how to do the scheduling in MoL\_PseudoEvolution. Turns out that this is not true. Here’s my current thoughts:
```text
rl=0 delta_t=4.0 timefac=1
rl=1 delta_t=2.0 timefac=2
rl=2 delta_t=1.0 timefac=4
it == 0:
rl=0 t_of_level=0.0
rl=1 t_of_level=0.0
rl=2 t_of_level=0.0
it == 1:
rl=0 t_of_level=4.0
rl=1 t_of_level=2.0
rl=2 t_of_level=1.0
it == 2:
rl=2 t_of_level=2.0
it == 3:
rl=1 t_of_level=4.0
rl=2 t_of_level=3.0
it == 4:
rl=2 t_of_level=4.0
if wanting to compute every Delta_t = 4.0 which is every 4 iterations, need to
evaluate at:
it = 0 for rl=0,1,2
and
it = 1 for rl=0
it = 3 for rl=1
it = 4 for rl=2
```
and I am not sure to to come up with an expression that involves only `iteration`, `timefac` and the `every` value that would give those iterations \(ie not referring to the unknown maximum `timefac` value in the simulation\).
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2629/mol_pseudoevoluti…
#2637: update Cactus COPYRIGHT file to LGPL
Reporter: Roland Haas
Status: new
Milestone:
Version:
Type: bug
Priority: minor
Component: Cactus
Cactus originally used GPL with a spelled out exception that linking is exempt from the GPL copyleft, ie the same exception that LGPL used. Some time ago \(around 2015, see eg commit logs in cactus\* repos\) we switched to explicitly using LGPL and updated the README files. Cactus itself however was forgotten \(has no README file, but has COPYRIGHT\). Pull request [https://bitbucket.org/cactuscode/cactus/pull-requests/144](https://bitbucke… addresses this.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2637/update-cactus-cop…