58 results about "Arbitrary-precision arithmetic" patented technology

The invention discloses a heterogeneous Chiplet hybrid precision Transformer acceleration system based on an MRAM-SRAM (Magnetic Random Access Memory-Static Random Access Memory), which adopts a Chiplet block integration technology, gives full play to the modularization advantages of high yield and low cost of the Chiplet block integration technology, and internally integrates two types of heterogeneous calculation cores of MRAM calculation Chiplet and SRAM-NPU calculation Chiplet and two types of scheduling management Chiplet and data interaction Chiplet control units. In a Transform calculation process, an MRAM (Magnetic Random Access Memory) calculates a Chiplet to carry out approximate token relationship prediction in a self-attention mechanism, an SRAM-NPU (Synchronous Random Access Memory-Network Processing Unit) calculates the Chiplet to realize accurate attention calculation by utilizing a high-speed read-write characteristic of an SRAM (Static Random Access Memory) and a flexible calculation capability of the NPU, meanwhile, a scheduling management Chiplet coordinates calculation task allocation, and a data interaction Chiplet ensures efficient data flow; according to the acceleration system, optimal combination of different process node devices is realized through a Chiplet heterogeneous integration technology, and three dimensions of calculation efficiency, energy consumption and precision are balanced by combining the high-density low-power-consumption storage and calculation integrated advantage of an MRAM and high-precision calculation of an SRAM-NPU.

The invention relates to the technical field of digital twinning, in particular to a high-performance digital twinning system rendering method based on dual-grid and neural network mapping, and the method comprises the steps: constructing a high-precision calculation grid and a low-precision rendering grid, and building a space mapping relation between the two grids; executing multi-working-condition numerical simulation to obtain a high-dimensional simulation result snapshot matrix; the dimension reduction projection operator is used for reducing the dimension of a high-dimensional simulation result of any working condition into an r-dimensional feature vector; generating low-dimensional reference field data by using the space mapping relation; training the neural network until a mapping network parameter is obtained through convergence; obtaining current working condition parameters of the equipment in real time, and obtaining a current r-dimensional feature vector by using the dimensionality reduction projection operator; and quickly generating physical field prediction distribution on the low-precision rendering grid and performing pseudo-color rendering. According to the scheme, the problem that a traditional method is low in rendering efficiency is solved, and the method has the advantages of reducing computing resource consumption and keeping visualization precision.

The invention relates to an MVP formula-based edge element magnetostatic high-precision calculation method for smooth edges, and the method comprises the steps: dividing grid data into a plurality of polyhedral smooth regions with corresponding edges as centers according to the edges in the grid data; the polyhedral smooth domain is formed by surrounding corresponding edges by nodes, body centers and face centers of tetrahedron units in the grid data; calculating the vector shape function gradient of the corresponding edge based on the vector shape function of each tetrahedron unit in the polyhedral smooth domain corresponding to the edge; constructing a static magnetic field equation based on an MVP formula and coulomb specifications, and rewriting the static magnetic field equation according to the vector shape function gradient of each edge and a Galerkin method to obtain a system equation of a static magnetic field; further rewriting the system equation of the static magnetic field into a Newton iteration linear equation; and iteratively solving the Newton iteration linear equation by adopting a Newton iteration method until convergence to obtain the magnetic vector position in the static magnetic field. According to the method, the requirements of high calculation precision and stability of the engineering static magnetic field problem are met.

The invention relates to the technical field of hydraulic engineering, and discloses a reservoir capacity curve high-precision calculation method and system based on a three-dimensional model, and the method comprises the steps: carrying out the self-adaptive sampling and optimization of original terrain data, constructing a three-dimensional terrain surface model, carrying out the self-adaptive discretization grid division, and carrying out the volume analysis error correction. And high-precision reservoir capacity estimation is realized. A data integrity detection and repair mechanism is introduced, the calculation efficiency is improved by combining a parallel processing technology, and the result reliability is ensured by adopting cross validation. Meanwhile, a three-dimensional visual model is generated, layered display of topographic features is realized, a topographic change adaptability analysis function is realized, and a capacity calculation result can be dynamically updated. According to the method, the accuracy and efficiency of reservoir capacity calculation are remarkably improved, and reliable data support and a visual visualization tool are provided for hydraulic engineering planning and management.

