Notas para contenedor de docker:
Comando de docker para ejecución de la nota de forma local:
nota: cambiar <ruta a mi directorio> por la ruta de directorio que se desea mapear a /datos dentro del contenedor de docker.
docker run --rm -v <ruta a mi directorio>:/datos --name jupyterlab_numerical -p 8888:8888 -d palmoreck/jupyterlab_numerical:1.1.0password para jupyterlab: qwerty
Detener el contenedor de docker:
docker stop jupyterlab_numericalDocumentación de la imagen de docker palmoreck/jupyterlab_numerical:1.1.0 en liga.
In [1]:
!pip3 install --user -q cvxpy
In [1]:
import os
In [2]:
cur_directory = os.getcwd()
In [3]:
dir_alg_python = '/algoritmos/Python'
In [4]:
os.chdir(cur_directory + dir_alg_python)
In [5]:
import math
import numpy as np
from utils import compute_error
from algorithms_for_cieco import path_following_method_infeasible_init_point
In [6]:
fo = lambda x: x[0]**2 + x[1]**2 + x[2]**2 + x[3]**2-2*x[0]-3*x[3]
In [7]:
const = {0: lambda x: -x[0],
1: lambda x: -x[1],
2: lambda x: -x[2],
3: lambda x: -x[3]
}
In [8]:
A= np.array([[2,1,1,4],
[1,1,2,1]])
In [9]:
b=np.array([7,6])
In [10]:
x_ast=np.array([1.1232876712328763,0.6506849315068493,
1.8287671232876714,0.5684931506849317])
In [11]:
x_0 = np.array([1,1,1,1],dtype=float)
In [12]:
x_0
Out[12]:
In [13]:
nu_0 = np.array([0,0], dtype=float)
In [14]:
p_ast=fo(x_ast)
In [15]:
p_ast
Out[15]:
In [16]:
tol_outer_iter = 1e-6
tol=1e-8
tol_backtracking=1e-8
maxiter=30
mu=10
In [17]:
[x,iter_barrier,t] = path_following_method_infeasible_init_point(fo, A, b,
const,
x_0, nu_0, tol,
tol_backtracking, x_ast, p_ast, maxiter,
mu, tol_outer_iter = tol_outer_iter
)
In [18]:
[x,iter_barrier,t]
Out[18]:
In [19]:
compute_error(x_ast,x)
Out[19]:
In [20]:
x_ast
Out[20]:
In [21]:
import cvxpy as cp
In [22]:
x1 = cp.Variable()
x2 = cp.Variable()
x3 = cp.Variable()
x4 = cp.Variable()
In [23]:
# Create two constraints.
constraints = [2*x1+x2+x3+4*x4-7 == 0,x1+x2+2*x3+x4-6 == 0,x1>=0,x2>=0,x3>=0,x4>=0]
# Form objective.
obj = cp.Minimize(x1**2+x2**2+x3**2+x4**2-2*x1-3*x4)
In [24]:
# Form and solve problem.
prob = cp.Problem(obj, constraints)
prob.solve() # Returns the optimal value.
Out[24]:
In [25]:
print("status:", prob.status)
print("optimal value", prob.value)
print("optimal var", x1.value, x2.value, x3.value,x4.value)
In [26]:
fo = lambda x: 2*x[0] + 5*x[1]
In [27]:
const = {0: lambda x: 6-x[0]-x[1],
1: lambda x: -18+x[0]+2*x[1],
2: lambda x: -x[0],
3: lambda x: -x[1]
}
In [28]:
A=np.array([0,0],dtype=float)
b = 0
In [29]:
x_ast = np.array([6,0], dtype=float)
In [30]:
x_0 = np.array([4,4], dtype=float)
In [31]:
p_ast=fo(x_ast)
In [32]:
p_ast
Out[32]:
In [33]:
nu_0 = np.array([0,0], dtype=float)
In [34]:
tol_outer_iter = 1e-3
tol=1e-8
tol_backtracking=1e-8
maxiter=30
mu=10
In [35]:
[x,iter_barrier,t] = path_following_method_infeasible_init_point(fo, A, b,
const,
x_0, nu_0, tol,
tol_backtracking, x_ast, p_ast, maxiter,
mu, tol_outer_iter = tol_outer_iter
)
In [36]:
[x,iter_barrier,t]
Out[36]:
In [37]:
compute_error(x_ast,x)
Out[37]:
In [38]:
x1 = cp.Variable()
x2 = cp.Variable()
In [39]:
# Create two constraints.
constraints = [6-x1-x2 <= 0,-18+x1+2*x2<=0,x1>=0,x2>=0]
# Form objective.
obj = cp.Minimize(2*x1+5*x2)
In [40]:
# Form and solve problem.
prob = cp.Problem(obj, constraints)
prob.solve() # Returns the optimal value.
Out[40]:
In [41]:
print("status:", prob.status)
print("optimal value", prob.value)
print("optimal var", x1.value, x2.value)
Referencias: