Code Monkey home page Code Monkey logo

delaunaysparse's People

Contributors

tchlux avatar thchang avatar

Stargazers

avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar avatar

Watchers

avatar avatar avatar

Forkers

gillette7 awvwgk daniellengyel ljuillen thchang

delaunaysparse's Issues

Extract original input points from DelaunaySparse output

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.

How to set pmode in example.py

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?

IERR size should be M x IR?

DelaunaySparse/python/delsparse.py

Line 86 in 6a80558

! IERR(1:M) contains integer valued error flags associated with the

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.

Weird outputs in interp_out

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.

Rank mismatch warning during make

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.

python bindings not creating delsparse_clib.so

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.

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.