The invention provides a comprehensive pipe gallery safety state evaluation method and system, and relates to the technical field of comprehensive pipe gallery safety state evaluation methods, and the method comprises the steps: employing a steady-state heat release method to calculate the simulation temperature of each three-dimensional grid node; a finite difference method is adopted to calculate a temperature field distribution coefficient of each three-dimensional grid node, and a temperature distribution coefficient of a three-dimensional grid center point is obtained through an arithmetic average method; a coupling damage method is adopted to construct a coupling coefficient matrix of any grid position, equivalent damage energy of any grid position is calculated through an inner product method, and the equivalent damage energy is mapped into a safety coefficient of any grid center point through normalization. By adopting three-dimensional grid discretization, complex geometric characteristics of the pipe gallery structure can be more accurately captured, the position of a potential fire source is clearly marked, the hidden danger recognition efficiency is improved through targeted monitoring, high-precision calculation of a simulated temperature field under the complex boundary condition is achieved, equivalent damage energy is mapped into a dimensionless safety coefficient, and cross-regional risk comparability is achieved.

The present invention discloses a mathematical model of a dual three-phase permanent magnet synchronous motor considering the electromagnetic-temperature coupling effect. Based on the finite element calculation results of the electromagnetic field and combined with mathematical methods such as multi-dimensional interpolation and inversion, a high-precision mathematical model of the electromagnetic module of the dual three-phase permanent magnet synchronous motor is established, which can fully consider the magnetic saturation, harmonic magnetic field and rotor position influence in the actual operation of the motor. On this basis, based on the finite element calculation results of the temperature field and using the linear system identification method, a fast and high-precision calculation matrix from various losses to the temperature rise of multiple motor components is established to realize the fast and high-precision calculation of the motor temperature rise, achieving the coupling from electromagnetic to temperature; further, the influence of temperature on the electromagnetic characteristics of the motor is analyzed and calculated, and it is fed back to the electromagnetic module of the motor to realize the coupling from temperature to electromagnetic. The model of the present invention has better fidelity and can provide a more accurate controlled object for high-quality motor control strategies.

The invention belongs to the technical field of in-memory computing, and discloses a multi-precision in-memory computing circuit with balanced output bit width, a memory and a chip, which are characterized in that through the collaborative change of a broadcast mode of an input side and an accumulation mode of an output side, under different computing precisions, the multi-precision in-memory computing circuit with balanced output bit width is obtained; the external interface bandwidth (the number of ADCs / the number of adders) of the physical array is kept constant. According to the invention, the contradiction between high bandwidth demand in a low-precision mode and low bandwidth demand in a high-precision mode in a multi-precision in-memory computing architecture is solved, the idle waste of the output bandwidth during high-precision computing is avoided, the relatively balanced bit width of an output interface is realized, and the area cost of a chip can be greatly reduced.

The invention discloses a method for optimizing calculation precision of floating-point instruction conversion based on extended control words. By extending an ARM floating-point calculation control word, pre-establishing a rounding mode mapping table, an exception type mapping table and a virtual exception instruction table, the conversion of an x86 floating-point control word loading instruction is realized based on the rounding mode mapping table and the exception type mapping table. For an x86 floating-point calculation instruction, the ARM instruction sequence to which the instruction should be converted is preliminarily determined according to whether the instruction is a high-precision calculation-related instruction. Then, according to whether the instruction exists in a virtual exception instruction table, it is determined whether to add an instruction group for simulating exception triggering to the instruction, so as to complete the conversion of the x86 floating-point calculation instruction, realize the conversion and execution of the x86 floating-point calculation-related instruction to the ARM system, reduce the precision loss of program cross-architecture execution, and improve the reliability of program cross-architecture execution.

The application discloses a lightweight method of river surface flow velocity field depth reconstruction model, relates to the computer science and fluid mechanics cross technical field, collects river surface PTV and PIV data, and constructs a flow velocity field matrix of the target river surface; based on the matrix block calculation strategy and the flow velocity field matrix, the target river surface is divided into a plurality of sub-matrix regions; based on the pattern recognition technology and the flow velocity gradient characteristics, each sub-matrix region is distinguished. The application divides the global flow velocity field into a plurality of sub-matrix regions through the matrix block strategy, adaptively determines the optimal block size according to the GPU parallel thread number, memory bandwidth and other constraints, maximizes the hardware utilization rate, and in the calculation process, adopts a differentiated resource allocation mechanism, concentrates high-precision calculation in a high-error clustering area, and adopts a simplified model in a low-error clustering area, significantly reduces the overall calculation load, and realizes efficient balance between calculation efficiency improvement and resource dynamic scheduling.

