# LinearProgramming accuracy problem

I am working with the LinearProgramming to solve my problem. I am using Simplex and RevisedSimplex methods. They are giving me drastically different results for the same parameters. Examples are attached as plots. Generally solutions are also different from expected. It is rather not a problem in my code, have checked it many times. It is all probably because of accuracy issues - I have constraints like:

... = ~ 10^(-62)

This is something I cannot avoid since I work on a strictly derived scientific model. Is there any way to overcome my problems? If not Mathematica, maybe there are some other packages for LP known for much higher precision? In one paper I have found analogous problem solved with Mathematica's LP which must have operated with the same range of precision as in my case (10^(-62)).

Thanks for suggestions.

Replay

Category: mathematical optimization Time: 2016-07-28 Views: 0

## Related post

• ArcGIS Buffer accuracy problem? 2012-06-07

Once again I encounter something very strange with ArcGIS. I perform a buffer analysis on one side of a polyline. There seems to be a serious accuracy problem and I can't figure out what the cause is. So I hope to find some expertise out there...What

• GPS accuracy problems outside vs inside with repeater 2013-12-19

I have an interesting problem testing GPS receivers: While working inside our lab, we get reasonable results. We have a GPS repeater, as shown in the picture. Inside the lab, I get good accuracy (between 3m to 5m error). On the roof, I get more than

• OpenLayers.Control.DrawFeature Accuracy problem 2011-12-01

When i draw a polygon using the OpenLayers.Control.DrawFeature then i export that polygon to KML using OpenLayers.Format.KML. Though, when I open it in google earth there is a big difference between the areas the polygon crosses even though the verti

• Open street map accuracy problem 2011-12-08

I want to locate shops with open street map. For that I send a query to nominatim, a webservice which geocodes coordinates from addresses. I draw a marker at the returned coordinates, and it´s inaccurate about 100 m. (~300 ft). To double check, I sea

• Circle vs. Rect Collision accuracy problem 2013-11-25

The blue ball is dynamic and is supposed to bounce off the walls, as well as the rectangle. Currently, the green rectangle is adding an invisible line to the scene, which doesn't let my ball go where I want it to. You can see the issue on the picture

• Android 4.1 touch accuracy 2013-03-02

I have upgraded my Galaxy S2 to Android 4.1.2 XXLSJ.Now I'm having accuracy problem in games as basketball touch and overkill.If I try to aim slowly, it starts to shake and phone can't understand the exact point that I want to touch. I remember that

• Point rotation around point, accuracy question 2014-04-19

I study OpenGL ES 2.0. But I think it's more a C++ question rather then an OpenGL one. I am stuck with the following rotation question: It is known, that rotation transformation can be applied using these equations (for simplicity suppose we rotating

• What is my problem with MeshFunction? 2014-08-25

I want plot the max value in a sine plot. We can use the following code Plot[Sin[x], {x, 0, 50}, Mesh -> {{.99}}, MeshFunctions -> {#2 &}, MeshStyle -> {PointSize[Large], Red}] Why can't Mesh -> 1 be used to plot the red point? But in a ma

• Scaling PID (Proportional Integral Derivative) Output 2010-10-12

I have implemented a PID function using the formula, correction = Kp * error + Kd * (error - prevError) + kI * (sum of errors) What should I do to keep my output between a certain range? say 0-255 If I disregard any value not between 0 to 255 it prod

• How to export 3D models that consist of several parts (eg. turret on a tank)? 2011-01-26

What are the standard alternatives for the mechanics of attaching turrets and such to 3D models for use in-game? I don't mean the logic, but rather the graphics aspects. My naive approach is to extend the MD2-like format that I'm using (blender-expor

• how to calibrate temperature sensor IC DS1621? 2011-03-28

I have made temperature indicator using microcontroller P89V51. My temperature sensing IC is DS1621. I'm using LCD for the display..n ckt is working very fine, Do I have to calibrate this IC? If yes, then how can I? --------------Solutions-----------

• Help on STC3100 calculations 2011-07-27

I would rephrase my question: What I'm trying to do is to measure the battery voltage. Based on STC3100 datasheet, the battery voltage value is located on address 0x08 and 0x09 where the value 0x00E6 was read from. The actual value of the voltage bas

• XST Verilog - Casting real to integer constants 2011-12-11

When I try to synthesize the following Verilog code using Xilinx XST, I get the error, "Unsupported real constant". If I try wrapping that expression in an \$rtoi function, XST gives a different error: "Unsupported System Function Call"

• Infinite expression 1/0. encountered - caused by precision? 2012-05-17

Planck's law dependent on frequency rho is as follows: B[T_, h_, rho_, k_, c_] := (2 h rho^3)/c^2 1/(Exp[h rho/(T k)] - 1) As you can see, the denominator can only be zero for rho = 0. All the rest are nonzero constants. Curiously if we try to plot i

• Can the solar system be accurately represented in 3d space using doubles (or longs)? 2012-05-22

I would like know to how to best manage coordinates in a 3D game whose aim is to realistically model the entire solar system, yet be able to handle the smallest movements in a "ship" (ie: perhaps we can consider 1cm to be the smallest acceptable

• pgfplots plots function not correct 2012-10-10

This code \documentclass{article} \usepackage{pgfplots} \pgfplotsset{compat=newest} \begin{document} \begin{tikzpicture} \begin{axis} \addplot [domain=15:30, samples=1000] {4*65*sin(0.25*asin(x/65))/x}; \end{axis} \end{tikzpicture} \end{document} plo

• Backsubstituting solution into FindRoot and Plot 2014-07-11

I want to find roots from a transcendental equation using FindRoot as a function of a real parameter and then plot the solutions vs. the parameter. The following commands work perfectly: equ[a_?NumericQ] := FindRoot[x == E^(-a x), {x, 1}][[1, 2]]; Pl

• Interpolating an Antiderivative 2014-09-17

I'd like to be able to make InterpolatingFunctions for antiderivatives of functions that can't be integrated symbolically. However, the following code returns several error messages: FunctionInterpolation[Integrate[Sqrt[1 + x^3], {x, 0, t}], {t, 0, 1

• How do I perform the following numerical integration 2014-09-17

I'd like to be able to make InterpolatingFunctions for antiderivatives of functions that can't be integrated symbolically. However, the following code returns several error messages: FunctionInterpolation[Integrate[Sqrt[1 + x^3], {x, 0, t}], {t, 0, 1