foscat 3.2.0__tar.gz → 3.3.0__tar.gz

{foscat-3.2.0/src/foscat.egg-info → foscat-3.3.0}/PKG-INFO

@@ -1,11 +1,8 @@
 Metadata-Version: 2.1
 Name: foscat
-Version: 3.2.0
+Version: 3.3.0
 Summary: Generate synthetic Healpix or 2D data using Cross Scattering Transform
-Home-page: https://github.com/jmdelouis/FOSCAT
-Author: Jean-Marc DELOUIS
 Author-email: Jean-Marc DELOUIS <jean.marc.delouis@ifremer.fr>
-Maintainer: Theo Foulquier
 Maintainer-email: Theo Foulquier <theo.foulquier@ifremer.fr>
 License: BSD-3-Clause
 Project-URL: Repository, https://github.com/jmdelouis/FOSCAT.git
@@ -22,12 +19,6 @@ Classifier: Programming Language :: Python :: 3.12
 Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENCE
-Requires-Dist: imageio
-Requires-Dist: imagecodecs
-Requires-Dist: matplotlib
-Requires-Dist: numpy
-Requires-Dist: tensorflow
-Requires-Dist: healpy
 
 # foscat
 
@@ -49,23 +40,27 @@ https://github.com/IAOCEA/demo-foscat-pangeo-eosc/blob/main/Demo_Synthesis.ipynb
 
 # FOSCAT_DEMO
 
-The python scripts *demo.py* included in this package demonstrate how to use the foscat library to generate synthetic fields that have patterns with the same statistical properties as a specified image.
+The python scripts _demo.py_ included in this package demonstrate how to use the foscat library to generate synthetic fields that have patterns with the same statistical properties as a specified image.
 
 # Install foscat library
 
-Before installing, make sure you have python installed in your enviroment.  
+Before installing, make sure you have python installed in your enviroment.
 The last version of the foscat library can be installed using PyPi:
+
 ```
 pip install foscat
 ```
+
 Load the FOSCAT_DEMO package from github.
+
 ```
 git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 ```
 
+## Recommended installing procedures for mac users
+
+It is recomended to use python=3.9\*.
 
-##  Recommended installing procedures  for mac users 
-It is recomended to use python=3.9*.  
 ```
 micromamba create -n FOSCAT
 micromamba install -n FOSCAT ‘python==3.9*’
@@ -75,8 +70,9 @@ git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 
 ```
 
-##  Recommended installing procedures  HPC users 
-It is recomended to install tensorflow in advance.  For [DATARMOR](https://pcdm.ifremer.fr/Equipement)  for using GPU ; 
+## Recommended installing procedures HPC users
+
+It is recomended to install tensorflow in advance. For [DATARMOR](https://pcdm.ifremer.fr/Equipement) for using GPU ;
 
 ```
 micromamba create -n FOSCAT
@@ -91,33 +87,38 @@ git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 # Spherical data example
 
 ## compute a synthetic image
+
 ```
 python demo.py -n=32 -k -c -s=100
 ```
-The *demo.py* script serves as a demonstration of the capabilities of the foscat library. It utilizes the Cross Wavelet Scattering Transform to generate a Healpix map that possesses the same characteristics as a specified input map. 
-- ```-n=32``` computes map with nside=32.
-- ```-k``` uses 5x5 kernel.
-- ```-c``` uses Scattering Covariance.
-- ```-l``` uses LBFGS minimizer.
-- ```-s=100``` computes 100 steps. 
+
+The _demo.py_ script serves as a demonstration of the capabilities of the foscat library. It utilizes the Cross Wavelet Scattering Transform to generate a Healpix map that possesses the same characteristics as a specified input map.
+
+- `-n=32` computes map with nside=32.
+- `-k` uses 5x5 kernel.
+- `-c` uses Scattering Covariance.
+- `-l` uses LBFGS minimizer.
+- `-s=100` computes 100 steps.
+
 ```
 python demo.py -n=8 [-c|--cov][-s|--steps=3000][-S=1234|--seed=1234][-k|--k5x5][-d|--data][-o|--out][-r|--orient] [-p|--path][-a|--adam]
 
 ```
-* The "-n" option specifies the nside of the input map. The maximum nside value is 256 with the default map.
-* The "--cov" option (optional) uses scat_cov instead of scat.
-* The "--steps" option (optional) specifies the number of iterations. If not specified, the default value is 1000.
-* The "--seed" option (optional) specifies the seed of the random generator.
-* The "--path" option (optional) allows you to define the path where the output files will be written. The default path is "data".
-* The "--k5x5" option (optional) uses a 5x5 kernel instead of a 3x3.
-* The "--data" option (optional) specifies the input data file to be used. If not specified, the default file "LSS_map_nside128.npy" will be used.
-* The "--out" option (optional) specifies the output file name. If not specified, the output file will be saved in "demo".
-* The "--orient" option (optional) specifies the number of orientations. If not specified, the default value is 4.
-* The "--adam" option (optional) makes the synthesis using the ADAM optimizer instead of the L_BFGS.
+
+- The "-n" option specifies the nside of the input map. The maximum nside value is 256 with the default map.
+- The "--cov" option (optional) uses scat_cov instead of scat.
+- The "--steps" option (optional) specifies the number of iterations. If not specified, the default value is 1000.
+- The "--seed" option (optional) specifies the seed of the random generator.
+- The "--path" option (optional) allows you to define the path where the output files will be written. The default path is "data".
+- The "--k5x5" option (optional) uses a 5x5 kernel instead of a 3x3.
+- The "--data" option (optional) specifies the input data file to be used. If not specified, the default file "LSS_map_nside128.npy" will be used.
+- The "--out" option (optional) specifies the output file name. If not specified, the output file will be saved in "demo".
+- The "--orient" option (optional) specifies the number of orientations. If not specified, the default value is 4.
+- The "--adam" option (optional) makes the synthesis using the ADAM optimizer instead of the L_BFGS.
 
 ## plot the result
 
-The following script generates a series of plots that showcase different aspects of the synthesis process using the *demo.py* script.
+The following script generates a series of plots that showcase different aspects of the synthesis process using the _demo.py_ script.
 
 > python test2D.py
 
@@ -131,7 +132,7 @@ python plotdemo.py -n=32 -c
 
 # compute a synthetic turbulent field
 
-The python scripts *demo2D.py* included in this package demonstrate how to use the foscat library to generate a 2D synthetic fields that have patterns with the same statistical properties as a specified 2D image. In this particular case, the input field is a sea surface temperature extracted from a north atlantic ocean simulation.
+The python scripts _demo2D.py_ included in this package demonstrate how to use the foscat library to generate a 2D synthetic fields that have patterns with the same statistical properties as a specified 2D image. In this particular case, the input field is a sea surface temperature extracted from a north atlantic ocean simulation.
 
 > python testHealpix.py
 
@@ -141,17 +142,19 @@ python demo2d.py -n=32 -k -c
 
 > python testHplot.py
 
-The following script generates a series of plots that showcase different aspects of the synthesis process using the *demo2D.py* script.
+The following script generates a series of plots that showcase different aspects of the synthesis process using the _demo2D.py_ script.
+
 ```
 python plotdemo2d.py -n=32 -c
 ```
+
 For more information, see the [documentation](https://foscat-documentation.readthedocs.io/en/latest/index.html).
 
 > mpirun -np 3 testHealpix_mpi.py
 
 ## Authors and acknowledgment
 
-Authors: J.-M. Delouis, T. Foulquier, J. Mangin, L. Mousset, T. Odaka, F. Paul, E. Allys
+Authors: J.-M. Delouis, P. Campeti, T. Foulquier, J. Mangin, L. Mousset, T. Odaka, F. Paul, E. Allys
 
 This work is part of the R & T Deepsee project supported by CNES. The authors acknowledge the heritage of the Planck-HFI consortium regarding data, software, knowledge. This work has been supported by the Programme National de Télédétection Spatiale (PNTS, http://programmes.insu.cnrs.fr/pnts/), grant n◦ PNTS-2020-08
 

{foscat-3.2.0 → foscat-3.3.0}/README.md

@@ -18,23 +18,27 @@ https://github.com/IAOCEA/demo-foscat-pangeo-eosc/blob/main/Demo_Synthesis.ipynb
 
 # FOSCAT_DEMO
 
-The python scripts *demo.py* included in this package demonstrate how to use the foscat library to generate synthetic fields that have patterns with the same statistical properties as a specified image.
+The python scripts _demo.py_ included in this package demonstrate how to use the foscat library to generate synthetic fields that have patterns with the same statistical properties as a specified image.
 
 # Install foscat library
 
-Before installing, make sure you have python installed in your enviroment.  
+Before installing, make sure you have python installed in your enviroment.
 The last version of the foscat library can be installed using PyPi:
+
 ```
 pip install foscat
 ```
+
 Load the FOSCAT_DEMO package from github.
+
 ```
 git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 ```
 
+## Recommended installing procedures for mac users
+
+It is recomended to use python=3.9\*.
 
-##  Recommended installing procedures  for mac users 
-It is recomended to use python=3.9*.  
 ```
 micromamba create -n FOSCAT
 micromamba install -n FOSCAT ‘python==3.9*’
@@ -44,8 +48,9 @@ git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 
 ```
 
-##  Recommended installing procedures  HPC users 
-It is recomended to install tensorflow in advance.  For [DATARMOR](https://pcdm.ifremer.fr/Equipement)  for using GPU ; 
+## Recommended installing procedures HPC users
+
+It is recomended to install tensorflow in advance. For [DATARMOR](https://pcdm.ifremer.fr/Equipement) for using GPU ;
 
 ```
 micromamba create -n FOSCAT
@@ -60,33 +65,38 @@ git clone https://github.com/jmdelouis/FOSCAT_DEMO.git
 # Spherical data example
 
 ## compute a synthetic image
+
 ```
 python demo.py -n=32 -k -c -s=100
 ```
-The *demo.py* script serves as a demonstration of the capabilities of the foscat library. It utilizes the Cross Wavelet Scattering Transform to generate a Healpix map that possesses the same characteristics as a specified input map. 
-- ```-n=32``` computes map with nside=32.
-- ```-k``` uses 5x5 kernel.
-- ```-c``` uses Scattering Covariance.
-- ```-l``` uses LBFGS minimizer.
-- ```-s=100``` computes 100 steps. 
+
+The _demo.py_ script serves as a demonstration of the capabilities of the foscat library. It utilizes the Cross Wavelet Scattering Transform to generate a Healpix map that possesses the same characteristics as a specified input map.
+
+- `-n=32` computes map with nside=32.
+- `-k` uses 5x5 kernel.
+- `-c` uses Scattering Covariance.
+- `-l` uses LBFGS minimizer.
+- `-s=100` computes 100 steps.
+
 ```
 python demo.py -n=8 [-c|--cov][-s|--steps=3000][-S=1234|--seed=1234][-k|--k5x5][-d|--data][-o|--out][-r|--orient] [-p|--path][-a|--adam]
 
 ```
-* The "-n" option specifies the nside of the input map. The maximum nside value is 256 with the default map.
-* The "--cov" option (optional) uses scat_cov instead of scat.
-* The "--steps" option (optional) specifies the number of iterations. If not specified, the default value is 1000.
-* The "--seed" option (optional) specifies the seed of the random generator.
-* The "--path" option (optional) allows you to define the path where the output files will be written. The default path is "data".
-* The "--k5x5" option (optional) uses a 5x5 kernel instead of a 3x3.
-* The "--data" option (optional) specifies the input data file to be used. If not specified, the default file "LSS_map_nside128.npy" will be used.
-* The "--out" option (optional) specifies the output file name. If not specified, the output file will be saved in "demo".
-* The "--orient" option (optional) specifies the number of orientations. If not specified, the default value is 4.
-* The "--adam" option (optional) makes the synthesis using the ADAM optimizer instead of the L_BFGS.
+
+- The "-n" option specifies the nside of the input map. The maximum nside value is 256 with the default map.
+- The "--cov" option (optional) uses scat_cov instead of scat.
+- The "--steps" option (optional) specifies the number of iterations. If not specified, the default value is 1000.
+- The "--seed" option (optional) specifies the seed of the random generator.
+- The "--path" option (optional) allows you to define the path where the output files will be written. The default path is "data".
+- The "--k5x5" option (optional) uses a 5x5 kernel instead of a 3x3.
+- The "--data" option (optional) specifies the input data file to be used. If not specified, the default file "LSS_map_nside128.npy" will be used.
+- The "--out" option (optional) specifies the output file name. If not specified, the output file will be saved in "demo".
+- The "--orient" option (optional) specifies the number of orientations. If not specified, the default value is 4.
+- The "--adam" option (optional) makes the synthesis using the ADAM optimizer instead of the L_BFGS.
 
 ## plot the result
 
-The following script generates a series of plots that showcase different aspects of the synthesis process using the *demo.py* script.
+The following script generates a series of plots that showcase different aspects of the synthesis process using the _demo.py_ script.
 
 > python test2D.py
 
@@ -100,7 +110,7 @@ python plotdemo.py -n=32 -c
 
 # compute a synthetic turbulent field
 
-The python scripts *demo2D.py* included in this package demonstrate how to use the foscat library to generate a 2D synthetic fields that have patterns with the same statistical properties as a specified 2D image. In this particular case, the input field is a sea surface temperature extracted from a north atlantic ocean simulation.
+The python scripts _demo2D.py_ included in this package demonstrate how to use the foscat library to generate a 2D synthetic fields that have patterns with the same statistical properties as a specified 2D image. In this particular case, the input field is a sea surface temperature extracted from a north atlantic ocean simulation.
 
 > python testHealpix.py
 
@@ -110,17 +120,19 @@ python demo2d.py -n=32 -k -c
 
 > python testHplot.py
 
-The following script generates a series of plots that showcase different aspects of the synthesis process using the *demo2D.py* script.
+The following script generates a series of plots that showcase different aspects of the synthesis process using the _demo2D.py_ script.
+
 ```
 python plotdemo2d.py -n=32 -c
 ```
+
 For more information, see the [documentation](https://foscat-documentation.readthedocs.io/en/latest/index.html).
 
 > mpirun -np 3 testHealpix_mpi.py
 
 ## Authors and acknowledgment
 
-Authors: J.-M. Delouis, T. Foulquier, J. Mangin, L. Mousset, T. Odaka, F. Paul, E. Allys
+Authors: J.-M. Delouis, P. Campeti, T. Foulquier, J. Mangin, L. Mousset, T. Odaka, F. Paul, E. Allys
 
 This work is part of the R & T Deepsee project supported by CNES. The authors acknowledge the heritage of the Planck-HFI consortium regarding data, software, knowledge. This work has been supported by the Programme National de Télédétection Spatiale (PNTS, http://programmes.insu.cnrs.fr/pnts/), grant n◦ PNTS-2020-08
 

{foscat-3.2.0 → foscat-3.3.0}/pyproject.toml

@@ -1,14 +1,14 @@
 [project]
 name = "foscat"
-version = "3.2.0"
+version = "3.3.0"
 description = "Generate synthetic Healpix or 2D data using Cross Scattering Transform"
 readme = "README.md"
-license = {text = "BSD-3-Clause"}
+license = { text = "BSD-3-Clause" }
 authors = [
-  {name = "Jean-Marc DELOUIS", email = "jean.marc.delouis@ifremer.fr"},
+  { name = "Jean-Marc DELOUIS", email = "jean.marc.delouis@ifremer.fr" },
 ]
 maintainers = [
-  {name = "Theo Foulquier", email = "theo.foulquier@ifremer.fr"},
+  { name = "Theo Foulquier", email = "theo.foulquier@ifremer.fr" },
 ]
 dependencies = [
   "imageio",
@@ -36,7 +36,7 @@ Issues = "https://github.com/jmdelouis/FOSCAT/issues"
 Documentation = "https://foscat-documentation.readthedocs.io/en/latest/index.html"
 
 [tool.setuptools]
-package-dir = {"" = "src"}
+package-dir = { "" = "src" }
 packages = ["foscat"]
 
 [build-system]
@@ -47,38 +47,33 @@ build-backend = "setuptools.build_meta"
 target-version = "py39"
 builtins = ["ellipsis"]
 exclude = [
-    ".git",
-    ".eggs",
-    "build",
-    "dist",
-    "__pycache__",
+  ".git",
+  ".eggs",
+  "build",
+  "dist",
+  "__pycache__",
 ]
 line-length = 100
 
 [tool.ruff.lint]
 ignore = [
-    "E402", # E402: module level import not at top of file
-    "E501", # E501: line too long - let black worry about that
-    "E731", # E731: do not assign a lambda expression, use a def
-    # ignore for now, fix in a later PR
-    #"E741", # ambiguous variable name
-    "E722", # bare except
-    "UP031", # percent-based string interpolation
-    #"F401", # unused imports
-    #"F811", # redefinition of unused name
-    #"F821", # undefined name
-    #"F841", # unused local variable
+  "E402", # E402: module level import not at top of file
+  "E501", # E501: line too long - let black worry about that
+  "E731", # E731: do not assign a lambda expression, use a def
+  # ignore for now, fix in a later PR
+  "E722",  # bare except
+  "UP031", # percent-based string interpolation
 ]
 select = [
-    "F",  # Pyflakes
-    "E",  # Pycodestyle
-    "I",  # isort
-    "UP", # Pyupgrade
-    "TID", # flake8-tidy-imports
-    "W",
+  "F",   # Pyflakes
+  "E",   # Pycodestyle
+  "I",   # isort
+  "UP",  # Pyupgrade
+  "TID", # flake8-tidy-imports
+  "W",
 ]
 extend-safe-fixes = [
-    "TID252", # absolute imports
+  "TID252", # absolute imports
 ]
 fixable = ["I", "TID252"]
 

foscat-3.3.0/src/foscat/CircSpline.py

@@ -0,0 +1,75 @@
+import numpy as np
+
+
+class CircSpline:
+    def __init__(self, nodes, degree=3):
+        """
+        Initializes the Spline1D instance.
+
+        Parameters:
+        - nodes (int): The number of nodes in the spline.
+        - degree (int): The degree of the spline. Default is 3.
+        """
+        self.degree = degree
+        self.nodes = nodes
+        
+        
+    def cubic_spline_function(self,x):
+        """
+        Evaluate the cubic spline basis function.
+    
+        Args:
+            x (float or array): Input value(s) to evaluate the spline basis function.
+    
+        Returns:
+            float or array: Result of the cubic spline basis function.
+        """
+        return -2 * x**3 + 3 * x**2
+
+
+    def eval(self,x):
+        """
+        Compute a 3rd-degree cubic spline with 4-point support.
+    
+        Args:
+            x (float or array): Input value(s) to compute the spline.
+
+        Returns:
+            indices (array): Indices of the spline support points.
+            coefficients (array): Normalized spline coefficients.
+        """
+        N=self.nodes
+        
+        if isinstance(x, float):
+            # Single scalar input
+            base_idx = int(x * (N))
+            indices = np.zeros([4], dtype="int")
+            coefficients = np.zeros([4])
+        else:
+            # Array input
+            base_idx = (x * (N)).astype("int")
+            indices = np.zeros([4, x.shape[0]], dtype="int")
+            coefficients = np.zeros([4, x.shape[0]])
+
+        # Compute the fractional part of the input
+        fractional_part = x * (N) - base_idx
+
+        # Compute spline coefficients for 4 support points
+        coefficients[3] = self.cubic_spline_function(fractional_part / 2) / 2
+        coefficients[2] = self.cubic_spline_function(0.5 + fractional_part / 2) / 2
+        coefficients[1] = self.cubic_spline_function(1 - fractional_part / 2) / 2
+        coefficients[0] = self.cubic_spline_function(0.5 - fractional_part / 2) / 2
+
+        # Assign indices for the support points
+        indices[3] = (base_idx + 3)%N
+        indices[2] = (base_idx + 2)%N
+        indices[1] = (base_idx + 1)%N
+        indices[0] = base_idx
+
+        # Adjust indices to start from 0
+        indices = indices - 1
+        # Square coefficients and normalize
+        coefficients = coefficients * coefficients
+        coefficients /= np.sum(coefficients, axis=0)
+
+        return indices, coefficients

{foscat-3.2.0 → foscat-3.3.0}/src/foscat/FoCUS.py

@@ -38,7 +38,7 @@ class FoCUS:
         mpi_rank=0,
     ):
 
-        self.__version__ = "3.2.0"
+        self.__version__ = "3.3.0"
         # P00 coeff for normalization for scat_cov
         self.TMPFILE_VERSION = TMPFILE_VERSION
         self.P1_dic = None
@@ -2193,7 +2193,7 @@ class FoCUS:
                         + ishape[axis + 2 :],
                     )
 
-            return self.backend.bk_reshape(res, in_image.shape+[self.NORIENT])
+            return self.backend.bk_reshape(res, in_image.shape + [self.NORIENT])
         elif self.use_1D:
             ishape = list(in_image.shape)
             if len(ishape) < axis + 1:
@@ -2265,7 +2265,7 @@ class FoCUS:
                         res, ishape[0:axis] + [res.shape[1]] + ishape[axis + 1 :]
                     )
 
-            return self.backend.bk_reshape(res, in_image.shape+[self.NORIENT])
+            return self.backend.bk_reshape(res, in_image.shape + [self.NORIENT])
 
         else:
             nside = int(np.sqrt(image.shape[axis] // 12))

{foscat-3.2.0 → foscat-3.3.0}/src/foscat/Softmax.py

@@ -1,4 +1,4 @@
-#import tensorflow as tf
+# import tensorflow as tf
 from tensorflow.keras.layers import Dense, Softmax
 from tensorflow.keras.models import Sequential
 

foscat-3.3.0/src/foscat/Spline1D.py

@@ -0,0 +1,92 @@
+import numpy as np
+
+class Spline1D:
+    def __init__(self, nodes, degree=3):
+        """
+        Initializes the Spline1D instance.
+
+        Parameters:
+        - nodes (int): The number of nodes in the spline.
+        - degree (int): The degree of the spline. Default is 3.
+        """
+        self.degree = degree
+        self.nodes = nodes
+        
+        
+    def cubic_spline_function(self,x):
+        """
+        Evaluate the cubic spline basis function.
+    
+        Args:
+            x (float or array): Input value(s) to evaluate the spline basis function.
+    
+        Returns:
+            float or array: Result of the cubic spline basis function.
+        """
+        return -2 * x**3 + 3 * x**2
+
+
+    def eval(self,x):
+        """
+        Compute a 3rd-degree cubic spline with 4-point support.
+    
+        Args:
+            x (float or array): Input value(s) to compute the spline.
+
+        Returns:
+            indices (array): Indices of the spline support points.
+            coefficients (array): Normalized spline coefficients.
+        """
+        N=self.nodes
+        
+        if isinstance(x, float):
+            # Single scalar input
+            base_idx = int(x * (N-1))
+            indices = np.zeros([4], dtype="int")
+            coefficients = np.zeros([4])
+        else:
+            # Array input
+            base_idx = (x * (N-1)).astype("int")
+            indices = np.zeros([4, x.shape[0]], dtype="int")
+            coefficients = np.zeros([4, x.shape[0]])
+
+        # Compute the fractional part of the input
+        fractional_part = x * (N-1) - base_idx
+
+        # Compute spline coefficients for 4 support points
+        coefficients[3] = self.cubic_spline_function(fractional_part / 2) / 2
+        coefficients[2] = self.cubic_spline_function(0.5 + fractional_part / 2) / 2
+        coefficients[1] = self.cubic_spline_function(1 - fractional_part / 2) / 2
+        coefficients[0] = self.cubic_spline_function(0.5 - fractional_part / 2) / 2
+
+        # Assign indices for the support points
+        indices[3] = base_idx + 3
+        indices[2] = base_idx + 2
+        indices[1] = base_idx + 1
+        indices[0] = base_idx
+
+        # Handle boundary conditions
+        if isinstance(x, float):
+            if indices[0] == 0:
+                indices[0] = 1
+            if indices[1] == 0:
+                indices[1] = 1
+            if indices[2] == N + 1:
+                indices[2] = N
+            if indices[3] == N + 1:
+                indices[3] = N
+            if indices[3] == N + 2:
+                indices[3] = N
+        else:
+            indices[0, indices[0] == 0] = 1
+            indices[1, indices[1] == 0] = 1
+            indices[2, indices[2] >= N + 1] = N
+            indices[3, indices[3] >= N + 1] = N
+
+        # Adjust indices to start from 0
+        indices = indices - 1
+        # Square coefficients and normalize
+        coefficients = coefficients * coefficients
+        coefficients /= np.sum(coefficients, axis=0)
+
+        return indices, coefficients