The invention provides a slow action generation method based on motion focus area analysis. The method comprises the following steps: determining a video of a slow action to be generated; detecting a motion focus area in the video, and grading the motion focus area; and for different grades of motion focus areas, different strategies are adopted to carry out motion modeling and differential interpolation processing, and a slow motion frame sequence is generated. According to the invention, through identification and grading processing of the motion focus area, high-precision calculation of the core motion object, simplified calculation of the secondary area and simplest processing of the static background, compared with a traditional scheme of uniform and dense calculation of the whole frame, the method has the advantages of high calculation complexity and lower calculation amount, and guarantees the visual effect of the core area.

The invention discloses a lightweight method for a river surface flow velocity field depth reconstruction model, and relates to the technical field of crossing of computer science and fluid mechanics, and the method comprises the steps: collecting the PTV and PIV data of a river surface, and constructing a flow velocity field matrix of a target river surface; dividing the surface of the target river into a plurality of sub-matrix areas based on a matrix partitioning calculation strategy and the flow velocity field matrix; and based on a mode recognition technology and flow velocity gradient features, distinguishing each sub-matrix region. According to the method, a global flow velocity field is divided into a plurality of sub-matrix areas through a matrix partitioning strategy, the optimal partitioning size is adaptively determined according to constraints such as GPU parallel thread count and memory bandwidth, the hardware utilization rate is maximized, and in the calculation process, a differentiated resource allocation mechanism is adopted, high-precision calculation is concentrated in a high-error clustering area, and the calculation efficiency is improved. A simplified model is adopted in the low-error clustering area, the overall calculation load is remarkably reduced, and efficient balance between calculation efficiency improvement and resource dynamic scheduling is achieved.

The invention provides a multi-level virtual-real coordinate mapping and drift correction method, system and platform, a dynamic floating original point mechanism is established, the original point position of a virtual coordinate system is adjusted in real time, so that a scene around a user is always in a high-precision calculation range of a single-precision floating-point number, and the user experience is improved. The dynamic floating origin mechanism comprises a coordinate mapping engine, an origin decision module and a seamless switching manager. A multi-level real world-virtual-real space-virtual environment conversion link is directly opened, large-scale development of augmented reality XR contents of cities / parks / outdoor / indoor / temporary spaces, automatic driving, unmanned aerial vehicles, intelligent robot navigation and the like is promoted, and the problems of large-space rendering jitter, low development efficiency and the like in the prior art are solved. In order to achieve the purpose, the invention provides an underlying OS operating system, and provides a cross-multi-level real world coordinate system and virtual environment coordinate system mapping method, a drift correction system and a developer operating platform.

The invention discloses an intelligent earthwork computing system and method based on multi-source data fusion. The method comprises the following steps: acquiring a multi-angle image and high-precision coordinate data of a construction area through unmanned aerial vehicle oblique photography and RTK positioning technologies; combining fractal noise reduction and a wavelet transform algorithm to eliminate vegetation noise and extract topographic feature point cloud; a dynamic three-dimensional terrain model is constructed based on BIM platform fusion point cloud, GIS geographic information and design parameters, and excavation and filling amount deviation and geological disaster signal triggering model updating are detected in real time; optimizing an earthwork transportation path by adopting an improved ant colony algorithm, defining a dynamic weight cost function, and calculating a pavement paving area in combination with a Green's formula; and generating a visual report through the AR technology, superposing the difference between the design model and the construction state, interactively adjusting the scheme, and issuing an instruction to the unmanned transportation equipment. According to the method, high-precision calculation, multi-parameter dynamic optimization and construction whole-process closed-loop management of the earthwork volume under the complex terrain can be realized, and the method is suitable for real-time refined management and control of large-scale projects.

The invention discloses a multi-scale space target relative pose estimation method and system using parallel double encoders, and belongs to the technical field of space target pose estimation. The invention aims to solve the problems caused by distance change, illumination difference and insufficient multi-scale feature expression of a space target in a dynamic environment. According to the technical key points, a parallel double-encoder architecture is adopted, and local and global features of a space target are extracted; constructing a feature fusion module to fuse features extracted by the CNN and the Transform; multi-scale features are fused through FPN, and angular points are predicted; training and optimizing the network; and solving the pose by using a random consistency perspective n-point method. A plurality of non-coplanar control points are selected, other points are expressed as a linear combination of the control points, a linear equation set is constructed, and a rotation matrix and a translation vector between a target and a camera coordinate system are solved by using a singular value decomposition method. The deep learning network algorithm provides an efficient and accurate solution for space target pose estimation, the algorithm can deeply mine rich information of image data by adopting an encoder-decoder structure, and finally high-precision calculation of the six-degree-of-freedom pose of the target is achieved.

The application discloses a kind of quantization bit width dynamic selection neural network acceleration hardware architecture and method, hardware architecture includes global storage module, data scheduling module, local storage module, dynamic quantization prediction controller and computing unit array;The feature map in network is divided into neurons with block as unit, in each block, continue to be divided into group with space dimension as unit, define group as the minimum unit of executing dynamic quantization operation;Each block in all feature maps of target network is configured with a trainable threshold parameter, and according to the given basic quantization bit width and target bit width total quantity constraint, the upper and lower bounds of the sparsity selectable by each block are determined;Establish the training and inference model of dynamic quantization neural network, divide inference into high-precision and low-precision calculation two parts, whether the low-precision calculation is executed is judged by the result of high-precision calculation;On the premise of ensuring that the accuracy is not lost, the inference time in actual hardware is reduced as much as possible.

Embodiments detailed herein relate to arithmetic operations of float-point values. An exemplary processor includes decoding circuitry to decode an instruction, where the instruction specifies locations of a plurality of operands, values of which being in a floating-point format. The exemplary processor further includes execution circuitry to execute the decoded instruction, where the execution includes to: convert the values for each operand, each value being converted into a plurality of lower precision values, where an exponent is to be stored for each operand; perform arithmetic operations among lower precision values converted from values for the plurality of the operands; and generate a floating-point value by converting a resulting value from the arithmetic operations into the floating-point format and store the floating-point value.

The invention discloses a symmetric magnetic field calculation method, which comprises the following steps of: analyzing the symmetry of magnetic field source distribution, and selecting a reference source in a group of symmetric sources; calculating the magnetic field intensity of the reference source at any field point in the space; based on symmetry, mapping the field points to equivalent field points corresponding to the symmetric sources through a rotation matrix, and calculating magnetic fields of the symmetric sources at the equivalent field points; restoring the magnetic field at the equivalent field point to an original field point direction through an inverse matrix of the rotating matrix; and finally, superposing the magnetic field contributions of all symmetrical sources at the original field point to obtain the total magnetic field of the point. According to the method, the symmetry of a magnetic field source is fully utilized, a rotating group processing method in the group theory is adopted, accurate calculation is performed in combination with the Bios-Savart law, the symmetry and direction correctness of the magnetic field can be strictly ensured, and high-precision calculation of the extremely low-intensity magnetic field is realized; the problems of errors and direction errors existing in symmetrical field and weak field calculation of a traditional finite element method are solved.

To solve the problem that a highly precise calculation cannot be performed efficiently since a time increment required from local spatial resolution becomes much smaller than an essentially required time scale and hence many repeated computations are required for the entire analysis period and computation time becomes enormous in a variation analysis of an environment fluid including a local component expressing a fast component of time variation and an average component expressing a slow component of time variation.SOLUTION: Computation considering a change in an actual time scale is performed by amplifying an average component of the amount of variation in a next time step by a constant on the basis of the amount of variation of an environment fluid in a current time step so as to reduce calculation time by calculating a short period obtained by reducing the time scale of an actual phenomenon as an object to 1 / constant, and a time change in various calculated variation amounts is extended by a multiple of a constant, thus speedily and precisely calculating a time change in the amount of variation in an actual time scale.SELECTED DRAWING: Figure 2

The invention belongs to the technical field of hydraulic structure numerical analysis, and particularly relates to a two-dimensional temperature field adaptive grid merging method. The method is based on an independent coverage manifold method, an error estimator is constructed through the heat flow continuity degree between adjacent grids to replace an empirical formula, full-domain automatic error estimation is achieved, the accuracy of judging whether grid merging is carried out or not is enhanced, the method comprises the steps of error estimation, grid marking, grid shape judgment, grid merging and the like, and the method is suitable for large-scale popularization and application. And by merging the local grids, the calculation precision and the calculation efficiency are both considered, and local high-precision calculation is realized.

The invention relates to the technical field of artificial intelligence accelerators, in particular to a configurable mixing precision calculation array circuit based on an outer product. According to the method, the basic int4 multiplier is taken as a reference, the results of the basic multiplier are fused according to different configuration modes to support high-precision calculation, calculation of 12 different specifications can be supported at most, the resource utilization rate of the multiplier is increased, and the resource utilization rate of the multiplier reaches 100%; by expanding a 1-bit sign bit, the array supports multiplying a signed number by a signed number, multiplying an unsigned number by a signed number and multiplying an unsigned number by an unsigned number, and various calculation requirements of the neural network are flexibly adapted; a matrix transposition switch is built in, so that the overhead of external transposition is reduced; the bit width of output data is limited through functions such as saturation truncation, and storage overhead is reduced. According to the method, the calculation efficiency is improved, at most four kinds of calculation precision and 12 kinds of calculation modes are supported, and the resource utilization rate is improved by 2-4 times; and application scene extension, perfect signed / unsigned calculation support and dynamic reconfiguration during operation are supported.

The application discloses a single-precision hardware implementation double-precision simulation method based on floating-point decomposition and reconstruction, including a 32-bit single-precision server; firstly, a prediction model is acquired to train floating-point 64-bit double-precision scalar or tensor original data, and the implementation double-precision simulation method comprises the following steps: step S1, data preprocessing; step S2, computation graph reconstruction; two floating-point 32-bit single-precision data are rewritten into a double-channel computation form by using a traditional computation flow; step S3, dynamic precision recovery; step S4, adaptive quantization control; an error feedback mechanism is introduced, when the error exceeds a threshold value, a dynamic decomposition result of computation is triggered, and the decomposition result is iterated to two floating-point 32-bit single-precision data, the application meets the high-precision computation demand of the single-precision server; the application is completely based on a single-precision instruction set and does not need special hardware support, and is suitable for an NPU cluster; compared with original double-precision storage, the memory occupation is reduced by 25%.

This invention discloses a high-precision approximation calculation device for the softmax function, comprising a first module, a second module, and a third module. The first module includes a first zero-value prediction unit, a first softmax function calculation unit, and a first selector unit. The second module includes a second softmax function calculation unit. The third module includes a second zero-value prediction unit, a third softmax function calculation unit, and a second selector unit. The first and second zero-value prediction units are used to predict whether the output data of the device is zero based on their corresponding input data, so as to skip subsequent calculations. The first, second, and third softmax function calculation units are used to implement high-precision calculation of the softmax function, and the first and second selector units are used to select the final output data. This invention can greatly reduce the power consumption and area overhead of the computing device, and improve the energy efficiency of the device.

The invention discloses a high-precision calculation method and device for a normalized arc tangent function. The method comprises the steps that S1, input floating point data is decomposed and preprocessed; the symbol, the order code and the mantissa of the input floating point data are obtained through data conversion, and the order code is compared with the single-precision offset code to obtain a shift code; obtaining input data for second-order polynomial calculation by using the shift code; s2, simplifying a calculation interval by using the characteristics of an odd function; s3, layering and segmenting the effective interval and selecting a polynomial coefficient; s4, carrying out second-order polynomial and fifth-order Taylor dual-path parallel calculation; step S5: result formatting; and S6, dynamically selecting and outputting. The device is used for executing the method. The invention has the advantages of high precision, low time delay, small area, low cost and the like.

The invention belongs to the field of numerical calculation of structural dynamics, and discloses a high-precision calculation method, application and equipment for structural dynamic response considering a load derivative. According to the equation, loads and load derivatives at the current moment and the next moment are considered, a load matrix at the current moment and a load matrix at the next moment are constructed respectively, and the load action part of a traditional discrete state equation is improved. The algorithm provided by the invention gives consideration to the requirements of rapid calculation and high precision, can be widely applied to the fields of structural vibration analysis, vibration control, health monitoring, response prediction and the like, and has a good application prospect.

The invention discloses an automatic CFD solving method for an internal flow field of a bimodal scramjet engine, and the method comprises the steps: constructing an AnsysFluent automatic solving environment through employing a Python script, loading a grid file of a bimodal scramjet engine model, configuring a physical model and boundary conditions of the bimodal scramjet engine through employing the Python script, and carrying out the automatic solving of the internal flow field of the bimodal scramjet engine. A Python script is utilized to complete setting of a numerical discrete format and a solving method and construction of a flow field initial value of a computational domain; taking the flow field initial value as an iteration starting point, and executing dynamic climbing calculation under first-order precision by utilizing a Python script to obtain stable initial flow field distribution; the stable initial flow field distribution is used as an initial calculation value, and a Python script is used for switching second-order high-precision calculation and carrying out convergence judgment. According to the method, the problems of tedious manual intervention, difficulty in shock wave capture and subjective convergence judgment in high-Mach-number strong compressible flow calculation are effectively solved.

