Solving Equation Patents (Class 708/446)
  • Patent number: 7337265
    Abstract: Disclosed herein is an organization of cache memory for hardware acceleration of the FDTD method. The organization of cache memory for hardware acceleration of the FDTD method provides a substantial speedup to the finite-difference time-domain (FDTD) algorithm when implemented in a piece of digital hardware. The organization of cache memory for hardware acceleration of the FDTD method utilizes a very high bandwidth dual-port on-chip memory in a particular way. By creating many small banks of internal memory and arranging them carefully, all data dependencies can be statically wired. This allows for a many-fold speedup over SRAM-based solutions and removes the burden of data dependence calculation that streaming SDRAM-based solutions must perform.
    Type: Grant
    Filed: March 24, 2004
    Date of Patent: February 26, 2008
    Assignee: EM Photonics, Inc.
    Inventors: John R. Humphrey, James P. Durbano, Fernando E. Ortiz, Dennis W. Prather
  • Publication number: 20080005217
    Abstract: The invention provides a method and apparatus for operations of addition, subtraction, multiplication, and division of fractions, decimals, percents, proportions and equations. It allows the students to find the transparent area models on their sheets with the correct parts shaded, (2) students can count the total number of parts and the total number of shaded parts in each model, (3) the relationship of the shaded parts to the total number in each model can be named, (4) by intersecting the vertical transparent area model for one whole with horizontal transparent area models for one whole to show examples of one whole being renamed several different ways.
    Type: Application
    Filed: May 23, 2006
    Publication date: January 3, 2008
    Inventor: Dorothy Strong
  • Patent number: 7296047
    Abstract: One embodiment of the present invention provides a system that solves an overdetermined system of interval linear equations. During operation, the system receives a representation of the overdetermined system of interval linear equations Ax=b, wherein A is a matrix with m rows corresponding to m equations, and n columns corresponding to n variables, and wherein x includes n variable components, b includes m scalar components, and m>n. Next, the system performs a Gaussian Elimination operation to transform Ax=b into the form [ T W ] ? x = [ u v ] , wherein T is a square upper triangular matrix of order n, u is a vector with n components, v is a vector with m?n components, and W is a matrix with m?n rows and n columns, wherein W is zero except in the last column, which is represented as a column vector z with m?n components. Next, the system performs an interval intersection operation based on the equations zixn=vi (i=1, . . . , m?n) and Tnnx=un to solve for xn.
    Type: Grant
    Filed: May 28, 2003
    Date of Patent: November 13, 2007
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 7236998
    Abstract: A method and system for solving a large system of dense linear equations using a system having a processing unit and one or more secondary processing units that can access a common memory for sharing data. A set of coefficients corresponding to a system of linear equations is received, and the coefficients, after being placed in matrix form, are divided into blocks and loaded into the common memory. Each of the processors is programmed to perform matrix operations on individual blocks to solve the linear equations. A table containing a list of the matrix operations is created in the common memory to keep track of the operations that have been performed and the operations that are still pending. SPUs determine whether tasks are pending, access the coefficients by accessing the common memory, perform the required tasks, and store the result back in the common memory for the result to be accessible by the PU and the other SPUs.
    Type: Grant
    Filed: September 25, 2003
    Date of Patent: June 26, 2007
    Assignee: International Business Machines Corporation
    Inventors: Mark Richard Nutter, VanDung Dang To
  • Patent number: 7222145
    Abstract: One embodiment of the present invention provides a system for solving a system of equations in fixed-point form. During operation, the system receives a representation of the equations in fixed-point form and stores the representation a computer memory. Next, the system reduces the dimension of the system of equations, when possible, by eliminating variables in the system of equations to produce a reduced system of equations. The system then performs interval intersections based on the Fixed Point theorem to reduce the size of a box containing solutions to the reduced system of equations. In a variation on this embodiment, the system additionally applies interval techniques to find solutions to the system of equations, when such solutions exist.
    Type: Grant
    Filed: August 2, 2002
    Date of Patent: May 22, 2007
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Ramon E. Moore
  • Patent number: 7219115
    Abstract: A method for obtaining a stable and accurate solution for an ill-conditioned system of normal equations associated with digital Weiner filter for a time invariant system and/or an autoregressive operator of an autoregressive model. A time invariant stochastic model, uses a Gram_Schmidt process of orthonormalisation to condition the coefficient matrix, a singular matrix associated with such a system of normal equations, to an identity matrix. The observed output of the digital Weiner filter and/or autoregressive operator is defined as a time advanced version of the input. The method has application in situations where digitized data at smaller sampling intervals are made available.
    Type: Grant
    Filed: March 27, 2003
    Date of Patent: May 15, 2007
    Assignee: Council of Scientific & Industrial Research
    Inventor: Joseph Bernard Selvaraj
  • Patent number: 7155420
    Abstract: In one embodiment, a system and method for solving linear programs includes a perceptron algorithm configured to move toward a solution to the linear program. A transform algorithm is configured to stretch portions of a vector space within which the linear program is defined. A decision module decides between continued application of the perceptron algorithm and application of the transform algorithm based on a rate at which the approximate solutions are approaching a satisfactory solution.
    Type: Grant
    Filed: April 30, 2003
    Date of Patent: December 26, 2006
    Assignee: Microsoft Corporation
    Inventors: John Dunagan, Santosh S. Vempala
  • Patent number: 7152023
    Abstract: An active set algorithm exploits a ‘hot start’ for the set of binding constraints at optimality along with efficient linear algebra to make rapid progress towards the solution. The linear algebra is designed to deal with degenerate constraints as the required factorizations are performed and as degeneracy emerges, and not via a mostly unnecessary pre-process step. Combined together, these novel approaches enable solution of the control problem in real-time.
    Type: Grant
    Filed: February 14, 2003
    Date of Patent: December 19, 2006
    Assignee: United Technologies Corporation
    Inventor: Indraneel Das
  • Patent number: 7133888
    Abstract: A method and apparatus for quantum computing. A computer-program source code, data, and unsubstantiated output variables are converted into a class of computable functions by a program compiler. The computable functions are encoded, and a continualization method is applied to the encoded functions to determine a first-order, time-dependent, differential equation. Variational calculus is employed to construct a Lagrangian whose minimum geodesic is the solution for the first-order, time-dependent, differential equation. The Lagrangian is converted into a quantum, canonical, Hamiltonian operator which is realized as an excitation field via an excitation generator. The excitation field is repeatedly applied to a quantum processor consisting of a lattice of polymer nodes to generate an intensity-versus-vibrational-frequency spectrum of the lattice nodes.
    Type: Grant
    Filed: February 17, 2004
    Date of Patent: November 7, 2006
    Assignee: Clearsight Systems Inc.
    Inventors: Wolf Kohn, Anil Nerode
  • Patent number: 7099851
    Abstract: One embodiment of the present invention provides a system that solves a global optimization problem specified by a function ƒ and a set of equality constraints q1(x)=0 (i=1, . . . , r), wherein ƒ is a scalar function of a vector x=(x1, x2, x3, . . . xn). During operation, the system receives a representation of the function ƒ and the set of equality constraints and stores the representation in a memory within a computer system. Next, the system and performs an interval global optimization process to compute guaranteed bounds on a globally minimum value of the function ƒ(x) subject to the set of equality constraints. Performing this interval global optimization process involves, applying term consistency to the set of equality constraints over a subbox X, and excluding portions of the subbox X that violate the set of equality constraints.
    Type: Grant
    Filed: December 13, 2001
    Date of Patent: August 29, 2006
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 7089275
    Abstract: One embodiment of the present invention provides a system that uses a block-partitioned technique to efficiently solve a system of linear equations. The system first receives a matrix that specifies the system of linear equations to be used in performing a time-based simulation. This matrix includes a static portion containing entries that remain fixed over multiple time steps in the time-based simulation, as well as a dynamic portion containing entries that change between time steps in the time-based simulation. Next, the system performs the time-based simulation, wherein performing the time-based simulation involves solving the system of linear equations for each time step in the time-based simulation. In solving the system of linear equations, the system factorizes the static portion of the matrix only once, and reuses the factorization of the static portion in performing an overall factorization of the matrix for each time step of the time-based simulation.
    Type: Grant
    Filed: January 29, 2003
    Date of Patent: August 8, 2006
    Assignee: SUN Microsystems, Inc.
    Inventor: Rajat P. Garg
  • Patent number: 7079145
    Abstract: A projected iterative descent method is used to resolve LCPs related to rigid body dynamics, such that animation of the rigid body dynamics on a display system occur in real-time.
    Type: Grant
    Filed: March 8, 2004
    Date of Patent: July 18, 2006
    Assignee: AGEIA Technologies, Inc.
    Inventors: Richard Tonge, Lihua Zhang, Dilip Sequeira
  • Patent number: 7076516
    Abstract: A method of identifying one or more regions of the domain of a function that do not contain solutions is described along with a related subdivision method. These methods may be employed in the context of branch and bound methods that use interval analysis to search for solutions of functions. The one or more regions of the function domain that do not contain solutions are identified using a cropping formula derived from one or more components (low order and high order) of a Taylor Form inclusion function. A Corner Taylor Form inclusion function is also described which might be used to identify the output range of a function.
    Type: Grant
    Filed: September 18, 2001
    Date of Patent: July 11, 2006
    Assignee: California Institute of Technology
    Inventors: Marcel Gavriliu, Alan H. Barr
  • Patent number: 7073175
    Abstract: A static schedule is selected from a set of static schedules for an application dependent on the state of the application. A scheduling system stores a set of pre-defined static schedules for each state of the application. A scheduling system learns the costs of predefined schedules for each state of the application on-line as the application executes. Upon the detection of a state change in the application during run-time, the scheduling system selects a new static schedule for the application. The new static schedule is determined based on schedule costs and exploration criteria.
    Type: Grant
    Filed: June 3, 2005
    Date of Patent: July 4, 2006
    Assignee: Hewlett-Packard Development Company, Inc.
    Inventors: James M. Rehg, Kathleen Knobe
  • Patent number: 7065542
    Abstract: A calculating apparatus comprises a memory and a calculating section. The memory is stored with table data for referencing a value of a variable x by using a value of a variable y and a value of a solution z. The calculating section inputs x_in and y_in, reads two values y_a and y_b corresponding to y_in from the memory, reads two values x_aa and x_ab corresponding to x_in from the memory in a line corresponding to y_a, reads two values x_ba and x_bb corresponding to x_in from the memory in a line corresponding to y_b, and calculates a solution z corresponding to x_in and y_in in an interpolation by using x_in and y_in and the read values y_a, y_b, x_aa, x_ab, x_ba, and x_bb.
    Type: Grant
    Filed: June 11, 2002
    Date of Patent: June 20, 2006
    Assignee: Oki Electric Industry Co., Ltd.
    Inventor: Kiwamu Yoda
  • Patent number: 7065545
    Abstract: A computer method of vector operations for calculating the inverse of a general square matrix and for solving linear equations systems. The invention comprises a new method of factorization and executing multiply-add operations useful for effecting dot-product operations of one-dimensional vectors. This new method reduces the computation time over computers programmed to use prior art methods.
    Type: Grant
    Filed: May 7, 2002
    Date of Patent: June 20, 2006
    Inventor: Raul Gerardo Quintero-de-la-Garza
  • Patent number: 7062524
    Abstract: A system that solves a global inequality constrained optimization problem specified by a function ƒ and a set of inequality constraints pi(x)?0(i=1, . . . , m), wherein ƒ and pi are scalar functions of a vector x=(x1, x2, x3, . . . xn). The system performs an interval inequality constrained global optimization process to compute guaranteed bounds on a globally minimum value of the function ƒ(x) subject to the set of inequality constraints. The system applies term consistency and box consistency to a set of relations associated with the global inequality constrained optimization problem over a subbox X, and excludes any portion of the subbox X that violates the set of relations. The system also performs an interval Newton step on the subbox X to produce a resulting subbox Y. The system integrates the sub-parts of the process with branch tests designed to increase the overall speed of the process.
    Type: Grant
    Filed: January 8, 2002
    Date of Patent: June 13, 2006
    Assignee: Sun Microsystems, Inc
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 7028065
    Abstract: One embodiment of the present invention provides a system that solves a global optimization problem specified by a function ƒ and a set of inequality constraints pi(x)?0 (i=1, . . . , m), wherein ƒ and pi are scalar functions of a vector x=(x1, x2, x3, . . . xn). The system operates by receiving a representation of the function ƒ and the set of inequality constraints, and then storing the representation in a memory within the computer system. Next, the system performs an interval inequality constrained global optimization process to compute guaranteed bounds on the minimum value of the function ƒ(x) subject to the set of inequality constraints. While performing the interval global optimization process, the system applies term consistency at various places in the process over a subbox X, and excludes any portion of the subbox X that violates term consistency.
    Type: Grant
    Filed: December 13, 2001
    Date of Patent: April 11, 2006
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6993548
    Abstract: One embodiment of the present invention provides a system that bounds the solution set of a system of nonlinear equations specified by the set of linear equations Ax=b, wherein A is an interval matrix and b is an interval vector. During operation, the system preconditions the set of linear equations Ax=b by multiplying through by a matrix B to produce a preconditioned set of linear equations M0x=r, wherein M0=BA and r=Bb. Next, the system widens the matrix M0 to produce a widened matrix, M, wherein the midpoints of the elements of M form the identity matrix. Finally, the system uses M and r to compute the hull h of the system Mx=r, which bounds the solution set of the system M0x=r.
    Type: Grant
    Filed: January 31, 2002
    Date of Patent: January 31, 2006
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6976048
    Abstract: One embodiment of the present invention provides a system that solves a problem involving an interval parameter p through an interval solution process. During operation, the system receives a representation of the problem, wherein the problem includes a number of variables x1, x2, x3, . . . xn and at least one interval parameter p. The system stores the representation in a computer memory, and then performs the interval solution process on the problem. During this interval solution process, the system splits the problem into sub-problems by splitting the interval parameter p into subintervals, and creating separate sub-problems for each subinterval. The system then performs the interval solution process on the sub-problems. By splitting the interval parameter p, the system can achieve a tighter bound on the solution set of the problem. The decision to split on any parameter p is made in exactly the same way it would be made if p were a variable of the problem.
    Type: Grant
    Filed: March 11, 2002
    Date of Patent: December 13, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6961743
    Abstract: One embodiment of the present invention provides a system that solves a global optimization problem specified by a function ƒ and a set of equality constraints qi(x)=0 (i=1, . . . , r), wherein ƒ is a scalar function of a vector x=(x1, x2, x3, . . . xn). During operation, the system receives a representation of the function ƒ and the set of equality constraints and stores the representation in a memory. Next, the system performs an interval equality constrained global optimization process to compute guaranteed bounds on a globally minimum value of the function ƒ(x) subject to the set of equality constraints. During this process, the system applies term consistency to a set of relations associated with the interval equality constrained global optimization problem over a subbox X, and excludes any portion of the subbox X that violates the set of relations. It also applies box consistency to the set of relations, and excludes any portion of the subbox X that violates the set of relations.
    Type: Grant
    Filed: January 8, 2002
    Date of Patent: November 1, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6957341
    Abstract: A method and system for secure computational outsourcing and disguise. According to an embodiment, a first set of actual arguments and a second set of actual arguments for an outsourced computation are determined. A first group of disguised arguments corresponding to the first set of actual arguments is prepared with a first computer. A second group of disguised arguments corresponding to the second set of actual arguments is prepared with a second computer. The first and second groups of disguised arguments are output from the first and second computers, respectively, for performance of the outsourced computation. A third computer performs the outsourced computation and returns a disguised result to the first and/or second computers. The first and/or second computers then unveil the actual result from the disguised result.
    Type: Grant
    Filed: March 18, 2002
    Date of Patent: October 18, 2005
    Assignee: Purdue Research Foundation
    Inventors: John R. Rice, Mikhail J. Atallah
  • Patent number: 6950844
    Abstract: One embodiment of the present invention provides a system that performs a procedure to solve a system of linear inequalities. During operation, the system receives a representation of the system of linear inequalities Ax?b, wherein Ax?b can be a linearized form of a system of nonlinear equations. Within this representation, A is an interval matrix with m rows corresponding to m inequalities, and with n columns corresponding to n variables, the vector x includes n variable components, and the vector b includes m scalar interval components. The system solves the system of linear inequalities Ax?b by performing a Gaussian elimination process using only positive multipliers so as not to change the sense of any inequality.
    Type: Grant
    Filed: March 11, 2002
    Date of Patent: September 27, 2005
    Assignee: Sun Microsystems, Inc
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6922713
    Abstract: One embodiment of the present invention provides a system that receives a representation of the function ƒ and stores the representation in a memory. Next, the system performs an interval global optimization process to compute guaranteed bounds on a globally minimum value of the function ƒ(x) over a subbox X. This interval global optimization process applies term consistency to a set of relations associated with the function ƒ over the subbox X, and excludes any portion of the subbox X that violates any member of the set of relations. It also applies box consistency to the set of relations associated with the function ƒ over the subbox X, and excludes any portion of the subbox X that violates the set of relations.
    Type: Grant
    Filed: January 8, 2002
    Date of Patent: July 26, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6920472
    Abstract: One embodiment of the present invention provides a system for finding the roots of a system of nonlinear equations within an interval vector X=(X1, . . . , Xn), wherein the system of non-linear equations is specified by a vector function ƒ=(ƒ1, . . . , ƒn). The system operates by receiving a representation of the interval vector X (which is also called a box), wherein for each dimension, i, the representation of Xi includes a first floating-point number, ?i, representing the left endpoint of Xi, and a second floating-point number, bi, representing the right endpoint of Xi. Next, the system performs an interval Newton step on X to produce a resulting interval vector, X?, wherein the point of expansion of the interval Newton step is a point, x, within the interval X, and wherein performing the interval Newton step involves evaluating ƒ(x) to produce an interval result ƒ1(x).
    Type: Grant
    Filed: September 13, 2001
    Date of Patent: July 19, 2005
    Assignee: Sun Microsystems, Inc
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6920473
    Abstract: A multiplicand is multiplied by a multiplier using a modulus. The multiplicand, the multiplier and the modulus are polynomials of variable. A multiplication look-ahead method to obtain a multiplication shift value is carried out. An intermediate result polynomial is shifted to the left by the number of digits of the multiplication shift value. A reduction shift value equalling the difference of the degree of the shifted intermediate result polynomial and the degree of the modulus polynomial is obtained in a reduction look-ahead method. The modulus polynomial is then shifted by a number of digits equalling the reduction shift value. In a three-operands addition, the shifted polynomial and the multiplicand are summed and the shifted modulus polynomial is subtracted. The modular multiplication are iteratively executed and processed progressively until all the powers of the multiplier polynomial have been processed.
    Type: Grant
    Filed: July 21, 2003
    Date of Patent: July 19, 2005
    Assignee: Infineon Technologies AG
    Inventors: Astrid Elbe, Norbert Janssen, Holger Sedlak, Jean-Pierre Seifert
  • Patent number: 6915321
    Abstract: One embodiment of the present invention provides a computer-based system for solving a system of nonlinear equations specified by a vector function, f, wherein f(x)=0 represents ƒ1(x)=0, ƒ2(x)=0, ƒ3(x)=0 . . . , ƒn(x)=0, wherein x is a vector (x1, X2, X3, . . . xn). The system operates by receiving a representation of a subbox X=(X1, X2, . . . , Xn), wherein for each dimension, i, the representation of Xi, includes a first floating-point number, ai, representing the left endpoint of Xi, and a second floating-point number, bi, representing the right endpoint of Xi. The system stores the representation in a computer memory. Next, the system applies term consistency to the set of nonlinear equations, ƒ1(x)=0, ƒ2(x)=0, ƒ3(x)=0, . . . , ƒn,(x)=0, over X, and excludes portions of X that violate the set of nonlinear equations. The system also applies box consistency to the set of nonlinear equations over X, and excludes portions of X that violate the set of nonlinear equations.
    Type: Grant
    Filed: January 8, 2002
    Date of Patent: July 5, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6915320
    Abstract: One embodiment of the present invention provides a system for finding zeros of a function, ƒ, within an interval, X, using the interval version of Newton's method. The system operates by receiving a representation of the interval X. This representation including a first floating-point number, a, representing the left endpoint of X, and a second floating-point number, b, representing the right endpoint of X. Next, the system performs an interval Newton step on X, wherein the point of expansion is the midpoint, x, of the interval X. Note that performing the interval Newton step involves evaluating ƒ(x) to produce an interval result ƒI(x). If ƒI(x) contains zero, the system evaluates ƒ(a) to produce an interval result ƒI(a). It also evaluates ƒ(b) to produce an interval result ƒI(b).
    Type: Grant
    Filed: August 9, 2001
    Date of Patent: July 5, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6895422
    Abstract: One embodiment of the present invention provides a system for finding the roots of a polynomial or a quadratic equation with interval coefficients. The system operates by receiving a representation of a polynomial equation, which can be a quadratic equation of the form F(x)=Ax2+Bx+C=0, wherein A=[AL, AU], B=[BL, BU] and C=[CL, CU] are interval coefficients. Next, the system computes intervals containing roots of the functions F1(x), F2(x), F3(x) and F4(x), wherein F1(x)=ALx2+BLx+CL, F2(x)=AUx2+BUx+CU, F3(x)=ALx2+BUx+CL and F4(x)=AUx2+BLx+CU. The system then places the computed intervals into a list, L, and orders the computed intervals in L by their left endpoints, so that for a each entry, Si=[SiL, SiU], SiL?Si+1,L. Next, the system establishes interval roots for F(x) from the interval entries in list L.
    Type: Grant
    Filed: September 13, 2001
    Date of Patent: May 17, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Patent number: 6859817
    Abstract: A computer-based system for solving a system of nonlinear equations specified by a vector function, f, wherein f(x)=0 represents ƒ1(x)=0, ƒ2(x)=0, ƒ3(x)=0, . . . , ƒn(x)=0, wherein x is a vector (x1, x2, x3, . . . xn). The system operates by receiving a representation of an interval vector X=(X1, X2, . . . , Xn), wherein for each dimension, i, the representation of Xi includes a first floating-point number, ai, representing the left endpoint of Xi, and a second floating-point number, bi, representing the right endpoint of Xi. For each nonlinear equation ƒi(x)=0 in the system of equations f(x)=0, each individual component function ƒi(x) can be written in the form ƒi(x)=g(x?j)?h(x) or g(x?j)=h(x), where g can be analytically inverted so that an explicit expression for x?j can be obtained: x?j=g?1(h(x)). Next, the system substitutes the interval vector element Xj into the modified equation to produce the equation g(X?j)=h(X), and solves for X?j=g?1(h(X)).
    Type: Grant
    Filed: November 16, 2001
    Date of Patent: February 22, 2005
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Publication number: 20040267850
    Abstract: A method of measuring temperature of a component capable of emitting thermal radiation and reflecting background radiation comprising the steps of: providing a pyrometer for measuring the radiation from the component, characterised by coating a part of the component with a first emissivity coating and a part of the component with a different and second emissivity coating, each with known emissivities EH and EL respectively, recording a first radiation measurement from the first emissivity coating RH and a second radiation measurement from the second emissivity coating RL, then calculating the true radiation RBlade from the component from the equation 1 R Blade = ( R H E H - R L E H ⁢ ( 1 - E H 1 - E L ) ) ( 1 - E L E H ⁢ ( 1 - E H 1 - E L ) )
    Type: Application
    Filed: June 1, 2004
    Publication date: December 30, 2004
    Inventors: Colin Bird, Colin J. Parrish
  • Patent number: 6826585
    Abstract: In a simultaneous-linear-equations solving method of calculating the numerical solutions of simultaneous linear equations having a coefficient matrix, all the elements of coefficient matrix elements including zero elements and all the elements of right-side vector elements are stored into an array. Next, a non-zero-structure-specifying index table is created which indicates the row number of a terminal-end non-zero element in each column and the column number of a terminal-end non-zero element in each row within the array. Moreover, a decomposition processing is executed toward the elements existing within a range indicated by the created index table. Finally, a forward/backward substitution processing is executed toward the coefficient matrix elements subjected to the decomposition processing and the right-side vector elements stored into the array, thereby determining the numerical solutions.
    Type: Grant
    Filed: October 19, 2001
    Date of Patent: November 30, 2004
    Assignee: Hitachi Software Engineering Co., Ltd.
    Inventors: Saori Mitsunaga, Shinichi Tanaka, Hiroki Kawamura
  • Publication number: 20040236811
    Abstract: A method of computing shortest paths in a weighted graph having vertices and an adjacency matrix with memory resources and a processor including (a) selecting integer weights; (b) carrying out a series of incrementations, an incrementation including finding a set of vertices to which one may arrive from a given set of vertices; (c) carrying out a series of decrementations, a decrementation including finding a set of vertices from which one may go to arrive to a given set of vertices; (d) causing the incrementations and decrementations to be carried out in any order; (e) transforming vectors of increments/decrements in paths, the paths making up a set E1 of the shortest paths in term of number of arcs or using a given number of arcs, Na; (f) selecting n-uple of paths C of lowest cost among set of paths E1; (g) calculating Nb=Na+1; (h) computing iteratively, while Nb≦W(C) the following steps: i.
    Type: Application
    Filed: May 19, 2004
    Publication date: November 25, 2004
    Applicant: Kode, a corporation of France
    Inventor: Michel Koskas
  • Patent number: 6823352
    Abstract: One embodiment of the present invention provides a system for solving a nonlinear equation through interval arithmetic. During operation, the system receives a representation of the nonlinear equation ƒ(x)=0, as well as a representation of an initial interval, X, wherein this representation of X includes a first floating-point number, XL, for the left endpoint of X, and a second floating-point number, XU, for the right endpoint of X. Next, the system symbolically manipulates the nonlinear equation ƒ(x)=0 to solve for a first term, g1(x), thereby producing a modified equation g1(x)=h1(x), wherein the first term g1(x) can be analytically inverted to produce an inverse function g1−1(x). The system then plugs the initial interval X into the modified equation to produce the equation g1(X′)=h1(X), and solves for X′=g1−1[h1(X)].
    Type: Grant
    Filed: September 13, 2001
    Date of Patent: November 23, 2004
    Assignee: Sun Microsystems, Inc.
    Inventors: G. William Walster, Eldon R. Hansen
  • Publication number: 20040220987
    Abstract: In one embodiment, a system and method for solving linear programs includes a perceptron algorithm configured to move toward a solution to the linear program. A transform algorithm is configured to stretch portions of a vector space within which the linear program is defined. A decision module decides between continued application of the perceptron algorithm and application of the transform algorithm based on a rate at which the approximate solutions are approaching a satisfactory solution.
    Type: Application
    Filed: April 30, 2003
    Publication date: November 4, 2004
    Applicant: MICROSOFT CORPORATION
    Inventors: John Dunagan, Santosh S. Vempala
  • Publication number: 20040220988
    Abstract: An interleaver parameter generator circuit used to calculate and generate on an as needed basis interleaver parameters for interleaving blocks of information of varying lengths in accordance with a pseudorandom pattern defined by the 3GPP standard. The interleaver parameter generator circuit calculates and generates the defined interleaver parameters based on an input parameter that represents the length of the block of information to be interleaved. At least one of the defined parameters is calculated and generated using a decomposed form of its definition. The interleaver parameter generator circuit uses well known circuit blocks such as multipliers, subtractors, Compare-and-Select circuits and other circuits to calculate and generate the defined parameters.
    Type: Application
    Filed: May 1, 2003
    Publication date: November 4, 2004
    Inventors: Mark Patrick Barry, Benjamin John Widdup
  • Publication number: 20040210612
    Abstract: Elements of a vector sequence A·&phgr;m are sampled to be stored in a memory. In this sampling, a combination of spatial sampling and local sampling on the basis of physical phenomenon is employed. A residual minimization coefficient &agr;1m (wherein l=1, . . . , L) used for obtaining a corrected approximate value &phgr;m is approximately obtained by using elements of a vector sequence A·&phgr;k (wherein k=m−L+1, . . . , m−1) stored in the memory.
    Type: Application
    Filed: September 26, 2003
    Publication date: October 21, 2004
    Applicant: VINAS Co., Ltd.
    Inventors: Atsuhiro Tamura, Akihiro Ida
  • Patent number: 6799151
    Abstract: Matrix element calculation carried out efficiently without the overhead of communication between a host computer and processor elements even in parallel calculation utilizing a low-cost communication device and multiple processor elements having memories of a small capacity. In a method for calculating molecular orbitals, for example, all elements F(I, J) of a Fock matrix are calculated where an outermost loop is a loop associated with combinations (RT) of contracted shell R and contracted shell T which satisfy relationships R≦Nshell and T≦R. A second loop is a loop associated with contracted shell S, and a third loop is a loop associated with contracted shell U. Alternatively, the second loop is a loop associated with the contracted shell U, and the third loop is a loop associated with the contracted shell S. The value of S ranges from 1 to R, and the value of U ranges from 1 to R.
    Type: Grant
    Filed: April 7, 2000
    Date of Patent: September 28, 2004
    Assignees: Taisho Pharmaceutical Co., Ltd, Honda Motor Co., Ltd.
    Inventors: So Yamada, Shinjiro Inabata, Nobuaki Miyakawa, Hajime Takashima, Kunihiro Kitamura, Shigeru Obara
  • Patent number: 6799194
    Abstract: In a preconditioning process for an iteration method to solve simultaneous linear equations through multilevel block incomplete factorization of a coefficient matrix, a set of variable numbers of variables to be removed is determined at each level of the factorization such that a block matrix comprising coefficients of the variables can be diagonal dominant. The approximate inverse matrix of the block matrix is obtained in iterative computation, and non-zero elements of a coefficient matrix at a coarse level are reduced.
    Type: Grant
    Filed: June 26, 2001
    Date of Patent: September 28, 2004
    Assignees: Fujitsu Limited, Australian National University
    Inventors: Lutz Grosz, Makoto Nakanishi
  • Publication number: 20040181565
    Abstract: Diagonal elements of a triangular matrix are stored in memories 12 and 17, a computation using an output from each of shift stages REG1 to REG(N-1) of a shift register 11 and a diagonal element output from the memory 12 is performed, a computation result is input to the shift register 11, computation processing using a new register output from each of shift stages REG1 to REG(N-1) of the shift register 11 and the diagonal element output from the memory 12 is cyclically repeated, thereby solving a simultaneous linear equation.
    Type: Application
    Filed: February 2, 2004
    Publication date: September 16, 2004
    Inventor: Tetsuya Ikeda
  • Publication number: 20040167951
    Abstract: A·U+B(U)=f, wherein A is a linear differential operator, B is a nonlinear differential operator, and f is an inhomogeneous term (source term) in a nonlinear partial differential equation to be satisfied by a physical quantity U, is solved by successive approximation. In calculation, (f−A·Um−B(Um)) is given as a nonlinear residual rr of an approximate solution Um, wherein m is the number of repeating times, and the approximate solution Um is repeatedly corrected so as to reduce a norm of a nonlinear residual rm+1 employed in a subsequent step.
    Type: Application
    Filed: February 19, 2003
    Publication date: August 26, 2004
    Applicants: Kazuo Kikuchi, VINAS Co., Ltd.
    Inventors: Atsuhiro Tamura, Kazuo Kikuchi, Akihiro Ida
  • Publication number: 20040148324
    Abstract: One embodiment of the present invention provides a system that uses a block-partitioned technique to efficiently solve a system of linear equations. The system first receives a matrix that specifies the system of linear equations to be used in performing a time-based simulation. This matrix includes a static portion containing entries that remain fixed over multiple time steps in the time-based simulation, as well as a dynamic portion containing entries that change between time steps in the time-based simulation. Next, the system performs the time-based stimulation, wherein performing the time-based simulation involves solving the system of linear equations for each time step in the time-based simulation. In solving the system of linear equations, the system factorizes the static portion of the matrix only once, and reuses the factorization of the static portion in performing an overall factorization of the matrix for each time step of the time-based simulation.
    Type: Application
    Filed: January 29, 2003
    Publication date: July 29, 2004
    Inventor: Rajat P. Garg
  • Publication number: 20040143617
    Abstract: The present invention provides a method and system that produces a near-optimum schedule in linear time by providing an optimal resource ordering scheme. The present invention is embodied in a scheduling computer program.
    Type: Application
    Filed: October 23, 2003
    Publication date: July 22, 2004
    Inventors: Wolf Kohn, Vladimir M. Brayman
  • Publication number: 20040122882
    Abstract: A method for solving a system of N linear equations in N unknown variables.
    Type: Application
    Filed: October 15, 2003
    Publication date: June 24, 2004
    Inventors: Yuriy Zakharov, Timothy Conrad Tozer
  • Publication number: 20040034678
    Abstract: The poor scalability of existing superscalar processors has been of great concern to the computer engineering community. In particular, the critical-path delays of many components in existing implementations grow quadratically with the issue width and the window size. This patent presents a novel way to reimplement these components and reduce their critical-path delay growth. It then describes an entire processor microarchitecture, called the Ultrascalar processor, that has better critical-path delay growth than existing superscalars. Most of our scalable designs are based on a single circuit, a cyclic segmented parallel prefix (cspp). We observe that processor components typically operate on a wrap-around sequence of instructions, computing some associative property of that sequence. For example, to assign an ALU to the oldest requesting instruction, each instruction in the instruction sequence must be told whether any preceding instructions are requesting an ALU.
    Type: Application
    Filed: June 27, 2003
    Publication date: February 19, 2004
    Applicant: Yale University
    Inventors: Bradley C. Kuszmaul, Dana Sue Henry-Kuszmaul
  • Patent number: 6694343
    Abstract: A computer-based method and system comprising three data structures: partially ordered data structure (or simply ordered data structure), contiguous list v, and vector p, is used for solving a large sparse triangular system of linear equations which utilizes only the non-zero components of a matrix to solve large sparse triangular linear equations and generates explicitly only the non-zero entries of the solution. A list of the row indices of the known non-zero values of x which require further processing is stored in the ordered data structure. Actual non-zero values of x are stored in the contiguous list v and the corresponding pointers to the location of these values are stored in the vector p. The computer-based method manipulates these three matrices to find a solution to an upper or lower sparse triangular system of linear equations.
    Type: Grant
    Filed: February 8, 2001
    Date of Patent: February 17, 2004
    Assignee: International Business Machines Corporation
    Inventors: John Joseph Forrest, Nimrod Megiddo, John Anthony Tomlin
  • Publication number: 20040015531
    Abstract: One embodiment of the present invention provides a system that computes interval parameter bounds from fallible measurements. During operation, the system receives a set of measurements z1, . . . , zn, wherein an observation model describes each z1 as a function of a p-element vector parameter x=(xi, . . . , xp). Next, the system forms a system of nonlinear equations zi−h(x)=0 (i=1, . . . , n) based on the observation model. Finally, the system solves the system of nonlinear equations to determine interval parameter bounds on x.
    Type: Application
    Filed: July 15, 2003
    Publication date: January 22, 2004
    Inventors: G. William Walster, Eldon R. Hansen
  • Publication number: 20030212723
    Abstract: A computer method of vector operations for calculating the inverse of a general square matrix and for solving linear equations systems. The invention comprises a new method of factorization and executing multiply-add operations useful for effecting dot-product operations of one-dimensional vectors. This new method reduces the computation time over computers programmed to use prior art methods.
    Type: Application
    Filed: May 7, 2002
    Publication date: November 13, 2003
    Inventor: Raul Gerardo Quintero-de-la-Garza
  • Patent number: 6636828
    Abstract: The coefficient matrix, corresponding to the simultaneous linear equations to be solved, is divided into a plurality of row sets. The row sets as divided are processed in a parallel fashion, and entries specifying the nonzero elements contained in the first to nth row sets are added to the entry sets E1 to En. Moreover, in regard to each row set, fill-ins which take place at the time of eliminating the ith variable are obtained in a parallel fashion, and entries specifying the fill-ins are added to the entry sets E1 to En. The coefficient matrix is compressed using those entry sets E1 to En.
    Type: Grant
    Filed: May 10, 1999
    Date of Patent: October 21, 2003
    Assignee: NEC Electronics Corp.
    Inventor: Koutaro Hachiya
  • Publication number: 20030172099
    Abstract: One embodiment of the present invention provides a system that performs a procedure to solve a system of linear inequalities. During operation, the system receives a representation of the system of linear inequalities Ax≦b, wherein Ax≦b can be a linearized form of a system of nonlinear equations. Within this representation, A is an interval matrix with m rows corresponding to m inequalities, and with n columns corresponding to n variables, the vector x includes n variable components, and the vector b includes m scalar interval components. The system solves the system of linear inequalities Ax≦b by performing a Gaussian elimination process using only positive multipliers so as not to change the sense of any inequality.
    Type: Application
    Filed: March 11, 2002
    Publication date: September 11, 2003
    Inventors: G. William Walster, Eldon R. Hansen