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

different outputs from Circuitscape Python and Julia in advanced mode #414

Closed
cbogdanc opened this issue Apr 12, 2024 · 2 comments · Fixed by #421
Closed

different outputs from Circuitscape Python and Julia in advanced mode #414

cbogdanc opened this issue Apr 12, 2024 · 2 comments · Fixed by #421

Comments

@cbogdanc
Copy link

Hello,

I have encountered an issue where I'm getting two different outputs when using the same input and configuration files in Circuitscape Python vs. Julia. The model setup has a rectangular shape resistance map, with source nodes around the perimeter (with discrete x,y coordinates in a text file) and ground nodes in a narrow band represented as raster cells (with the same node ID). In the advanced mode, the cumulative current output from the Python version has values in the 0-1 range, and in the Julia version has values in the 0-46 million range, and a different spatial distribution.

Is this a known issue? Is there a way to obtain the same (or very similar) output from the two Circuitscape versions?

Thank you,
Bogdan

@cbogdanc
Copy link
Author

cbogdanc commented May 10, 2024

As a comparison, the outputs in pairwise mode are exactly the same in the Python and Julia version. The test was done using the same resistance map and the same source nodes file (discrete points) as for the advanced mode, and a set of ground nodes represented as discrete points (along the narrow band described above).

Here is a test dataset with input and output files for the advanced mode simulations.
I appreciate any insight that I can get into the model behavior differences between the two Circuitscape versions.

Thank you,
Bogdan

Circuitscape_TestData.zip

@ranjanan
Copy link
Member

I didn't support use_unit_currents just yet. I've added support now through #421. I've verified that it produces the same output as the python package.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants