scipy curve_fit doesn’t like math module

by

Be careful with numpy-arrays, operations working on arrays and operations working on scalars!

Scipy optimize assumes the input (initial-point) to be a 1d-array and often things go wrong in other cases (a list for example becomes an array and if you assumed to work on lists, things go havoc; those kind of problems are common here on StackOverflow and debugging is not that easy to do by the eye; code-interaction helps!).

import numpy as np
import math

x = np.ones(1)

np.sin(x)
> array([0.84147098])

math.sin(x)
> 0.8414709848078965                     # this only works as numpy has dedicated support
                                         # as indicated by the error-msg below!
x = np.ones(2)

np.sin(x)
> array([0.84147098, 0.84147098])

math.sin(x)
> TypeError: only size-1 arrays can be converted to Python scalars

To be honest: this is part of a very basic understanding of numpy and should be understood when using scipy’s somewhat sensitive functions.

More Related Contents:

  1. How to do exponential and logarithmic curve fitting in Python? I found only polynomial fitting
  2. How to apply piecewise linear fit in Python?
  3. python numpy/scipy curve fitting
  4. Import Error: No module named numpy
  5. Installing scipy in Python 3.5 on 32-bit Windows 7 Machine
  6. Error “Import Error: No module named numpy” on Windows
  7. Fitting a closed curve to a set of points
  8. How to share numpy random state of a parent process with child processes?
  9. Interpolating a closed curve using scipy
  10. Computing the correlation coefficient between two multi-dimensional arrays
  11. Calculate average of every x rows in a table and create new table
  12. How does condensed distance matrix work? (pdist)
  13. ImportError in importing from sklearn: cannot import name check_build
  14. How do I fit a sine curve to my data with pylab and numpy?
  15. binning data in python with scipy/numpy
  16. Transform a set of numbers in numpy so that each number gets converted into a number of other numbers which are less than it
  17. Calculating Pearson correlation and significance in Python
  18. Calculating the area under a curve given a set of coordinates, without knowing the function
  19. TypeError: only integer scalar arrays can be converted to a scalar index with 1D numpy indices array
  20. scipy.curve_fit() returns multiple lines
  21. How to calculate cumulative normal distribution?
  22. Concatenate sparse matrices in Python using SciPy/Numpy
  23. How to replicate array to specific length array
  24. How to display progress of scipy.optimize function?
  25. Fitting a Weibull distribution using Scipy
  26. Is the order of a Python dictionary guaranteed over iterations?
  27. Is there a numpy/scipy dot product, calculating only the diagonal entries of the result?
  28. sampling random floats on a range in numpy
  29. Fast interpolation of grid data
  30. In sklearn.decomposition.PCA, why are components_ negative?

Leave a Comment