Re: COMPUTE_POTENTIAL_ENERGY

From: Andrew Benson <cond0025_at_herald.ox.ac.uk>
Date: Mon, 16 Jan 2006 09:12:44 +0000

Hi Volker,

> I have made these changes in the version of the code that can be
> downloaded from gadget's web-site. So I'd suggest to refetch the code
> and try it again. (The only files that have changed are potential.c,
> pm_nonperiodic.c, and proto.h) If I'm right about my suspicion, it
> should work now - if not, please let me know.

Thanks for tracking this down.

potential.c failed to compile initially. This was due to a missing ";" at line
289 - once added the code compiled just fine.

I tried re-running a test case in which the potential is re-computed
frequently. The code still crashes unfortunately. As far as I can tell, the
crash occurs during the call to:

force_treeevaluate_potential_shortrange(i, 0);

which is made in potential.c. Could this be a similar problem - i.e. the
allowed region fo space needs to be recomputed before
force_treeevaluate_potential_shortrange is called?

Cheers,
Andrew.

-- 
--------------------------------------------------------------
|                  |                                         |
| Andrew J. Benson | e-mail:   abenson_at_astro.ox.ac.uk        |
| Astrophysics     | Phone:    +44 (0)1865 273296            |
| Keble Road       | Cell:     +44 (0)7793 058893            |
| Oxford           | Fax:      +44 (0)8712 476655 (personal) |
| OX1 3RH          | Fax:       +1 501 638 7383 (personal US)|
| U.K.             | Fax:      +44 (0)1865 273390     (work) |
|                  | Homepage:                               |
|                  |     www-astro.physics.ox.ac.uk/~abenson |
|                  |                                         |
--------------------------------------------------------------
Received on 2006-01-16 10:38:57

This archive was generated by hypermail 2.3.0 : 2022-09-01 14:03:41 CEST