The present invention discloses a high-precision estimation method for the angle of arrival, which successively includes the following steps: A: First, establish an off-grid DOA bilinear system model, and then factorize the off-grid DOA bilinear system model using the Bayesian formula; B: Estimate the sparse incident wave signals S ; C: Estimate the off-grid error parameter vector β ; D: According to the obtained sparse incident wave signals S and the off-grid error parameter vector β , perform high-precision calculation of the angle of arrival through the off-grid DOA bilinear system model obtained in step A. The present invention can solve the contradiction between complexity and robustness in off-grid angle of arrival estimation, effectively improve the calculation accuracy of the angle of arrival, and has obvious performance advantages in scenarios with fewer snapshots such as airborne radars.

The present disclosure relates to an information processing apparatus, an information processing method, and a computer-readable storage medium. The information processing apparatus according to the present disclosure includes processing circuitry configured to: obtain a sound field in a virtual space, by using a sound field synthesis model, based on sound source information about a sound source in a real space and virtual space information indicating an object present in a virtual space. The sound field synthesis model is pre-obtained through machine learning by using a sound field in a virtual space obtained from high-precision calculation as learning data.

The application provides a construction safety visual calculation and resource scheduling method, system and electronic equipment, which comprises the following steps: deploying and initializing the calculation and resource scheduling system; the edge end collects video streams in real time and performs preprocessing, light-weight model inference and local alarm; the edge end performs three-level screening on the video data, intelligent compression of the region of interest (ROI) and end-to-end encryption before uploading to the cloud; a multi-objective optimization decision model is constructed based on FAHP, and the allocation strategy of the task in the edge end and the cloud is dynamically adjusted according to the real-time collected network state parameters, computing power state parameters and event characteristic parameters; the cloud adopts a high-precision model for inference, links a knowledge base to generate a rectification scheme, and realizes full-process closed-loop management; an incremental data set is constructed to fine-tune the light-weight model and the high-precision model. The scheme provided by the application can balance the contradiction between the real-time alarm demand of the construction site and the high-precision calculation demand of the complex model.

The invention discloses a training and reasoning system, a training method and a reasoning method, and relates to the technical field of artificial intelligence, the training and reasoning system comprises a resistive random access memory storage layer, a dynamic random access memory storage layer and a data scheduling module; the resistive random access memory storage layer and the dynamic random access memory storage layer are heterogeneously integrated; the resistive random access memory storage layer is used for storing the weight after training stabilization and executing multiplication and addition operation; the dynamic random access memory storage layer is used for storing intermediate results obtained in training and reasoning processes; and the data scheduling module is used for storing the weight to be updated in the storage layer of the dynamic random access memory, storing the weight after training stabilization in the storage layer of the resistive random access memory, and storing an intermediate result obtained based on multiplication and addition operation in the training and reasoning process in the storage layer of the dynamic random access memory. The system can meet the requirements of high bandwidth, low power consumption, high capacity and high-precision calculation at the same time.

The application relates to a high-precision calculation method of an edge-smoothed edge element static magnetic field based on an MVP formula, which divides grid data into a plurality of polyhedral smoothing domains with corresponding edges as centers according to each edge in the grid data, and the polyhedral smoothing domain is formed by nodes, body centers and face centers of tetrahedral units in the grid data around the corresponding edge; the vector function gradient of the corresponding edge is calculated based on the vector function of each tetrahedral unit in the polyhedral smoothing domain corresponding to the edge; the static magnetic field equation is constructed based on the MVP formula and the Coulomb norm, and the static magnetic field equation is rewritten according to the vector function gradient of each edge and the Galerkin method, so that the system equation of the static magnetic field is obtained; the system equation of the static magnetic field is further rewritten into Newton iteration linear equations; the Newton iteration method is used to iteratively solve the Newton iteration linear equations until convergence, so that the magnetic vector potential in the static magnetic field is obtained. The method realizes the high calculation precision and stability requirement of an engineering static magnetic field problem.