vtopt / delaunaysparse Goto Github PK
View Code? Open in Web Editor NEWInterpolation via a Sparse Subset of the Delaunay Triangulation
Home Page: https://vtopt.github.io/DelaunaySparse
License: Other
Interpolation via a Sparse Subset of the Delaunay Triangulation
Home Page: https://vtopt.github.io/DelaunaySparse
License: Other
I'm now encountering this error when trying to run example.py
from the python subfolder. The first part are slatec
errors, which are explained in Issue # 8. The issue here is that (possibly because of that?) it's not creating the file delsparsee_clib.so
, which is required for the python bindings to work. Here's the error; system info at the bottom of the issue:
$ python example.py
Running command
gfortran -fPIC -shared -O3 -fopenmp -std=legacy delsparse_src/blas.f delsparse_src/lapack.f delsparse_src/real_precision.f90 delsparse_src/slatec.f delsparse_src/delsparse.f90 delsparse_src/delsparse_bind_c.f90 -o /Users/gillette7/Desktop/dataqual - results/delSp-repo/DelaunaySparse/python/delsparse_clib.so
delsparse_src/slatec.f:953:24:
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
953 | CALL DCOPY (N, WS, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
delsparse_src/slatec.f:965:25:
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
965 | CALL DCOPY (KRANK, W, MDW+1, WS(N2), 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-2)
delsparse_src/slatec.f:3040:25:
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
3040 | CALL DCOPY(N,X,0,X,1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
delsparse_src/slatec.f:3955:28:
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
3955 | CALL DCOPY (NSOLN, Z, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
delsparse_src/slatec.f:4123:25:
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
4123 | CALL DCOPY (NSOLN, Z, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
delsparse_src/slatec.f:595:26:
595 | 160 CALL DCOPY (KRANKE, W, MDW+1, WS(KRANKE+1), 1)
| 1
......
860 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-2)
gfortran: error: -E or -x required when input is from standard input
Traceback (most recent call last):
File "example.py", line 3, in <module>
import delsparse
File "/Users/gillette7/Desktop/dataqual - results/delSp-repo/DelaunaySparse/python/delsparse.py", line 51, in <module>
delsparse_clib = ctypes.CDLL(path_to_lib)
File "/opt/miniconda3/envs/gbcNN/lib/python3.7/ctypes/__init__.py", line 364, in __init__
self._handle = _dlopen(self._name, mode)
OSError: dlopen(/Users/gillette7/Desktop/dataqual - results/delSp-repo/DelaunaySparse/python/delsparse_clib.so, 6): image not found
MacOS BigSur 11.7.6
intel processor
$ gfortran --version
GNU Fortran (Homebrew GCC 13.1.0) 13.1.0
$ python --version
Python 3.7.11
Running in a conda env
Likely this is a gfortran issue as I am able to run on a linux machine that has
$ gfortran --version
GNU Fortran (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0
Any fix welcome, even if it's just advice on how to compel the use of an older version of gfortran
.
When I pass pts
as an input to DelaunaySparse (with the python bindings) it overwrites pts
with the shifted and scaled version after the run. I need to look up information from the pre-scaled pts
data set after running DelaunaySparse. Is there an easy way to do this, or do I have to make a deep copy of pts
before passing it in? Thanks.
When I run example.py
, it gives me error code 90: PMODE not set. (Although it does compute something). It looks like the pmode
variable is never set in the example.py
script, but no matter what I try (setting pmode
at different parts of the script, changing the value to 1 or 2, etc), I keep getting the error code 90. Any thoughts on what I'm missing here?
I'm running some examples where the input data points form a regular lattice in 2D and the points for interpolation sometimes lie on the edges or diagonals of this grid. I've found that the output values from interp_out
in this case are sometimes things like 1e271
or nan
. (I'm using the python bindings). Is this a known issue, perhaps due to the non-uniqueness of the Delaunay triangulation on a lattice? Or is this (more likely?) some other bug in my code? Thanks.
DelaunaySparse/python/delsparse.py
Line 86 in 6a80558
When trying to use the python wrapper on a case with output dimension (IR) larger than 1, I got error code 21: The size of the error array IERR does not match the number of interpolation points M. When I changed the size of IERR to M x IR, it worked. I think the readme has to be changed (in two spots) to reflect this.
I've started getting this error when building DelaunaySparse:
src$ make
gfortran -c -fopenmp delsparse.f90 -o delsparse.o
gfortran -c -fopenmp -std=legacy slatec.f -o slatec.o
slatec.f:955:24:
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
955 | CALL DCOPY (N, WS, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
slatec.f:967:25:
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
967 | CALL DCOPY (KRANK, W, MDW+1, WS(N2), 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-2)
slatec.f:3046:25:
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
3046 | CALL DCOPY(N,X,0,X,1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
slatec.f:3963:28:
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
3963 | CALL DCOPY (NSOLN, Z, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
slatec.f:4131:25:
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
......
4131 | CALL DCOPY (NSOLN, Z, 1, X, 1)
| 1
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-1)
slatec.f:596:26:
596 | 160 CALL DCOPY (KRANKE, W, MDW+1, WS(KRANKE+1), 1)
| 1
......
862 | CALL DCOPY (N, 0.D0, 0, WS, 1)
| 2
Warning: Rank mismatch between actual argument at (1) and actual argument at (2) (scalar and rank-2)
gfortran -c -fopenmp lapack.f -o lapack.o
gfortran -c -fopenmp blas.f -o blas.o
gfortran -fopenmp samples.f90 delsparse.o slatec.o lapack.o blas.o -o samples
gfortran -fopenmp samplep.f90 delsparse.o slatec.o lapack.o blas.o -o samplep
gfortran -fopenmp test_install.f90 delsparse.o slatec.o lapack.o blas.o -o test_install
./test_install
The installation of DELAUNAYSPARSE appears correct.
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.