神楽坂帕琪 mudream4869

## raytracing.py
import numpy as np
import matplotlib.pyplot as plt

w = 400
h = 300

def normalize(x):
    x /= np.linalg.norm(x)
    return x

## ipm.cpp
#include <Eigen/Core>
#include <Eigen/Sparse>
using namespace std;
using namespace Eigen;

void sinitialize(VectorXd &s,double mu,unsigned int sizex,VectorXd &x)
{
  for (unsigned int i=0;i<sizex;i++)
    if (x[i]!=0)s[i] = mu/x[i];
    else s[i]=0.1;

## CHANGELOG.md

      
        
          
            
              
              1 file
            
          
          
            
              
              0 forks
            
          
          
            
              
              0 comments
            
          
          
            
              
              0 stars
            
          
        
        
          
              
          
          
            
                mudream4869
                / CHANGELOG.md
            
            
              Created
              March 2, 2022 04:04
                — forked from juampynr/CHANGELOG.md
            
              
                Sample CHANGELOG
              
          
        
      
        
  
      
    Change Log

All notable changes to this project will be documented in this file.
The format is based on Keep a Changelog
and this project adheres to Semantic Versioning.
[Unreleased] - yyyy-mm-dd

Here we write upgrading notes for brands. It's a team effort to make them as
	import numpy as np
	import matplotlib.pyplot as plt

	w = 400
	h = 300

	def normalize(x):
	x /= np.linalg.norm(x)
	return x
	#include <Eigen/Core>
	#include <Eigen/Sparse>
	using namespace std;
	using namespace Eigen;

	void sinitialize(VectorXd &s,double mu,unsigned int sizex,VectorXd &x)
	{
	for (unsigned int i=0;i<sizex;i++)
	if (x[i]!=0)s[i] = mu/x[i];
	else s[i]=0.1;